{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# HiC link size empirical distributions for datasets in ALLHiC paper\n",
    "\n",
    "We plot the link size (x-axis) vs link density (y-axis) for each of the four datasets used in the ALLHiC paper. We then try to infer the coeffiencients in the power law model.\n",
    "\n",
    "$$Y=A \\times X^B$$\n",
    "\n",
    "Based on `allhic extract`, the coefficients for the four datasets are:\n",
    "- Ler0 (AAGCTT): $Y=0.00162 \\times X^{-0.8135}$\n",
    "- Rice (GATC): $Y=0.000347 \\times X^{-0.8042}$\n",
    "- AP85 Chr1A (AAGCTT): $Y=1.62e^{-5} \\times X^{-0.8164}$\n",
    "- Molokai Chr01A (AAGCTT): $Y=0.000325 \\times X^{-0.8844}$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Dependencies"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import matplotlib.pyplot as plt\n",
    "import seaborn as sns\n",
    "sns.set_style(\"darkgrid\")\n",
    "\n",
    "import pandas as pd"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Load data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>#Bin</th>\n",
       "      <th>BinStart</th>\n",
       "      <th>BinSize</th>\n",
       "      <th>NumLinks</th>\n",
       "      <th>TotalSize</th>\n",
       "      <th>LinkDensity</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>2048</td>\n",
       "      <td>90</td>\n",
       "      <td>45072</td>\n",
       "      <td>119653414</td>\n",
       "      <td>0.000004</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>2138</td>\n",
       "      <td>95</td>\n",
       "      <td>49670</td>\n",
       "      <td>119652784</td>\n",
       "      <td>0.000004</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2</td>\n",
       "      <td>2233</td>\n",
       "      <td>99</td>\n",
       "      <td>44248</td>\n",
       "      <td>119652119</td>\n",
       "      <td>0.000004</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>3</td>\n",
       "      <td>2332</td>\n",
       "      <td>103</td>\n",
       "      <td>44274</td>\n",
       "      <td>119651426</td>\n",
       "      <td>0.000004</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>4</td>\n",
       "      <td>2435</td>\n",
       "      <td>108</td>\n",
       "      <td>45417</td>\n",
       "      <td>119650705</td>\n",
       "      <td>0.000004</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   #Bin  BinStart  BinSize  NumLinks  TotalSize  LinkDensity\n",
       "0     0      2048       90     45072  119653414     0.000004\n",
       "1     1      2138       95     49670  119652784     0.000004\n",
       "2     2      2233       99     44248  119652119     0.000004\n",
       "3     3      2332      103     44274  119651426     0.000004\n",
       "4     4      2435      108     45417  119650705     0.000004"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def load_hic_data(filename):\n",
    "    \"\"\" Load tabular datasets on `distribution.txt`\n",
    "    File can be loaded directly into dataframe as first line contains header.\n",
    "    \"\"\"\n",
    "    df = pd.read_csv(filename, sep=\"\\t\")\n",
    "    xf = df[df[\"NumLinks\"] != 0]\n",
    "    return xf\n",
    "\n",
    "df = load_hic_data(\"ArabidopsisLer_RE_AAGCTT.distribution.txt\")\n",
    "df.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAESCAYAAAAIfCk9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3Xd4VFX6wPHvZNIbqfQS6qH33nuRKqAiIKur+7OvupZV11VXXVddC+5a0bWiIogo0qRJ71IklEMNoSW09J6Z+f1xhzEhJBlIJpNM3s/z5JFb5t43J+O8c8o9x2Sz2RBCCCGuxsvdAQghhKi8JEkIIYQoliQJIYQQxZIkIYQQoliSJIQQQhRLkoQQQohiebs7AFG5KKV8gHhgt9Z61HW8/jMgVmv9+lWO7QYGaq2Tr9j/GNBWa337dQVdfCz3AGFa61ecPP92YLLWekwZ79sDeAcIBs4A07XWZ4s5917AorWepZSKA3KALMAG+AFW4DGt9TKlVAxwFNh7lUt111rnOhFbFHAS+Exrfe+1HFdKeQMPA9MBHyAQ2AP8VWutC5w3BngMCLOfF2v/HU4qpb4DmtlP7WA/ZgGSgJ+BKfZjzYDzQIp9exJwH/CT1npNab+nKD+SJMSVJgK7ga5KqVZa6wPldWGtdcfyupaT9/ugIu8HoJTyBb4DpmitN9qTwP+AG65ybiPgdqBngd3TtNY7CpwzGfgUqGPflVXGcrwT+BGYqpT6m9b60jUc/wIIAIZorS/a45sKrFJKtdRap9u3nwHGaa2PKKVMwJPAL0qpNlrryQV+NxswSGt9ocA9XrEfWwO8o7X+rsD5LwAblFLdtdZZZSgDcQ0kSYgr3QvMwfjG+hBwj1JqIPA2kIHx7bgb8BrGh1sIYALu0lpvtF+jr/3DLRRYjvEtMt/+oRCN8e3wP8Aw4ByQaN+HUqo+8D4QY7/u51rrf9u/Ra8FlgE97Mce0FqvV0q1xPgg9rfv/1hr/Z5S6nkgSmv9gP3D+h4gF8gG7tZa73e2UJRS9TBqBw0xvh3P0Vq/bI9rPXDAHvMdQGqBsvgfMFMpFXn5g7WAp4AvtdZXfaLV/gHbGLjyg7y4GGcAf7nKodu01nuVUl7A3RjfyIOB/8P+oWx/fbHHlVJdgCFAY6115uXXaK2/VkpFYPyt04F/Av+ntT5iP25TSr2CUTv1w6gpXRetdYpSaqM9rrev9zri2kifhHBQSrUGegHzgM+BGUqpSPvhtsCtWuv2QGegLtBLa93afu6TBS5VH+MDpSNGk8KfrrjVfUALoDVGomhY4NhXwC9a63ZAH2C6UupyE0RDYK39m/STwLf25rHHMZohumB8Y+9v/8C7/HuZgZnASK11N2AW0Pcai+dL4BP7PboDQ5VSNxf4fV/UWrcAGmA01wBgbwI6D9QreDF7ApgELLriPl8ppfYopU7Zr9MSGFvgeIBSavcVP+/a7/WF1rrjVX4uN0+NxGgiWonxN3vA3oSEE8f7ARsKJogCv+M7Wusz9vdKDLDxiuM2rfVXWuvUIqV67ZZj1HZFBZGahCjoXmCR/RvvRaXUcYxvbZuBk1rrEwBa681KqWeAu5VSTYGBQFqB63yptc4AUErNBkZj1A4uGwp8bf8AzVVKfQW0V0oFYSSG4fb7pNj7OEYBW4AkrfXX9mNLlVIWoD2wAPhCKdUd4wPuz1prq1IK+7kWpdQ8YJNSajFG2/fXzhaKPa4BQIRS6kX77mCMJLgNyLeXERhfvK6sGZgw2t0LisToL4m7Yv80rfUOpVRj+++yW2t9rMDxYpubSqtJYPx9v7LX6hYCHwI3Ad/YzyvpuKng76WMwv3WvlkDo2Y5p0AZuMpxQLnw+uIKUpMQgOOD8DaMpqI4eydqHeABjOaV9ALnjgYW2zd/BD7A+BC5rOAHoheQd5VbFjw/v8C5pivO87Lfv+B5BY9ZtNaLgObAXKATsNfebOWgtZ6O8Y38CEYt5BucZ7bH1fvyt3OMpraX7cdztNaXY4vHqGUBjoEAkcDpK65pA0wFazxXxHsc4+/xuj35laqkmoS9/+MGYIr9b6sxviQ+Yo+zxOPAJqC3/fdBGy6XxVogSGudBByicB/L5XKYq5Tq4MzvUYo8iiZc4UKSJMRl04CLQF2tdYzWOgZogvGNOfqKc4dhNO+8D+wAJmB8kF42RSnlp5TyB/4ALL3i9UsxmrL87efcAqC1TsOoMdwPoJSqAcwAVthfF62UGmk/NhbjA2OvUupr4Bat9RyMpqxUoOnlmymlopRSJ4GLWuuZGB2r3ZwtGHszyRbs39KVUmEYTSrjr3L6ViBSKdXbvv1HYPOVI7rstbUkoFEJ992E0Vn8XnHJ5BrcjdFcVK/A37cL0Nkea4nHtdZbgV+AL+0joABHE2VHfv/g/gfwtlKqmf242V7r7AgcLOPvAEYfTXlcRzhJmpvEZfcCb2qtHd/StNbJSqn/8Pu3ycs+AL5RSu3FeA8tByYV+CA7jtGZG4LRFPT5Fa//EGOIYyxGYjpc4Ng04F2l1B2AL0az0GcYH6bZwG1KqVcxholOsDclvQh8rJS6G+PDagGwDhhk/z0uKKVewhiFk4VRI7myn+SykUqp9ALbyVrr+sBU4B377+wLfKO1/srece2gtc5TSk20nxtk//1mFHOv+Rj9AO8XcxyMzm1tj/dn7H0SVznvdq311fZfHnF1J0bCKhjrYaXUNxh9Or1LOP4IRk3iNowEvsReo4jA6Dd5D2ME1uWObBPG+8MHYzDBTmCw1vq6O60LGInRZyYqiEmmChdVgf3DOFZrHezuWMqLvd/hO6BrcSOcxO+UUqEYyaqr1jrb3fFUF9LcJISb2PsdPsdo6hGlex54WBJExZKahBBCiGJJTUIIIUSxJEkIIYQoVpUf3WS1Wm0WizSZmc0mpBxKJmVUOimj0nlKGfn4mC9QdHh7EVU+SVgsNpKTi8wUUO2EhQVKOZRCyqh0Ukal85Qyio4OOeHMedLcJIQQoliSJIQQQhRLkoQQQohiVfk+CSE8lcWST1LSefLzS11wrtwkJpqQZ6dKVtXKyNvbl/DwaMzm6/u4lyQhRCWVlHQef/9AgoJqYzJdOTmua5jNXlgs1gq5V1VVlcrIZrORkZFKUtJ5oqLqlP6Cq5DmJiEqqfz8XIKCQissQQjPYzKZCAoKLVNttFonCZvNxsWkFPLyrlymQIjKQRKEKKuyvoeqbXNTZnYOS9Zu4sTpBMJCgrltwkh8fXxKf6EQQlQj1bImYbVamf/zL5w4nQBAclo6h+NOlvIqIYSofqplksjJy+PcxaRC+/YfjXNPMEIIUYlVy+amAD8/hvTqyr7Dx0i4cAmA+DOJpGdmEhwY6ObohBCeLCsrizfeeAUfHx86derC8OGjHMcSEhJ4661XqVEjjAYNGjF48FC++OITMjLSeeml1wCIizvOvHnfkJKSTJcu3WnUKIaPP/6Axo2bMGTIcDp37lqu8VbLmgRAx1bNmTZuBA3r1HLsW7xmExeSUtwYlRCVS1zccW66aRxWqzHk02q18sgj97N06aJrvtaWLZu49daJ3HLLBL788rNrPqe4Yy+//A/GjBnGbbfdXOR6aWlpPPPME0ydOolp0yYTG/sbiYkJPPjg3UybNpnp029m7txvrvl3KWjBgu94/fVXHNuzZr3Hiy/+vdjz165dzcCBQ/jrX59hw4Z1hY6dPHmCXr368vTTzxEXd4x69erz1FPPFjonJqYxjz/+NC+88AoHD+7HZDIREBBIbm4uNWvWorxVy5pEQW1bNCH+bCIApxLO8/mCJTSsU4vhfbtTI8RjVsoU4rrExDSmYcMYNm1aT9++A/jww3dp2LARo0aNuabrWCwW3nzzVd56611q1qzFXXfNoG/f/jRu3MSpc0o6dsMNY5k06RZeeunZIvd9++3X6dGjNy+99Bp5eXlkZ2eTk5PDAw88glItyczM4I9/vI1u3XoUiiUp6RJ+fn4EBgY59p06dZL69RsUuceoUWOYOnUSaWn389tvu9m8eQMffPAJR48e4cMP3yl07lNPPcv58+do2rQZAF5ehb+nt2ih+OKLT1i9egUjRtxQbHlu2LCW2bM/Z+LEm+nQoROdOnXh0qWL/Pe/b/Hccy+V8te4NtW2JnFZyyaN6NmxbaFhYvFnE/nyx2XSmS0EcMstU1mwYD5r1qxi7949PPjgX675GgcO7KN+/QbUq1cfHx8fhg4dzoYNa50+p6RjHTt2JjQ0tMg9MzLS2bNnF2PGjAfAx8eHkJAQoqKiUKolAIGBQcTExHDhwrlCr921aydPPvkoOTk5ACxcuICZM/991d/N39+foUOHM2vWe8yc+TovvfQafn7+NG3ajNdem1noJzw8gujompw7Z9zPZiv8UN7ixT9x551385//fMDmzRuKLc++fQfwwQefsGLFUkeiCQkJJS+v/J/Or9Y1Ce9ftxP00vP0f/VNmowZxo69Bzh84hQ2m42c3DwWrt5Az45t6NO5vbtDFcJtunfvyTvvvMUHH7zLO+/Mwtu78MfGfffdRWZm0amz77//Ibp16wHA+fPnCjWFREfXZP/+2ELnl3SOM6+/0pkzpwkLC+Pll//BkSOHUKoVDz30GAEBAY5zzp49w6FDmtat2xZ67eDBQzl79jTPP/80gwYNZfHihbz11rvF3mv06PFMmzaZV155g3r16pcY14ABg3nrrdfYvHkDffr0Jzb2Nw4d0kyceBM9e/bik09msWLFMmrXrktKSjKzZr3HoUOaL7/8lNtuu4OdO3ewbt0v5Obm0rNnH9auXc3WrZtJT09n4sSiTW5lVX2TRF4eIY8+hPf+WMIH9cb7P+9TZ9LNnDl3gcVrNpGangHAlt37aBHTgOiIcDcHLKq7wNdeJsje9p3x2JNkPvF0oeNBzz5N4AdG80b68/8k674HCx0PfvTPBNjb8tNef5vsGXc4fe+2bdvTooUiKiqqyLH33vu41NdfbaqjKx/yKukcZ15/JYvFwqFDmocffoI2bdoyc+brzJ79GX/6070AZGZm8re/PcFDDz1KUFDRpuVp0/7Ac889xeuvv8K33/5AYAmDWj799CPCwsKxWCwlxgQQEBDA008/V2hf27bGF9EmTZo5Oqgve/zxwn/nzp27FumcHjBgcKn3vV7VtrnJZ8c2zIcOGhve3uR17Q5A3ZpR3DZ+JPVr/75g0469B90RohCVRlzccZo1U1c9dt99d3H77VOL/GzfvtVxTs2aNTl3LtGxff78OaKiCi+KVtI5zrz+StHRNYmOrkmbNkYtYdCgIRyy/z+fn5/PM888wfDhI4v9gN2zZxfHjh2lf/+BfPrprGLv8803s8nNzeWFF/7FvHlzSoypKqq2SSKvVx+SVq4nr0s3Mh57CmujGMcxfz9fBnTr5Ng+eOwEKWnpbohSiMrh+PFjNGnS9KrH3nvvYz777OsiP5ebmgBatmzNyZMnOXPmNHl5eaxcuZw+ffoXuk5J5zjz+itFRkZRs2Yt4uPjANixYxsxMU2w2Wz8618v0KhRY6ZMmX7V1x46dJBXX32JV155k6effo7U1FRmzXqvyHm//rqdJUsW8swzz9O5c1cyMtI5fFiXGFdVY6pKU95eTV6exVampQStVuPninZW32VL2LF2AyuatQWTiZqR4dxyw5BKO3WHpyyp6EpVrYwSEk5Qu3ajCr3n1WY4TUxM4IEH7mbevB/LdO3Nmzfw9ttvYrVaGD16HH/4w50APPbYn3nyyb8TFRVd7Dklvf65555m9+5fSU5OJiIikjvv/D/GjJkAwOHDmldeeYn8/Dzq1q3HU089x/Hjx7j//rto2rQZJpPxPfnuu++jV6++jnv99ttugoKCHaOQ8vPzWbLkJ8aNu9FRRgkJxlDaf/97JjExjQFYsuQndu36lb/97fkylVV5u9p7KTo65Feg1IcqJElchSk1hfA+3TAnJnCsSQu+vfkOUsIiaFSvNhOG9Mfb21yu9ysPVe0D0B2qWhlVliQhCquKZVSWJFFtm5tKEvDu25gTjXmd6iVfIsfPH4ATpxNYuHoDefkya6wQonqQJHEVmX9+lMz7H8JmNpPz6ut06t3dcez4qTPMW7qarOwcN0YohBAVo/oOgS1JUBAZz71I1m23Y23chF6A1Wpj6559AIRsWMui1BTGTh6Pv5+ve2MVQggXkiRRAqt9NIcJ6NulPYH+fuxeupzpsz8k19eXjbG7aP/CC3j7SDEK17DZbLLwkCiTsvY7S3PTNejcRnHnhmX45uUSnJFO+8ULWLFuU5VaFF1UHd7evmRkpMr7S1y3y2tce3tff4uHfAW+Rj73PUDW0cP4nUtkwY3TiD9xGsuaTYzs3xNvc+Ub9SSqrvDwaJKSzpOenlxh9zSZTJKUSlHVysjb25fw8JIfPCzx9eUYS7WQO3wUuZv7suODD4mPqAeAPh5PXn4+44b0w+dkfKEH84S4XmazN1FRdSr0nlVtmLA7VLcykuam6xEcQou//IUOLZs7dh07eYafXn2TiB4dCXzoPkxJl9wYoBBClA+pSVwnLy8vhvTqgq+PN9v3HsDLYmH0nE8wWa0EfTMbU3Y2GR9+4u4whRCiTKQmUQYmk4l+XTvQtV1L/PJyuRRhtPvl+vjwRee+xJ9JcHOEQghRNlKTKCOTycSAbp3o1q4127p24te53xKcnkqcTwAnfl7DqP49adU0xpjrOD8fKuncT0IIcTVSkygngf5+DOzeiXp33cWePoMAY/jZ0nVbWLFxO9bv5xE+uA/eWza7OVIhhHCeJIly1qFlM+6YNJrIsBqAkSj0nr0E/vUxvPVBwseNwG/+XDdHKYQQzpEk4QLBgYHcPGowDWrXBKDe6Xh8s7MAyImMInfocHeGJ4QQTpM+CRcJDPDnplGDiTt1lqX+frz5l+cZ/+McdnXqTu7OffTr2oEaIUWXTBRCiMpE1pOoAClp6fywch0XklKMDmyTCbOXF53bKHp3bkfwZx9jysgg694Hr7tju7o94HM9pIxKJ2VUOk8pI1lPohKpERLMrWOG0bJJI7BP1maxWtm+9wCz3/4A/+efIfil5wkb0g+vUyfdG6wQQhQgSaKC+Pr4MHpgb6aMHkqd6EjH/l7rV+KTY6xNkZyWzqWAIHeFKIQQRUiSqGD1akVz65hh9O/WEYAlN0xi0ejJZPv5M3f8rXy5ZBWbd8WSkpbu5kiFEKIS9kkopVoDDwEW4D2tdWxJ51eFPoninDiTQNzpsySlpHH20BEy/QMcx7xMJjq3VQzdtxPLoCGlThroKe2kriRlVDopo9J5Shk52ydRGUc33QucAeoCce4NxbUa1a1No7q1AUjs2JbFazaRlJoGgNVmI3HxUmq8/yo2P38yH3+KrPv/DF5S+RNCVBy3Jwml1MPAUPvmZqAR8BzQGZgBvOem0CpUragIZkwYxeETJ9l94DBnEs8z4YevMNlsmLKzSP5pIXroaJo3aUiAn5+7wxVCVBNuTxJa65nAzMvbSqmPgQzgEtWsz8Tb20yrpjG0bNKI3QcO8+Opadw4fzaRFxP5aviNJG3ewZrtu+nQshld2rYkODCg9IsKIUQZuDRJKKV6AK9qrQcqpbwwagUdgBzgLq31kau87APgIyAPeMKV8VVWJpOJTq1bEPbH2/mqdRv89+8jKSIKgLz8fHbEHmTnPk3PtAvUnTSJRvUqdmEaIUT14bKOa6XUE8BtQIbWuqdSaiIwTmt9u1KqJ/CU1np8We9jtVptFkvl6nwvTzabjbPnL3I8/iw79x3i3EVjKcu2e39l+lezONCyHcf/+je6jxspT3CXwmz2wmKxujuMSk3KqHSeUkY+Pma3d1wfBSYCX9q3+wLLALTWW5RSpQbnDIvF5hEjDUoS6BtIm2ZNad20CUdPnmbPlh2MW/gtAK0O7iVp1ge8nphCjeAgmjSsR6fWLQgPDXFz1JWPp4xKcSUpo9J5ShlFRzv3GeGyNn+t9XyMJqPLQoGUAtsWpZTb+0SqEpPJRLOG9Zk8ajDWseOxmUykhoSyfLhRIUtJz2DX/kN88t0iZi/8mSPxp9wcsRCiqqvID+lUoGDq8tJa51fg/T2GLSSUvJnvkDT1Ns4fj6NOeC3izyY6qsBeFgvJp07z44VLtGraiI6tWlAnOhKTfUoQIYRwVkUmiY3AWGCuvU9ibwXe2yNZuvcgonsPJgLBIX7EHohj14FD1P/6C/qsX8HCcbewz2bjwNETRISF0qtjW2P+KCGEcFJFDjFdAGQrpTYBbwGPVOC9PZ632UxM/TpMUjGMXL2YGqnJ3Db7Q9rt3QnApeRUFq/ZxOI1mzh3McnN0QohqgqX1iS01nFAT/u/rcA9rryfAPPJeAgOhqxMsporGH8jPifPkJdvtOwdPHaCg8dO0LxRfQb36kJwYKCbIxZCVGbScexh8vr049LG7QS9+DzZU6YyvFsPBuTmsWrzdg4cPeFYz+LwiVPEnU6gY6vmNKpbm/q1ozGbze4OXwhRyZT6nIRSKhzoA0QC54ANWuu0CojNKVV5gr/yVNqwPJvNRvyZRAKefpyk1DSWDx9Prp+/43igvx+dWiu6t2+Fl4fOD+UpQxddScqodJ5SRmWe4E8pFQ28ArQENHAW6AE8q5TaC/xda51YPuEKVzOZTDRNOEnYsoWYbDba79/Nf+/9K2mhYQBkZuewcedvnDl3ngHdOhERFiqjoYQQJTY3PYcxpcahKw8opVoBzwL3uyowUf785s/FZK85+nXszMDxo4k/m8jR+NOkZ2YBcPzUWY6fOkvDOrUY0rsrETVC3RmyEMLNKt16EtdKmpsMTlWBbTb8vp9H0L9eJHn+T441KixWKxt27GFH7MFCp5vNXvTs0IZu7Vp5RH+FpzQTuJKUUek8pYycbW5ypk+iD/AuUBs4BfxJa72rPIIsD5IkDNf0xs3PB+8rKpE2G3lP/IVtqj1bzQEUfF9EhtVgSK8u1K9ds0o3QXnK/9yuJGVUOk8pI2eThDM9lP8BpmqtawO3YyQMUZVdmSAA30U/Uvfz/zH+mUd49MA2akWGO45dTE5h7tLVfP3TCnYfOOwYTiuE8HzOJIlkrfV+APtSolU/hYrCLBaCX3gWAJPVSpDNxtSxwxnUozM+BRJKwoWLrNq8gzmLV5KTm+uuaIUQFciZ5qZvMBYBWg10AToBcwC01rNcHWBppLnJUNYqsPnYEYIffwSzPkjSxu3YahijnlLTM9i8O5YDR+KwWH+fHrlerWiG9u5KVHhYmWOvKJ7STOBKUkal85QyKs/mpoPASaAZxiR9a4E6GH0UwkNYmjQj5buFJC9d5UgQAKHBQYzo2YVHfHLp3kY59p9OPM/nC5by48r1nL+U7I6QhRAVoNQnrrXW/1BKjQVaALFa659dH5ZwC5MJa4OGRXYHzHqf4Of/xthOnYm++2EWJ2c5jh2JP8XxU2cY1LMz7VWzKt2xLYQoqtSahFLqXWAakA/cpZR6zeVRiUrD62Q8Qa/9EwCfXTvpdOoYU8cOp1nD+o5zLFYrKzftYO7S1SSlVJqH8YUQ5cCZ5qYOWuspWuu3tdY3YawwJ6oJa81aZD7wMDZfX/JbtSHr3gepEx3J+KH9mDFhJNEF+iROJZxj9sJlHDx2wo0RCyHKkzNJIl4pVR9AKVULo39CVBd+fmQ+9iRJazaT+u4s8PFxHIqOCGfqwB70qf37gka5efksXrOJH1euJzMr211RCyHKSUlzN50FbIA/cKNS6gTQALhQQbGJSsTSrPlV94e99jKj586hy8OP83nNGFIyjFEfR+JPkXDxEmMH9aFuzaiKDFUIUY5kWg4P4Y5hed47dxA2aohjPqgLX8zh5+Ao9h466jjHy2RiQPdOdGrdwu2d2p4ydNGVpIxK5yllVJ5DYIW4OosFS+MmAOQMGYZtxCiG9+3OxGED8Pf1BcBqs/HL1p0sXrOJ3Lw8d0YrhLgOkiTEdcvv1oOkNZvJePwp0l95A+w1hcYN6jJ9/AjqhoY4ztXH4/n6pxUy+kmIKsaZIbCPVUQgoory9yfz8accM8peViM4iP/7/nPuWTSH4LQUwJgD6osflrJ1z36sBZ7eFkJUXs7UJG5QSlX9eaJFhfJdtBD/lT8Ts+EX/vqflwjOMR7Ay7dY2PDrHuYtW016ZtVv1xXC0zmTJKKAM0qpLUqpzUqpTa4OSlR9Pr9ud/w7f9JNTLhpAtERBZ+pOM83i1ZK85MQlZwzE/w1unKf1rrSPC0lo5sMlXHEhc/aXwh841VSv5yDrUYYVquVrb/tZ/OuWMd6FV4mE6pJQwZ270xggH8pVyybylhGlY2UUek8pYzKvMZ1AfnAq0A08B3wG1BpkoSovPIGDCJlwCDHtpeXF706tqVWjVCC75zB5h79ONK8NQeOniApJZ2bbxhcaGpyIYT7OdPcNAv4BPAF1gFvuzQi4fHaLP2RNrE7uet/bzNx/peAsVbFDyvWkZwqzU9CVCbOJAl/rfVqwKa11oDMtSCuX24uAbPec2yGdujg+Hf82UQ+/2Ep8WcS3BGZEOIqnEkSOUqpEYBZKdUTSRKiLHx9SVq+luxJN5Pfui0Rz/+DHh3aOA7n51tYsGId+48cp6rPBiCEJ3CmAfj/gNcxRjk9Btzr0oiEx7PVrEna+x9DejomX1/6dmlP04b1WLhqPXmXLtFr8y+syM5i76GjjOjbg7ACD+UJISqWU3M3KaU6AQrYp7Xe6/KoroGMbjJ4woiLpJQ0Uu68nW7rVnAhsibfT5xOfIvWtG/ZjN6d2uJnn+rjenlCGbmalFHpPKWMym3uJqXUS8A7QA/gQ6XU42UPT4iiok6doOv6lca/L54jOD2VfIuFnfs03y9fK81PQriBM30SI4F+WutHgH7ATa4NSVRXllatSf/3TKyhNUjrO4CEQcMcx86cu0Ds4WNujE6I6smZPolTQAiQAvgAiS6NSFRfXl5kz7iDnBE3YLLkM71OXVZt3sGeg0cA2L18FXUvnieyV083BypE9eFMTaIucEgptRw4ArRRSm2S6TmEq9hq1cJatx4mk4kB3ToREhQINhuj5n5Os4mjSHn8UWxZWe4OU4hqwZmahDQvCbfx8fHmhgG9OfTGm7SAJnCMAAAgAElEQVQ6aIyZaPb5R/zQtCXtb5tmJBAhhMuUmiQq0zxNonqqXzuayFsmcmb9cuoeO8yWHv3Z4hvMrgVLmDxyELWjIt0dohAeSxYdElVCQNdumNZtYdt9f2HZyBsByMnNY8HytSSnpmG6cAFk9JMQ5e6akoRSSpKKcBtvXx8aP/884yeNdSyPmpmdw9yflhM8cQw1Jo7BfOSwm6MUwrM485zETUqpKUqpPwAJslKdcLf6tWsyfmh/zGbj7dtxxWL8D+7Hd+N6wkYNgfR0N0cohOdwpmbwGLACmA40AMa6NCIhnFC/djQThw3Ax9sbv5xsLF7GW3nfxCnYgoLcHJ0QnsOZ0U2XJ/RL01rnKKVkIh1RKTSsW5tbbhjCokB/9rXpSN8Nq5jfrD1NVm2gb9cORIaFGifabGAyuTdYIaooZ2oSx4EdwCdKqeeAra4NSQjn1YqK4LbxIwno0YN5N9+O1WzmSPwpZv+4jJNnEzGlphA+qA9+334tHdtCXIdSk4TW+nagndZ6EfAh8LSrgxLiWvj6+DBucF/aNG/s2JdvMaYctzz9V7z3xxL64D2Y/zDDjVEKUTU503H9X6315Z7ADkhNQlRCZrOZkf16Mm3scIICAoydmRn4rPjZcY51jHSnCXGtnOmTSFVKvQIEA20xJvwTolKqHR3JTSMHMW/ZajKAmQ8/y7AVP1Ez6QJBg4bg5+4AhahinF1P4t8YTU6VLkHIehIGT5njvrxkZmWzfOM2jsafBsBktRJaI4RbbhjqmMrD62Q8gf99i4yn/o4tPMKd4VYa8j4qnaeUkbPrSRSbJJRSZ4HLB01ALSABQGtdt3zCLDtJEgZPeeOWJ5vNxpH40yxbt5ncvHwAggMDGDe4L3WiIwmdfjN+K37GGhVF2tvvkTus0n0HqnDyPiqdp5SRs0mi2OYmrXWdco3ISUqph4GOQHPgK631e+6IQ1R9JpOJ5o3q4zO4LwuWr8Vqs5GemcW3S1YxMchMF3t/heniRaxh4W6OVojKqdQ+CaVUR4x1rv0v79Na/9FVAWmtZyqlfDFWw/vAVfcR1UdMvTrcOHwAS9ZsIisnF4vVyrw0KxlP/oM+X3xI7vBR5Hfr4e4whaiUnOm4/gzjA/ukKwKw1xyG2jc3a63/CdwKfK+1trrinqL6ialXh3umjWf2ghWcT0oGYElYbWKffpURfbsTesX53nv3QGYW+T1kgSNRvTmTJBK01h+7KgCt9Uxg5hW7+wF/ctU9RfUUERbKrWOG8fOGrejj8QDEp2fy5epNTBoxkNrR9inH8/MJ+fN9eO/bS9Ztd5Dx9+exSXOUqKacSRJxSqkngV3YO7K11sudubhSqgfwqtZ6oH0G2fcwnrXIAe7SWh8pLi6ttTweK8qdj483owf2pnZ0JOu378Zqs5Gdm8u8Zb8wecRA6tSMIuCjD/DeZyxw5P/dHDL//IgkCVFtOZMk/ABl/wEjUZSaJJRSTwC3ARn2XRMAf611L6VUT+ANYPzVXmt/ytspZrOJsDBZncxs9pJyKEXBMhratzOqSX1m/7CczOwccvPymL98DdMnDCNs2hSs2zbhtXgR1r8/S2iH1m6OvOLI+6h01a2Mik0SSilvrXU+cPd1XvsoMBH40r7dF1gGoLXeopQqdeiVMywWm0cMRysrTxmW50pXllFIYDCTRw5i3rJfyMrOISc3j//NXULrZjEMev9TQjeuJbf/ILiiXM2HD2GJaQw+PhX9K7icvI9K5yllFB3t3FytJU3L8YX9vxo4aP+5/O9Saa3nA3kFdoUCKQW2LUopZ2oyQrhMdEQ4N48aTID/789i7z8Sx1c/LedE+y5FEoEpNYUaE8cQPrQf3ttlhhrh+Up6TmKq/b+NizvnGqUCBVOXl72mIoRbRYWHMW3scNZs3cWR+FMAJKel8/VPy+nStiX9u3bAy75eRdDLL2BOTIDEBELvnMGl7b+Bn0z2ITxXSc1Nm/n9ietCtNa9r+NeGzEWLJpr75PYex3XEMIlaoQEM35oP/TxeH5ev5W8fOP7y6+xB0lJTeeGAb3w8fHGEtMYW2AgpsxMMl78lyQI4fFKau6ZUs73WgAMU0ptwpjm445yvr4QZaYaN6R2VAQrN20n7nQCAEfiT/HlwmWMGdiHmvc8QM7ocfjP/YaccTcWvUBWFlyehVYID+DUBH+VmczdZPCUzjRXupYystlsrN2+m19jf++CM3t50bdrBzq1boHZq2h3ntfJeMJHDCTzvofIuvu+KtmxLe+j0nlKGTk7d5MzK9MJUe2YTCYGdu/EiH498PE2KtwWq5W123bx+fdLSLxwqfALbDaCn3wUrwsXCH7h74Q8eL2DAoWoXIpNEkqpf9r/e9VnGYSoDto2b8L08SOoGfn7w3RJqWnMW7aauNNnHftMKcmYT592bGfdKUlCeIaS+iTGK6XOAA8qpWoVPKC1nuXasISoPCJqGNN57Nyn2bpnP7l5eeTk5jH/5zXUqxXNmEG9CQ4LJ2nFWgI+fA+vc4kyYaDwGCU1N90BhGE8cV3nih8hqhVvs5nu7Vtzyw1D8Pfzdew/nXieect+ITMrG3x8yHrgITJeeLno63fuIHTGFLxOxldk2EKUWakd10qpbhhPTzcFjmutL1REYM6SjmuDp3SmuVJ5lVFyahqbd8ey/0icY19IUCDD+3SnUb3amEymwi/Izyd82AC89+3FFhhI6rsfkTu6cq63Le+j0nlKGZVnx3VjYAvwN2CLUmp6GWMTokoLCw1hVP9ejB7Y25EQ0jIymb98DZ8vWML5S0mFzvfZsA7z/lhjw2Yjv227ig5ZiOvmTJJ4BOistZ4AdAIecm1IQlQNLZs0YvyQfvj7/t78dDE5lblLV3Mq4ZxjX97AwSQvWk5+q9ZkPPok1kYxbohWiOvjTJKwaq3TAbTWaUC2a0MSoupo2rAeM24cRYeWzRzPTmTn5PLtklUsXL2BSympAOR360HSyvVk3ftAkWv4rliG76KFUMWfWRKeyZkJ9o4qpd4A1gH9MfonhBB2IUGBDO3djXYtmvLdz7+QnZMLwOG4kxw9cYq+XTvQtW1LTFd5uM6UmkLwX/6MOTGBnJE3kPbvt7HVqlXkPCHcxZmaxB+BY8Aw+39lxTghrqJWVATTx41ANW7o2Ge12Vi3fTfzlq7myIlTXDlQJHDmG8aEgYDPrzvAX+aCEpWLTMvhITxlxIUrVWQZJZy/yOotv3L2/MVC+5s0qMuQXl0JDQ4CwJScRNCLzxPw5aekzvqUnAmTKiS+4sj7qHSeUkbOjm6SJOEhPOWN60oVXUYWi4Vftu5iz8HDhfZ7m810at2C7u1bO565MO/fh6VVa7hi+KzPhnXkd+yELdi5BWLKSt5HpfOUMiq3IbBKqcCC/xVCOMdsNjO0d1fuumksHVo2d+zPt1jYvvcA//vuJ3bEHiTfYsHSuk2RBOF1Mp4a028mvG93fJcurujwhQBKnrvpZaXUBODyu3NRxYQkhGepERLM0N5dueWGIdSKinDsz87JZe22XXzy3SJiDx8r3F9hnzDQlJmJ+cxpgl55CfJljS5R8UqqScwHagPtlFIrgJZKqQfsCwYJIa5R/do1mTZ2OKMH9qaGvU8CjAfxfl6/lSVrN5Ob9/uKvzk3TsYaFYXNZCLt9ZngLav9iopXUpIIAmYDv2mthwHHgXhgREUEJoQnMplMtGzSiDsmjWZQj84EFhjNdPDYCT6dv5jf9FHyLBZyJt/CpY07SHt31lUnDPRKOFtknxDlraQk0Rv4HGivlHofY7K/A1rrf1RIZEJ4MLPZTOc2irtuGke7Fk0d+9Mzs1ixcRsfzV3IgaNxWMPCyZl8S5HXe+/cQUSXtgQ99zdIT6/I0EU148wEf79gPCuxGFgFNNdaj6yA2Jwio5sMnjLiwpUqcxkdOBrHmm27jNlkC6hbM4qeHdoQU7/O7xMH5uURPnwg3vuMZeKzp0wj7T/vl0sclbmMKgtPKaMyj24qMJrpUa31cWC61vrBywlCRjsJUX5aNY3hj5NG079bR4IDf18j+8y5C3y/Yi1fLVzOCfua26bUVKwRRge4LSCAjEf/6paYRfVQUk/Yu0qpHcAcAK31TgClVBQwHWOyvz+4PEIhqgk/X1+6tWtFh5bN2PDrb+w5eASr1QpA4sVLfPfzLzSuX4c+ndtjnfcj/t99iykjo+iEgTab8XOVdbiFuFYlNjcppW4GHgQaABeAUOAs8J7W+tsKibAU0txk8JQqsCtVtTJKTc9gR+xB9uqj5FsshY6FBgfRtkUTOqhmBAb4Fzrm+9OPBL73Nmn/fhvLNU5LXtXKyB08pYzK9YlrpZQ/EA5c1Frnlj288iNJwuApb1xXqqpllJaRyaade4k9fKzIMX8/X0b260nThvUAY8LA8D7dMCcmYDObSf34i2ta4KiqllFF8pQycjZJODXwWmudjVGDEEJUsJCgQEb060HnNoptv+3n+Kkz5OQaz1Nk5+Tyw8p1NKpXm86tFc0P7sUr6RIA1sgo8vr2c2fowgPI0zlCVBHREWGMHtgbi8XCobiTrNu+m/TMLABOnE7gxOkEAvz96DbzIwZ98SF5d/4fthphbo5aVHXOzN3UpSICEUI4x2w206ppDDMmjKJN88aF1tTOys5hXWo2L0/6I8vrxDgWPbrM/9OP8f/0Y7B3iAtRGmeek5gDxGA8fT1ba51cAXE5TfokDJ7STupKnlpGl1JS2XvoKAePnnDULAqqVyuafl070CA3i4h+3TFlZpLXpRupn32FtVbtQud6ahmVJ08po/LuuA4HpgITgHPAR1rrNWWMsVxIkjB4yhvXlTy9jKxWKydOJ7Bux24uJKUUOX7zorl03rAKgNyWrUlZtR6uWC3P08uoPHhKGZXbVOF2tYCGQBTGUNiblFKfXH94Qojy5uXlReMGdblt/EhuHNafZo3q41WgKer7ETeycsho8rx9+HTIeJZv3cX5S5WqYUBUQs40N20FMoGPgPla6xz7/p+11m6f7E9qEgZP+XbjStWxjFLTM1i5aTvHT/0+ODEoPZWM4FDHdv3aNWmvmtJu6waCJ44j2T/0apcSdp7yPiq35ialVDet9fYC2wO01mvLHmL5kCRh8JQ3ritV5zJKy8gkMyub0+cuEHvoaJEaRP2Tx7n/vVexBIeQ+ewL5My4o8giSMLgKe+jMj8noZTqB7QGHlFKvWnfbQbuB9qWR5BCiIoREhRISFAgtaIi6NSqOacTL7Brv+bwiVNgsTDx+68w2Wx4p6Vy7pP/saVxa5rG1CemXh28zWZ3hy/cqKTnJJIwFh3yA+rY91mBJ1wdlBDCdUwmE/VrR1O/djRpGZnsP3KcXyZOZcScTwlNTWL+2CkkHTlO7JHj+Pn60L19a7q2bYmXzAVVLTnT3FRHa11pn7aW5iaDp1SBXUnKqHi5eXls3fIrqevWcLBOo8IHbTbapF0iauhQOrRshk81XyHPU95H5dHc9J3WejKwUyl1OZOYAJvWum75hCmEqAx8fXzo168noWMGcfDwSY7Gn+LgsXhS0zNoG7uL6V99yK4lP/DdrbfTc9hgYurVKfQQn/BcxSYJe4JAa12nuHOEEJ7Fy2Sibs0o6taMonfn9vy6aSv9/mm0MHfavY2MoBC+t5mpGRlO93atadqonvRZeLhS641Kqf5AIMYzFf8F/q61/trVgQkh3Mvs5UWP1i0wDxkKP8wnLaQGK4eOAeDcxSQWrdmIj7c3rZvF0LdLB/z9fN0csXAFZxoXXwOmAe8CfYC5gCQJIaoBW0QkmbM+JW/KNLIyMmgZVY/YQ8cc61tYcnLYt/cAh+JO0rFVczq0bEZQQEApVxVViTNJIgtIBPK11glKKT8XxySEqGTyBg/FGxgC9OrYll37D3Hw2AnaLf6eHlvXs+DGqWzOzmHbnv20ahpDr05tCQ0OcnfYohw4M6YtFVgJzFVK3Q/EuzYkIURlFhjgT58u7bmreztGrl5C1MVz/OnjmbTavweL1Urs4WN8Mn8Ra7ftIjMr293hijJypiZxM9BUa71fKdUWY3oOIUQ153PkECZfX8jOIqNpc1L7DgD7k9wWi5UdsQfZuf8QLWIa0K9rB6lZVFHOJImawBil1OQC+15wUTxCiCoid8hwkjZuJ+jvT5L7p3uZ0q0HJ88msm77HhIuXASMmWkPHjvB0fjTdGzVnE6tWxASFOjmyMW1cCZJzMNobjrp4liEEFWMtVZt0mZ95thuUKcWU8cO43DcSfz//hRJFitrBo4iD9i+9wC7DxxmcM8uRRZLEpWXM0kiTWv9jMsjEUJ4BJPJROtLiYStWITJZqNz7E7e/dNfyAgOJS8/n583bOXoydMM692NwAB/d4crSuFMkohVSk0BdgE2AK31IZdGJYSo0vy//AyTfcqfoNZtGDr+Btbt2ENSShoAR06cIv5MAt3bt6FzmxbVfqqPysyZv0xH+89lNmCwa8IxpiYH7gBCgTe01rtcdS8hhGukv/Ef8tt3JPDN10h/7U2aNWpAo3p1WLttN3sOHgYgNy+fDb/uYc/Bwwzv252YejK5Q2Xk7PKlNYBGwDGtdborA1JK3QNMsW9O1lpfKOl8meDP4CmTjrmSlFHpyr2McnLAr/CjVXGnzpD7zNNsbtWBC9G/r7HdpY2ie/vWlb4JylPeR2We4O8ypdQk4Bn7uXOVUjat9UtlD9Fx/YeBofbNzcAK4BOM2spo4PPyupcQooL5FX32tsWubdRY8j29VvzE2qFj+Ln/CDCZ+HWfZu+hYwzr042WTRpd5WLCHZxpbvoL0BNYBrwE7LD/t1xorWcCMy9vK6XmANOB8xiLHgkhPEVeHsHPPg2AV14eXSJDOVC3NvFnEwFjyvLFazZx6Hg8fbt2IKKGLKXqbs48cW21r2tt01rbgAxnL66U6qGUWmP/t5dS6gOl1Gal1BqlVLNiXjYH+BL4MzDb2XsJIaoAHx9SvphDXqfOWGrWIve5F5k8chDjBvelRkiw47TDJ07xxYKlbPttP840iQvXcWbRoZeBxkAXYDWQobV+tLQLK6WeAG6zn99TKTURGKe1vl0p1RN4Sms9vqy/gNVqtVks8iYym72wWKzuDqNSkzIqXYWVkcUCx45B8+aOXdk5uSxesYG8779nX5uOjjW227ZozKSRA/D2rhxTknvK+8jHx1w+fRJa66eVUiOBncBBrfVPTsZwFJiIUSsA6IvRZIXWeotSqtTgnGGx2DyiE6msPKUzzZWkjEpXoWUUXQ+uuNeY37YTPPsD4lu149vRN3Exqhaxh46TcP4S/bp2pEkD96935invo+joEKfOK7G5SSnVXin1AjAO8AEOOhuA1no+kFdgVyiQUmDbopSSwdFCCAC8TsYT9No/AWh4YC/Dzv0+ycOFpBQWrFjLpl17pfmpghWbJJRSN2GMMjqJUQNIA75TSl1vE1EqUDB1eWmt86/zWkIID2MNjyBrxh+xeXmR36oNdV/+FwO6dSz0oN3mXbF8u2QV5y4muTHS6qWkb/IPAQO01o6OaqXUZ8BC4MfruNdGYCzGMNqewN7ruIYQwlMFB5PxwsvkTL4ZAJOvL13btaJ1s8YsWbuZs8fjMFssnE6ErxctZ3if7rRu1tjNQXu+kpqb8gsmCACtdRpguc57LQCylVKbgLeAR67zOkIID5bfviP57X+f5CEwwJ8bhw9gxvY1PPrmc3TYvQ1LvoWl67awZXesND+5WEk1ieK6750ZNguA1joO4xkLtNZW4B6nIxNCCDu/3TuptfgHTDYbt875HxlBIRxp3oqNO/eScOESfTq3JzoizN1heqSSkkQbpdSVa1mbkAfchBAVzJSejrV2Hcxnz5A9cAi5/QfA2XMAHI0/zdH407Rs0oghvbri7+fr5mg9S0lJ4uZi9n/gikCEEKI4ef0HkrRxO4Gv/YusO/+PG+s3YNWmHcQePgaAl8XCwWMnSMvIZPKIQZXmmQpP4NQEf5WZTPBn8JSx264kZVS6qlZG5y4msWnnb/R56RnSQkJZNvJGopo0ZkjvbkSGuWZKj6pWRsUptwn+hBCisqoZGc7N2cnUOLAHgFYH9vDGoy/w5fkLjB3Ul6YN67k5wqrP6U5oIYSojHzXrHb8e3+bTuT4B2CxWFm4aj37Dh93X2AeQpKEEKJKS399JimfzCavQydqvP1fx0SBVpuNZeu3sGHHHhkmWwaSJIQQVZvJRO6YcSQvX0NU4ximjB5CZFgNwOjQjnn6UQ7MvnKgpnCWJAkhhGewzxobHBjIrWOG0aRBXfpsXE2n3dvo99h9JN19l9QoroMkCSGEx/Hz9WFcry4M3rACAC+bjaNZeSxbt4X8/OudNKJ6kiQhhPBI5uBg0lauI75TN87Wrsf6fsPYfzSOuUtXkZGZ5e7wqgxJEkIIj2WKaYzfkuVsePW/WM3GA3Znz19k9k/LOR8XR8BH70NeXskXqeYkSQghPJqX2cyAkUMY2KMTJnu/RXpGJimPPEzw3/5K+PCBeO/61c1RVl6SJIQQHs9kMtGlTUsmDhuAn68PdU/H02PTLwB479uL+UScewOsxCRJCCGqjZj6dbh1zHCym7dg2cgbyfXx4XDLdsT16u/u0CotmZZDCFGtRIaFMnH0ML6y2tjbrgtWLy/Sl65mVP9eqCYNAfA6cxosFqwNGro5WveTmoQQotqpERLMhKH9yapdl5SwCCxWK0vWbebEmQSw2Qh+/GEi+vUg4P13IN89qyxbLBZS0zPcPmRXZoH1EJ4yM6UrSRmVrrqVUUpaOt8vX8ullFQAfLy9ucOUTZNH73eck7R8DfkdOzu2XVlG5y4msWbrThIvJpFrH3UVEhTIjAmjyn2dDGdngZWahBCi2qoREsykEQMJDgwAIC8/nxVxZ8hs1gKArD/cWShBuMrZcxdYum4zX/20nJMJ5xwJAiAtI5NTiedcHkNxpE9CCFGthQYHMXnkYOYuXUVmVjbH6zXkn398mDvOHCXsoYeKviA5GSifb/VWq5VNu2LZumdfiedl5+SWy/2uh9QkhBDVXmRYKDePGkJIUCAAFi8zH9dvwb7zSYVPzM/He+gQQmfcitfpU2W+79rtu4skiNpRkUwfP4IOLZs59uVIkhBCCPeKDAtl6tjhREeEOfYtXbeF7XsPOCYGDJj1Pqbf9uC3bDHhwwZA1vVP73H4xCl27tOO7YZ1ajF1zDCmjh1GrcgIAvz9Hceyc933VLgkCSGEsAsODOCmkYOJCq/h2Ldu+25+3rAVi8WC19nTjv2Z9zwAAQHXdZ+z5y+ydO1mx3bThvWYPHIQdWpGOZ4K9/fzcRzPyZWahBBCVAoB/n7cNHIwdWtGOfbtO3yc71esJeX5f5L/yxqyx91I1r0PFH2xE6NFLyQl8/3yNeTZh9aGBAUyol8PR3K4zN/3934P6ZMQQohKJDDAn5tGDaZ1sxjHvvgziSxduxlLr96kffw5+PgUeo0pNYWwG4bgu2xJsde9mJzKvKWrHR/6/n6+TBo+kAA/vyLn+hVIElKTEEKISsbbbGZkv5706tjWse9Q3ElWbdp51fOD/vUiPr/uoMaMKQQ//kiR48mpaXy3bDWZ2TkA+Pp4M2n4QCILNG0VVPC5CKlJCCFEJWQymejVqS2dWrdw7Fu3bQ8Hj50ofF5aKr4//ejYzuvTt9Dx1PQM5i5dTbp9HQsfb28mDh9I7ejIYu8tNQkhhKgCTCYTA7t3IqZeHce+peu2EHf6rGPbFhJK0vqtZE2bQc6QYeSMn+g4lpaRydylq0nLMJ7S9jabmTCsP/VqRZd4Xz/fgh3XMrpJCCEqLS8vL0YP7E1EjVDAeAjuh5Xr2H/kuOMcW3gE6W+9Q+oXcxzrbWdkZTFv2WpM8ScYuXQ+AXm5jB/Sj4Z1apV6T2luEkKIKsTfz5dJIwZSIyQIAIvFytJ1W/hx1XpHLQFwdGhfTErhm0UrSUpOZfyP3zBw7XKefPdlmscddup+Pt7eeNmTTb7FQr7FPRP9ybQcQgjhpNDgIO68eTSfz/+Zi8kpABw5cYrEC5cYO6gPYHzrP5+UzJbd+8jLz6fpUU2rg3sB8Es4S+ZVRjJdjclkws/Plyx7R3dOTi7egdf3XEZZSJIQQohrEF4jhFvHDOOXrb+y77DR3JSWkcnXi1Zc9fx41YbDz/+Lpv/5NzljbyS/Ww+n7+Xv+3uSyM7NI0iShBBCVH5+vj6M7NeTFjENWbhqPRar9arnRYbV4IYBvQiLDOfSLVPAp+hHrjnWqGVY2ra76n0uc9cIJ0kSQghxnZo0qMuEof35ZetOrFYr/n6++Pv54ufrS92aUbRv2QxvsxkAW+RVhrvm5xP64D2YD+4n654HyHjsSQgKchz2qwSd15IkhBCiDGLq1+GO+qOv67UBs97He59Rkwj4ZBZZt9+JtUCS8K8Ez0rI6CYhhHCT3GEjyO1tPHiX8dhTWBvFFDruVwnmb5KahBBCuImleQtSFizG96cfyB01psjxiEvnMVmt2Ly8pE9CCCGqJZOJ3HE3Ft2dmsLAJx6gVVAo30+cTnZuSzcEJ81NQghRKQW9/AL+Fy/QKP4Yd3z6X5LOX8TihgfqJEkIIURlY7NhjYjE6m009iwePZljieeZvfBnTieer9BQJEkIIURlYzKR+cTTnF+5jm2jJ7K3XRcALiSlMGfxSlZs3E52WnqFhCJJQgghKimv1m1p9L9PGNCjM97eZsf++C1bqdG5DeeWFL/AUbnF4PI7CCGEuG5eXl50bduSOyaOpkmDemCzMf7Hb8j38mL+6UsuHxorSUIIIaqA0OAgJgztx8TObah1IZEfJkzFGhiA71Wm+ihPMgRWCCGqCJPJRONOHcjavoc2l1IYGBGGl5drv+tXuiShlJoCjAEuAs9ordPcHJIQQlQqvkFBNCkwfYcrVcbmpnHA7cAXwAz3hiKEENWb22sSSqmHgT4H7lEAAAe6SURBVKH2zc3Af4GPgJOAe5ZiEkIIAVSCJKG1ngnMvLytlJoI3AUMAhq4Ky4hhBAuThJKqR7Aq1rrgUopL+A9oAOQA9yltT5ylZclA59h1CLudmV8QgghSuayJKGUegK4Dciw75oA+GuteymlegJvAOOvfJ3WejWw2tn7mM0mwsICyyHiqu3/27v3GDvKMo7j32271QKtG2SFJqQ2UfjZGG8obZU2bFRaSAONaEIFL9VEIUWlQaNoYuQiJKDxD0FR1AhqoS3gpWpaTKgLSmlARBQbHpPW4o22rrXUS6GyXf9435Mdpjuc3e6eS8/+Pv+cuc97nt2Z57xnzjwzdeoUx6EOx6g+x6i+yRajRvYktgPnA9/N44uATQARsVXSmyZiJ4ODQ+zb99+J2NRRrafnGMehDseoPseovk6JUW/vzFEt17BfN0XE3cD/CpNmAU8XxgcltfyaiJmZVWvmSXo/UExdUyLiufFutLt76kBv78wnx7udTjDaTwaTmWNUn2NUX4fE6OWjWaiZSeIB4Fxgfb4m8bsJ2m7vBG3HzMxKmpkkfgCcJWkL0AV8oIn7NjOzI9A1NDTU6jaYmVmbaseyHGZm1iacJMzMrJKThJmZVXKSMDOzSr6Zrc2V6l+9klTXagh4HLgUmAOsjYiFeflFwLeBdwKX53mbWtH2RqsXm4g4JOlzwDLgOWB1RDwk6VZyXPINnbcDA8AXgDtqsexUpbjdSiqZc2JEPJvnnwY8QiqyORd4VURc0aLmNs0oj7U/AldExPWF9TYAs/J6O0nxeqbJzW8Y9yTaWK5/9U3gxXnSl0gPYlpM+hnx8tLyfcAtwLKI+G0Tm9p0o4lNPtmdCSwAVgBfKW2jG1gP7IiIVaQTQkcbIW4ATwHnFMYvAnY0s12tNoZjbTvwrsJ6xwOnNLGpTeck0d5q9a9q3gjcl4c3MvwcDiS9nVRl9+yI+ENhnVWS7pV0X/501ClGE5tFwM8iYigi/gRMk1S7+fJFwPeB35Q+JfdK2iBpq6TPNvYttEQ5bgB3AO8GyNWaTwMeLsx/c/4feljSsuY0s+lGe6wNAHskzcvjFwB3lrb19Xy8rZc0o1ENbhYniTY2Qv2rroiofdr9F/CSPPwK4FrSp6DyP+WWiHgbcD1wQwOb21SjjE25XlgxZl8GjgVOLm36OFL14jOAcyS9boKb3lIjxA3gIUCSjgXeCvy8NP8/pJPkMuCmnEg6yhiONUhJdUUeXg78sLS5myPiTGAn8KGJb21zddwfu8MdKgzPJD17A+AA6euCj5HKnhQTxf35dQughrewdUaKTbleWDFmNwJLgNdIek9hmcci4umIGCSdPE9tXJPbygbSCe9CYE1p3i9zb2wPKem+tNmNa4GqYw1SUjhP0lxgF1AsCXswIrbm4Y445pwkji6P5usOkJLCL/Lw3yJib0T8JE+7qbDO/Py6mHQBrlONFJsHgKWSpkiaQyoqOZCXeTwXmLwIuKHw9cE8ScflC9oLgN837y201BrSM+VnR8T20rzTASSdROppDdD5qo41IuLfQJB65reX1psu6fV5uCOOOSeJo8vHgaskPQhMB+4aYZlPAPMlvT+PL5S0GVgNfLI5zWyJw2ITEY+QDu4HgbtJv1B5nojYQYrLncAxwF5gHelT4F0Rsa05zW+tiAhSscwfjzB7Rv4f2gBcXPgappPVO9bWkJLAvaXpzwIflXQ/8DLSxfCjmms3mZlZJfckzMyskpOEmZlVcpIwM7NKThJmZlbJScLMzCo5SdikIKlP0toRpq+VNP0F1ts1xv2cJOmrR9LGwjYWS7qsqs0V68yQdJukrvHs26zMVWBtUouIFfWXGtP2dgGrjnT9fJK/knQD11vGsN8D+fnx7wNuO9L9m5U5SdikVivtDHyNdCPUXGA2sDIifl1Y7jpS/Z6P1G4my8UC15F65N3AJaQ6P2tJZbY35tWnAQtJJRpOJtXZGiQVlbs4Ioo1g84CtkXEQUkAp0i6h1QK4+aI+JakfuCJ3O4u4IKcnNYDm3CSsAnkr5vMhj0ZEUtJdZ0+XJso6YvAtIi4tHS38XxSLaNa3axZtRkRcSAi+kjJYiepd7ED+AZwfi4A91dgZakNfUCxzHs3cC7p7t5PFarYbsnbXwd8Ju/zn8AJkorF6MzGxUnCbNij+fXPDD9X4ETgtaSaRWUbSeWkfwRczfOLwtXcSKp6cQup7MVsUhHGflKBwTml5U8AdhfGt0bEwYg4AGwj9XQANufXchG53cDxle/QbIycJMyGjVSjZjewFHi1pLNL8/qApyJiCfB54LriTEnXkEpOX5MnDQB/AZbnXsC1HF6Wew/QUxh/g6RpuYz3PNJXVJCedwCppHmxCGEP8PcXeI9mY+JrEjaZLJH0q8L4haNZKSKGJH0QuEfSgoj4R571GLBO0mrSNYara+tIOh34NNCfew3k+ZcBP83PZNhPutBc1A+8A/hOHn+G1GPpAa6MiL35WsVKSZeTnvXw3rzPHmBfrlJqNiFc4M+sjeTksRlYEhEHK5bpBy6JiCdK01cB+yPiew1vqE0a/rrJrI1ExCHgKsb4M9r8oKkzOPz5Bmbj4p6EmZlVck/CzMwqOUmYmVklJwkzM6vkJGFmZpWcJMzMrNL/Adk7VtNYhnZbAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "#---- Copied from jcvi.graphics.base\n",
    "from functools import partial\n",
    "import matplotlib.ticker as ticker\n",
    "\n",
    "# human readable size (Kb, Mb, Gb)\n",
    "def human_readable(x, pos, base=False):\n",
    "    x = str(int(x))\n",
    "    if x.endswith(\"000000000\"):\n",
    "        x = x[:-9] + \"G\"\n",
    "    elif x.endswith(\"000000\"):\n",
    "        x = x[:-6] + \"M\"\n",
    "    elif x.endswith(\"000\"):\n",
    "        x = x[:-3] + \"K\"\n",
    "    if base and x[-1] in \"MK\":\n",
    "        x += \"b\"\n",
    "    return x\n",
    "\n",
    "human_readable_base = partial(human_readable, base=True)\n",
    "human_formatter = ticker.FuncFormatter(human_readable)\n",
    "human_base_formatter = ticker.FuncFormatter(human_readable_base)\n",
    "#------\n",
    "\n",
    "def plot_distribution(A, B, title, label, filename):\n",
    "    df = load_hic_data(filename)\n",
    "    regression = lambda A, B, X: A * X ** B\n",
    "    ax = sns.lineplot(x=\"BinStart\", y=\"LinkDensity\", data=df, lw=3, color=\"lightslategray\")\n",
    "    tx = df['BinStart']\n",
    "    ty = [regression(A, B, x) for x in tx]\n",
    "    ax.plot(tx, ty, 'r:', lw=3, label=label)\n",
    "    ax.legend()\n",
    "    ax.set_title(title)\n",
    "    ax.set_xlabel(\"Link size (bp)\")\n",
    "    ax.set_ylabel(\"Density (# of links per bp)\")\n",
    "    ax.set_xscale(\"log\", nonposx='clip')\n",
    "    ax.set_yscale(\"log\", nonposy='clip')\n",
    "    ax.xaxis.set_major_formatter(human_base_formatter)\n",
    "\n",
    "A, B = 0.00162, -0.8135\n",
    "title = \"Arabidopsis Ler0 (RE=AAGCTT)\"\n",
    "label = r\"$Y=0.00162 \\times X^{-0.8135}$\"\n",
    "filename = \"ArabidopsisLer_RE_AAGCTT.distribution.txt\"\n",
    "plot_distribution(A, B, title, label, filename)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAESCAYAAAAIfCk9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3Xd4lFXawOHftPRAAkmoCZ1D70jvvShFwd57WVf97Gtby65dVHTVtS0qKoIiIFJFeui9HFroHdL7lO+PGWISMskEkpkk89zXlYuZtz5zmJlnznvec47B4XAghBBCFMXo6wCEEEJUXJIkhBBCuCVJQgghhFuSJIQQQrglSUIIIYRbkiSEEEK4ZfZ1AEJ4QinlALYDNsABhAApwP1a6/VKqfuACK3162V0vnrAf4AxwIvAg8Ax12oDUA34Bfg/rbVDKfUn0ABILnSoV7XW0z04X1/gGaAJYAcygXe01t8W2u5KYBZwndb6R9eyCOBP1yZhQD1Au54v1Fo/oZQKAv4BjHbFbwK+Bd50xf8IcF5rPaWkWIV/kSQhKpMBWuuzF54opR4HPgR6aK0/KeNz/Rd40fUFCvCj1vqhfOeOBLYC811/AE94khAKU0qNAD4Drtdar3AtawjMV0plaK1/zrf5A8B3wKPAjwBa6ySgg2u//sBkrXWHfMc3ADOBPTjLKkspVRP4DWdSeR5nOa5TSi3QWp8s7WsQVZckCVEpKaXMQBxw3vX8JSBKa/2QUqo58CkQg/NX+ata6x9dtYPJrv0swA9a638VcexuQIzWel0xIdTCWZtJ9CDW/L/08/tJa/0a8Cbw6IUEAaC1PqiUugsIzXecxkB/nDWWXUqp7lrr+JLOD/QFWgKjtNY21/HPKaVuBhq6ntuUUtOAp3AmICEASRKiclniuuwUBWQBc4Dbi9juB+BzrfXHSqlY4E+l1O/AN8B7WuvZrssvc5VS+7TW0wrtP8F17PyuVUr1xvmlXQPYCNyrtV6bb5u3lFLPFdpvkNb6HK5f+oW5Ekgb/qqN5NFaLy+06D7gN631aaXUDzi/zK8t6riFdAHWXEgQ+Y6/F9ibb9ECnJfQJEmIPJIkRGUyQGt9VinVCZgLLNFan86/gVKqBtAe+BxAa30EaKKUCgX6ATWUUq+4Ng/D+eVdOEm0wJlo8vvRVUsJwFkbaQXMLrRNkZebiqtJAB+7Hjvybf8joIAA4LTWur9SKhBnQrzDtdn/gJVKqVjXayyOHc9uUkkA4pRSQVrrLA+2F35AkoSodLTWG5VSjwJfK6U2aa0P5lttdf2b/0tXASdxNtj21FpnuJZfqJEU5sDNl6rWOkcp9RCwAXgLeNiDePPaDIqilNqF8zLSHNf217qW98eZkAAmApHAZKXUh/ni/BvwZAkhxAOPKKVM+WsTSqmuwMNa65tdi3Jcx7SX9JqE/5BbYEWlpLX+HlgLvFdoeQrOL/BbAVyXm1YCwTi/LB9zLY9wLR9T1OFx3mXk7tw5wP3AA0qpjpf7WlwxfaCU6nlhgVKqGs47kS58qd8PvKa1bqC1bqi1bojz8tPdrlqSW1rr1cBu4F3XZTaUUrVwNlYn5Nu0MZDgen1CAJIkROX2EDBSKTWs0PIbgIlKqS04Lwnd5bpj5wagu1JqG7AG+F5r/V0Rx50ODC/uxK5G5u+Aj1x3D4GzTWJzob8XSnoRWut5wPXA00qprUqp3cA6nDX90Uqp9jhrIh8W2nUKzobz20o6B3A1zprUBle5LAZm4Ly994LhOC+BCZHHIEOFC3ExpdR84PlCDdNVllLKhLMxfqjW+pSv4xEVh9QkhCjavcAL+WoJVd3DwCRJEKIwqUkIIYRwq8Ld3aSUagX8HWeD3cda6+0+DkkIIfxWRbzcdD9wHOfYMgd9G4oQQvg3n9ckXAOLDXY9XY1zyIEXgU7ALfzV2ahIdrvdYbP57yUzk8mAP7/+4kjZFE3KxT1/KhuLxXQWiC5pO58nCa31JGDShedKqc+BdJxj8pRY07HZHCQlZZRfgBVcRESIX7/+4kjZFE3KxT1/Kpvo6PBDnmxXrknCNVDaG65hBYw4awXtgWyc967vK2K3T3COwJlLyT1JhRBClKNySxJKqSeBm3HWCgDGAkFa6x5Kqe7AOxTR21VrvR7nZSaPmEwGIiJCyiDiyslkMvr16y+OlE3RpFzck7K5WHnWJPYD43GOvAnQG5gHoLWOV0p1KYuTyOUm/6kel5aUTdGkXNzzp7KJjg73aLtySxJa6xmuiVMuqEbBWbtsSimz1tqKEH7KZrOSmHgGq9V7wyWdOmVA+kcVrSqWjdkcQGRkNCbTpX3de7PhOgXIn7qMkiCEv0tMPENQUAihobUxGLzTudtkMmKzyUCvRalqZeNwOEhPTyEx8QxRUXUu6Rje7CexEhgJ4GqT2ObFcwtRIVmtOYSGVvNaghD+xWAwEBpa7bJqqt6sSfwCDFFKrcI5GmVRM4p5n8MB8gEVPiQJQpSny31/lWuScE0G09312I5z/PsKw7xhHWEvPUfqOx9ga658HY4QQlQ4FXFYDu/IySH80YewrFlN5MBeBM6c4euIhBCiwvHbJGFZtwbTfldfPrOZ3I6dfRuQEEJUQD4flsNXcnv1IXHhMsIffZDsa67F3qChr0MSQlRxmZmZvPPO61gsFjp27MzQoSPy1p08eZL33nuD6tUjiI1twDXXXFvktpmZmTz44N3ceee92Gw2Vq9eQWLiecaPn8gVV3Qv85j9tiYBYGvVmqS5i8m8456L1ln+WIRl9UofRCWEdx08mMCECVdhtztv/bTb7Tz66IP8/vucUh8rPn4V118/nmuvHcs333xd6m3crXO3PDs7m7vvvoVbb72em26ayBdffFrgeDabjdtvv4Enn3wEgMOHD3LbbTfk/Q0d2o9p06aW+nVe8Msv03n77dfznn/22ce88srzbrdfuvQP+vcfxFNPPceKFcsKrDty5BA9evTm2Wdf5ODBA263/e67/zFwoHNM1L59+/PUU8/x7LMvsXjxgkt+HcXx25pEHpPpokWGpETCH74f0+lTZN52J+kvvIwjzLPeiUJUNg0bNiIuriGrVi2nd+9+fPrpR8TFNWDEiNGlOo7NZuPdd9/gvfc+IiamFnfddQu9e/elUaPGHm3jbl1cXAO3+wQEBPD++58QEhKC1Wrl/vvvpFu3nrRp0xaAn376ngYNGpGR4RwdKC6uIV9/PTUvlnHjRtK374CLXkti4nkCAwMJCQnNW3b06BHq148tsN2IEaO54YarSU19kK1bN7N69Qo++eRL9u/fx6efTi6w7TPPvMCZM6dp0qQpAEZjwd/ozZsrpkz5kj/+WMiwYSOL3HbdujU0bNiInJyCt7T+739fMH78RA//p0rHb5NEdk4OSalpREVGYCr0nxX62suYTjtncQz4/TfS//FiUYcQosq49tob+PHHqVitVrZt28IHH3xS6mPs2rWD+vVjqVevPgCDBw9lxYqlBZJEcdu4W9exY2e3+xgMBkJCnGMtWa1WbDZr3i2fp0+fYvXqldxyyx38+ON3F8W7YcM66tWrR+3aF3cy27RpIzNnTuett94nMDCQWbN+YdmyJbz99gcFtgsKCmLw4KF89tnHxMevYtKkjwgMDKJJk6a8+eaki44bHR3D6dOnadZM4XAU7LT322+zufPOe+nQoRPPPfckvXv3u2jbjRvXk5WVxcGDBwgMDKR79558+ulHdO/eE6VaePC/VHp+mSTsdjtTZs4jJS2dqMjqDOjWidg6tfLeXBmPPo7xxDECF8wj7Y13cVSP8HHEQpSvK67ozuTJ7/HJJx8xefJnmM0FvxoeeOAuMjIuHtPowQf/Tteu3QA4c+Y0MTG18tZFR8ewc2fBiSWL28bdupKOa7PZuPPOmzl27Ajjxk2gdes2AHzwwTvcf//DebWIwhYtms/gwcOKXDdw4GBOnDjGSy89y4ABg/ntt1m8995HRW47atQYbrzxGl5//Z28ROZOv34Dee+9N1m9egW9evUFYPv2rezZo+nevQdffvkZCxfOo3btukVue++9DwIwd+5sqleP4Oeff2L9+rWkp6dx7NgRxo69ptjzXwq/TBLZubmkZ2QCcDYxmZ/mLaF6WCixdWvRubUiqm49Ur75Ecua1eR273nR/sYD+7E3aiyd8ESZC3nzX4S6rnGnP/40GU8+W2B96AvPEvKJ8zJG2kuvkfnA3wqsD/u/hwl2XbNPfft9sm7xvM9qmzbtaN5cERUVddG6jz/+vMT9ixryqHBHruK2cbeupOOaTCa+/noqqampPPvs4xw4sI8TJ04QEVGDFi1asnHj+ov2z83NZeXKZdx330NuX8+NN97Kiy8+w9tvv86PP87Mq7EU9tVX/yUiIhKbzeb2WBcEBwfz7LMFr0y0adOONm3aAfDqq28WWFd42wtGjrwy7/GECdeVeN7L4ZdJIjgwkCG9rmDR6nVYrc7/2OS0dJL3HGDH3gS6tm1Jr05ti04QBxOoMbAXOX37k/bme9iLqKoKURkdPJjAyJFXFbnOk5pETEwMp12XacFZM4iKKjjxWXHbuFvnyXEBwsPD6dixM/Hxq0lJSWblymXEx68kJyeH9PQ0Xn75eV544RUA4uNX0rx5C2rUqOm2PLZs2cSBA/vp27c/X331GY899tRF23z//bfk5OTw8sv/5ssvP6Nfv4Fuj1dZ+WWSAGjdrBH1a0ezdusudh84SE6uc6xBh8PB2q072XvwCIEBARiMBpo1qE/n1gqjwUD4//0dQ0YGgfPmYkxMJGn2fB+/EiHKRkLCARo3blLkOk9qEi1atOLIkSMcP36M6OgYFi1awIsvvurxNu7WxcU1cLtPYmIiZrOZ8PBwsrOzWL9+LTfeeCu9evXJqyVs3LieH374Ni9BQPGXmgD27NnNG2+8yltvvU+dOnV5+eXn+eyzj7nnngfyttmwYR1z587i00+/IiQklA8/fJe9ezXNmlWt0Rv8NkkAVA8PY0ivrgzo3okTp8+yatN2jp48DUBiSmredidOn2X73gM0qRVF39p1iAYcRiNpL//LR5GLqirjyWcvusSUX/rL/yK9mPdd2jsfkPbOB27Xu3Pq1EnCwsLdXlLxhNls5rHHnuCxx/6G3W5j1Kir8pLO448/zNNPP09UVLTbbYrb393yc+fO8tprL2K327Hb7QwcOIRevfoUG2dWVhbr1q3liSf+Uew2r7zyRl4bw3PP/ZO5c2fnrT958iSvv/4qb701Ke8OqAkTrmfatO/5xz9eurQCrKAMlX3s9Nxcm6OsJglxOBys3ryd+M07ih1TvvXxQ3TNTafa8y9iKdTA5+0BA/1pkpTSqgxlc/LkIWrXbuDVc1a14bDLUlUtm6LeZ9HR4RuAEid/8+uaRGEGg4GeHdvSvkVTklLSMBoMHDh6nPjNOwpst6NuA3YAYTN+o3v7VtSvFUN4WAghWzYR9sKzzgEDW7byzYsQQogyJEmiCKHBwYQGBwNQJyaKts2bcPLMOY6cPI1OOExmVjYAaekZLFrlvHMi2ABPfvoWlgP7iBzch9QPPyF7/ASfvQYhhCgLkiQ8UC0slGphoTRvFEffrh3YvucA8Vt2kJGZlbdN7f17sBxKAMBmMpPTuauvwhVCiDLj12M3XQqL2UzHVs2585rR9OzYlpoR1QBIaNycDx5+jsOxjZg75Cp+3X+Uc0kpPo5WVAaVvV1QVGyX+/6ShusykJSSxqzFyzmTmITBNUiawzXUR6P6deneoTUNtm+CwCByS7jzorQqQ+Osr1SGsjl79oRrjmvvTWFaVRtny0JVK5sLc1xnZWVcNMe1pw3XkiTKiMPhIDk1nRUbtqATDhdYF5yRzhPvv0xIchKZN99O+kuv4AivVibnrQxfhL5SGcrGZrOSmHjmsuYgLi1nL+bK/bkvL1WxbMzmACIjozGZCrYuyN1NXmYwGIioFsao/j1p1bQh2/QB9h0+CsDQBb8SkpwEgH32TBLuuI+4lmEXjQIp/I/JZL7oF155qwzJ01ekbC4mSaKMGQwGGsfWo3FsPc4lJbNmy06WDBpFtZRkWu/czIwrr2Pnms1U37GX4X27U792jK9DFkIIt+RykxckpqSyZvN20hf/wcEGjQusaxxblz7Vg4nq2vWSOuHJLx/3pGyKJuXinj+VjaeXm+R6hxdEVgtneN8eDH7873Rp24LAAEveuqTNW2gybgS5o4dzZuu2Knc9VAhRuUmS8KLq4WH069qRW8eNpHFsPXA4uHrGNwTk5lJ33WqC77qVWYtXkOHqrCeEEL4mbRI+EB4awrghfUk8fYb0hb/CAY3dYGDO6IkcO3yUI9NP0bGVQjWKJSpSJjwSQviOJAkfioyJxj71e7bMnEnmqlUcq+8cgCs7J5f4zduJ37ydji2b0b9bJ7kTSgjhE/LN42NGo5G648fT5O23GTekLxHVwvLWxR3aT4+H7mDBfz5nz8EjHs18JYQQZUlqEhVI49h6NKhXh70Hj7B92y7Gv/sNtU6f4IZXnmSa3sWi7n3o1akd7VQTr/XOFUL4N6lJVDAmo5EWjRtwXXQYUYlnAbCaLRxq0ITMrGwWrVrHT7//UWBSJCGEKC+SJCooa8/eJP+5mqzOXdl11wNY68fmrTty8jRTfvmd1Zu2c+psog+jFEJUddKZrqJzDRiYa7cTv3k767btxuFwoHZvw2Yys69ZSxrH1mP0gJ4Xz5Ln5/ypY1RpSLm4509lI2M3VRWuu5osRiN9unSgecM4ls5fzDXT/0d4WirruvRizugJ/GqzMbJfD0KCg3wcsBCiKpHLTZVMraga3L55FeFpzjaJ5nucU6seOn6SL2fMYcGKtZxPlnkshBBlQ5JEJZTx+NNkXzkWgP1PP092kHOq1eycXLbt2c/U2Qs4fOKUL0MUQlQRcrmpEnLUqkXKF1Mwb1hH20H9sG3fz+L49SSlpAFQ7egRfpq7mGrhYVzRrhVtmzeWznhCiEsiDdeV3IWGNofDwYEjx1n786/c99bzJDRqzs/jbyI5ogY1I6rT/4qONKzv3XkLfM2fGiFLQ8rFPX8qGxkF1s8YDAaaxNblzkW/EpCbi9qzg+t++AKAc0nJzFjwJ78uXk5KWrqPIxVCVCaSJKqS3FwMHTrgMBhwGI0cefwfBW6L3XfoKF/9/BsrNmwlIzPLh4EKISoLudxUyRVVPTavXYNly0Yy776f9MxMlq/fwo69Cc6VDgcYDJhNJlo3a0SXNi0LjBdVlfjTpYPSkHJxz5/KxtPLTZIkKjlP39THTp1h5zdT6TP1c34efzMn6jp7cBsMBjq2ak7frh0wVbHGbX/6wJeGlIt7/lQ20iYhCqgXUY2Jv3xL7NFD/G3yv+iwKR4Ah8PBxh2aqbMXsGv/QeyuHt5CCAGSJPyGZdMGTEcOOZ8Eh9Dq1luIrR2Tt/70uUTmLl3Nd7MXcOzUGR9FKYSoaCRJ+Inc7j1JXLKSnO49SX/uJep06cyEEQPp06V9gT4Up88l8sNvi/h92WoyZRpVIfyetElUcqW+hnrhclK+xJCekcmJ//0PfeQ4u5u2zFseHBRIt3ataNeiaaUcPNCfri+XhpSLe/5UNtImIYpmNBZIEABhWRl0+/Atbvt8EnfO/YmgTOeHJDMrmz/XbuLbWfM5fU6GJBfCH5X481ApFQn0AmoCp4EVWmuZ8aYKCX39VYxnne0QjffuZNTL/2L+Nk1aujNZnE9K4Ztf59Ekrh5d2rSkXq0omRlPCD/h9nKTUioaeB1oAWjgBBAJdAS2Ac9rrct8FDml1CNAB6AZ8J3W+uPitpfLTZdfPTacOUPYc08S9MsMkr/5kZxhI7BabWzV+1i+YQtWa8G5tZs1jGVorysICgy4rPOWN3+6dFAaUi7u+VPZXHY/CaXUZOADrfWeIta1BB7SWj94uYG6OXcAMBm4T2td7D2ZkiTK7k1t3rwRa4dOBZadT05h68+/ssFugny1h2phoYwe0Is60TXL5NzlwZ8+8KUh5eKeP5VNpelM56o5DHY9Xa21fk0pdStwSms9r6T9JUmU75vaeGA/Nfr3IKNzV+bdcBfrUv8azsNoMNC5TQuuaNeqQtYq/OkDXxpSLu75U9mU2cx0SqlewEdAbeAocLfWetNlR+iitZ4ETCq0uA9wd1mdQ1wiu53wx/6GISuL0JXLuSo3l9qTv2DBijVk5+RidzhYt20Xm3buoVXTRvTq1FZmxhOiivHkvsYPgBu01juVUm2Az4CenhxcKdUNeENr3V8pZQQ+BtoD2cBdWut97uLSWlfue3OrAqsVa7sOWOJXgcFA2r/epHnDWGJqRDJveXxepzurzdl+oRMO0U41pUvbloQEBfo4eCFEWSjxcpNSarHWelC+54u01oOL28e13ZPAzUC61rq7Umo8cJXW+jalVHfgGa31mMuMH7vd7rDZ/DefmExGbLbyHUrDsHYNhvXrsT/wVxOUw+Fg6+79rFi3lZNnC94eWz08lBvHDvF5e4U3yqYyknJxz5/KxmIxlc3lJuC0Uupz4A+gM2BUSt0DoLX+rJj99gPjgW9cz3sD81z7xSulSgzOEzabw2+uIRbFK9dQm7d1/hU6T5MjR2j//qvsfuoF5qbkkJjivDM6OTWdT779lbaqKX26tCcwwFK+8bnhT9eXS0PKxT1/Kpvo6HCPtvOkM91u4AjQFEgBlgJ1cLZRuKW1ngHk5ltUDUjO99ymlKp83XiFU1YW4Y8+iGXbVtrcei335SYxqn9PAizO/1K7w8GW3Xv59td57Ek4jM1mK+GAQoiKqMQvaa31P5VSVwLNge1a6/mXeK4UIH/qMmqtrZd4LOFj5i2bMR09AoAjKBhrn360qB9LdGQEi1av4+hJZ3tFUmoas5esJCw0hF6d2tK6aSPpiCdEJVJiTUIp9RFwI2AF7lJKvXmJ51oJjHQdszvODnmikrJ26875P1eT07sv6c//E3t95/wUNSOrM3HEIEb261FgvKe09AzmL1/DvOXxUqsQohLx5HJPe611b9fj95VSqy7xXL8AQ1z7G4DbL/E4ooKwN2pM8ozZztnu8jEYDLTft4tmkYGsqavYqveR4RpRdue+g5xLSmFwzy7Ujqq4HfGEEE6eJInDSqn6WuujSqlaONsnPKK1Pgh0dz22A/ddUpSi4jIYCvTEBjCcO0f4ow9S/exZhl89kSte/jd/7D7A9r0HADh19jzfzVpAneiaDOvTnZoR1XwRuRDCA24vNymlTiiljgPDgb1Kqd3AAVxf+kK4E/rGqxjPngXAsnolAUGBDO19xUVTpJ44c47v5yzgyIkyHwJMCFFG3NYktNZ1vBmIqDrSn3oOQ2oqQTOmkfbWezjCq2EAurZtSePYusRv3sGehMPYHQ6yc3L5ecFShvfpTvNGsdKoLUQF4/Oxmy6XjN1Uce/rNm/djLVdh4uWm3bu4Hh0bX5ZtJz0zMy85ZHVwundpT3NGtQvk2RRkcvGl6Rc3POnspFJh4TPFZkgDuwjcvgAmt97Kze1aUJYaEjeusSUVGb/sYI5S1aSk5t70b5CCO/z5BbYx70RiPADdjthjzoHDAxYsYz6Tz/GDaOH0L5F0wK9svccPMJ3sxdw4vRZHwYrhADPahIjlVKmco9EVH1WK9bOXXEYjThMJtL+/RbhoSEM7tmVuyZcRTvVNG/T80kpfP/bIrbq/T4MWAjhyS2wUcBxpVQC4AAcWmuPRoEVooCAANJfeJnsq8Zi3rihwOWooMAAhvTqSp2oSP5Ys4lcqxWHw8HClWs5eOwE3dq1olZUDR8GL4R/8iRJXFnuUQi/Yu3Q6aIZ8ADMa+Lp89RjNH31Daafy+DUufMA7D14hH2HjtK3S3s6tVYYjdKUJoS3ePJpswKv4ZxHYiglDOwnxCVxDRho3rmd+tdcya2pJ2kcWzdvtcPhYOm6zXw7az4nz57zYaBC+BdPksRnwJdAALAMeL9cIxJ+ybxlM6bjxwBwhITiGDCIcUP6cfOY4QXmpThzPokf5ixi+54DvgpVCL/iSZII0lr/gbMtQgNZJe0gRGlZu3Xn/NJ4cvoOIP2Fl7HXqw9ATM1Irh05iB4d22A2O++fsNntzF+xhk0791DZ+/kIUdF50iaRrZQaBphco7dKkhDlwt6gIck/zbxowECTyUT/k4fpGm7ie8I5k5gEwB/xG1i/fTc9OraRIciFKCee1CTuwTliaxTwOHB/uUYk/JvBAIUapg3nzhH+2EPUvfd27p/7A3FBf/WpSElLZ/7yNcxespLcXJmeRIiyVmKS0FofBd4AXgH+qbVOKPeohMgn9I1XMZ5zNlYHrlvL6CH9adW0YYH5KvYePMK03xeTkpbuqzCFqJI86XH9KjAZ6AZ8qpR6otyjEiKf9KefI2vi9QCkvT2J4OgoRvTtwX3Xj6V9i7864J08e55vZs7j4NETvgpViCrHk8tNw4E+WutHgT7AhPINSYiCHDVqkjr5U87/sZKcwcPylgdYLAzu2ZUxUeEYXe0YWTk5/LxwKRt3aOx2u69CFqLK8CRJHOWvuaktgAz+L3zC1qbtRctM+/fS/aE7eHbGV8SmpwDOPhVL1mzkvz/MITXdP0b0FKK8eJIk6gJ7lFILgH1Aa6XUqsuYxlSIsmG3E/bYwxiysghbv4Y7504jpkZE3uqjJ88wbe5iTp9L9GGQQlRuntwCK5eXRMVks2Ht2g3L2ngAMt98l+tatmbd1l2s2bIDu8NBUmoa3/w6j7bNmzCoR2dMJhmrUojSkEmHKjl/miTFHfPWzZg3biDrtjvzlu07fJQ5S1Ziz7XicN1SG1e3FqP79yI4KNBXoVYI8p5xz5/KRiYdEn7D2q5DgQQB0DSuPg83ieGJj18n9rBzCI/Dx08xZebvHJY5tYXwWKmShFJKkoqoHLKyiHnqCWocPcQD/3mTLutWApCWkclPv//Byo3bZEgPITzgST+JCUqp65RStwInZaY6URmYt22BCwMGhobR+I7bC1xmit+8nUWr1pFrlV7aQhTHk5rB48BC4CYgFplfQlQC1q7dsG7eSs6AQaS/+AqxXTtzy9gRxNWtlbfNVr2fr2b8Jp3vhCiGJ0niwoB+qVrrbP7qMyFExdZ/R5bZAAAgAElEQVSgAck//EzWLbcDEBYSzPih/WnRuAGtdmym9fZNpKZnMGPBn6zetF0uPwlRBE9ugU0A1gN/U0q9CKwp35CEKEOFRoY1GY2Mat2M6nd/T0ByElvbdmbm2BtYtWkb2bm59OvaQUaTFSKfEpOE1vo2pVSY1jpNKbUeyPZCXEKUm7C3/01AsnO48UbHDmF1DRS4YftuzpxPZETf7oSFhPgyRCEqDE8arj/UWqe5nrZHahKikkt/5nkyr78JgNwPPiKu+V+DBB4+for//fI7+w4d9VV4QlQonlxuSlFKvQ6EAW1wDvgnRKXliIgk7f2Pybz3QWytWnOl3c7KjdtYu3UnADUO7GNWZhbtWjWnf7dOmKWXtvBjnswn8Q/ABDTVWvfXWsvkwqJKsLVqDYDRaKRPl/ZMHDGIuPRk7v/PG9zz2TscW76C6fOWkJEpkzEK/+W2JqGUOgFcuN3DANRSSh0H0FrX9UJsQnhVbK0oWs3/GYvVSqOD+7hm+hQ+iXqCqbMXMG5IP2pGVvd1iEJ4ndskobWu481AhPA5mw1b7z44Nq7H4XDw65jrwWAgOS2dqXMWcuWAXjSsLx8L4V9KHOBPKdUB5zzXQReWaa3vKOe4PCYD/PnPgGSldallY9q+DcvG9WzvP4S5f67O65VtMBgY2LUDHdq0KOtQvUreM+75U9mU5QB/XwMbgR/z/QlRZdnatCXrlttpGlef60YNJizUeTtsgwN7aHfdGFZ/9AnJqWklHEWIqsGTu5tOaq0/L/dIhKiAYmpGcuOVQ5kzdyFXz5hC9NnTjH75KWbt3En044/TvFGcr0MUolx5kiQOKqWeBjbhasjWWi8o16iEqEDCQoK5Ljaa6unO2kNOQCA7m7YgZclKhmTn0K5F0xKOIETl5UmSCASU6w+ciUKShPAv3XqQumo9locfYGW9xqRUjwRg4ap1nElMokeHNoQEB5VwECEqn+JugTVrra3AvV6MR4gKy16vPtnTf6VNdg67F/zJqbPnAdi8ay+m2b/Ss11LAq67wcdRClG2imu4nuL6VwO7XX8XHgvhnwwGgoICuXpof+rVigYgNC2VkVM/p97D92G5YQKGc+d8HKQQZae4fhI3uP5t5L1whKgcgoMCuXbkIPYcPILl0b8RmpHuXLF+HbuOHKdFzZq+DVCIMlLc5abV/NXjugCtdc9yi0iISsJgMKAaxXHq7ffY8H9/p/Pqpfw87kb2bNzOscxs+nRpT4DF4uswhbgsxTVcX+e1KISoxGo1bUzyt9/z1ZSp7Al2zsm1edde9h8+xsh+PWh46hjWNu1ABgoUlVCJPa4rOulx7T89REvL22WTk5vL3KWr2X/4WN6y6LOneOT9V7G170DapI+wNVfFHME75D3jnj+VTVn2uBZCeCDAYmHMoD6M7NeDoMAADHY746dPwZSbQ8D6tYQ98iBU8h9lwv+4TRJKqddc/47xXjhCVG4Gg4GWTRpyy9gRxMbUZF/TFlhNJmxGI9NHX0t6lgw7LiqX4tokxriGBv+bUqpW/hVa68/KKyClVFfgdqAa8I7WelN5nUuI8hIeGsLVo4eysm4dPmzTidgjCWwNCmP/zPmMHtCT+rVjnBtarWD2pE+rEL5R3OWm24EInD2u6xT6K0+dgVZAfeBIOZ9LiHJzYTKjlmNGs75rbwDSMzOZ9vsfrNmyE/OqFUT2uQLLqhU+jlQI9zwZKrwrsB9oAiRorc+WZQBKqUeAwa6nq4GFwGZgIFBLa/2/4vaXhmv/aWgrrYpUNgePnWDu0tVkZmUDYM7N4fGP/k3EyeMApL79Plm33O6VWCpSuVQ0/lQ2Zdlw3QiIB/4BxCulbrrM2ArQWk/SWo92/b0GPAbYgTOA9EgSVULDenW4ecxw6sZEAVDnxDGCEp3DetjDw8kZMsyX4QnhlicXQx8FOmmt05RS4cAfwLeeHFwp1Q14Q2vdXyllBD4G2gPZwF1a631F7PYD8A2QAzzhyXmEqAzCQ0OYOHIQK9ZvYT3w7mMvMXbm9+xr04GIjBxa+jpAIYrgSZKwa63TALTWqUopj27PUEo9CdwMuMYrYCwQpLXuoZTqDrwDXHTnlNZ6JjDTk3MAmEwGIiJCPN28yjGZjH79+otTUctmzNBe1KlVgzmLVzHllvudC5euJjUjnUG9OmMyGjH8/DPk5uKYOBEMhjI9f0Utl4pAyuZiniSJ/Uqpd4BlQF+c7ROe2A+Mx1krAOgNzAPQWscrpUq8FuYJm83hN9cQi+JP11BLqyKXTfMGDbhmeCALVq7Lm+Vu+bqt7Dt4jDFtFXEP3Ifx/Hmyv/uO1Hcn44iKKrNzV+Ry8TV/Kpvo6HCPtvOkTeIO4AAwxPXv3Z4cWGs9A8jNt6gakJzvuU0pJff+Cb8VV7c2N48ZTqP6f90weOLMOc49/hjG8872CvP2bRAU6KsQhSi5JuGaU+KjMjhXCpA/dRldxxbCbwUGWBg7uC/rtu1m1aZt2O12Zg0cBQ47nZcvJvXtSTjCPPvFJ0R58OawHCuBkQCuNoltXjy3EBWW0WikW/tWXDWwN0ajkeygYH4aNZFJT77C4vBocq0Ff0uZN20Am81H0Qp/U2KSUEqF5P/3MvwCZCmlVgHv4bxrSgjh0iSuHmMG9cFsdo4We7JGDOu27WL6vCVkZecAYNq7h4irhhMxajCmXTt9Ga7wE2470yml/gWsBf6utR6glPpDaz3Qq9F5QDrT+U9DW2lV1rI5dfY885ev4UxiUt6yyOrhXNmvJ81uuw7L2ngAcrt2I2nOglLf/VRZy8Ub/KlsyqIz3QygNtBWKbUQaKGUesh1qUgIUU5qRdXg5rHD6de1Q96yxORUfpg1nyPtOuEICMBhNpP65ntlfnusEIUVlyRCcXaa26q1HgIkAIcB6RoqRDkzGAx0aduSkf165F1+ygE+a9qeaS+/y+lX38DWus3FO1rlXhBRtopLEj2B/wHtlFL/wTnY3y6t9T+9EpkQgpZNGnLTVcOoGVE9b9kmLHwcEMmBI8cKbGtZuZzIXl2wrFjm7TBFFeY2SWitX9daX43zLqQ3AQPwsFJqnreCE0JAzYjq3HjVUDq0bJa3LDsnl18WLmPesnhycnMhI4Owx/6GOeEAEeNHE/T1Fz6MWFQlxU06dOFupv/TWicAN2mt/6a1Hl5ovRCinFnMZgb16ML1o4cQFvrXR2/HvgSm/f4H1i2b8zrg2atVJ2f4SF+FKqqY4i43faSUehA4BKC13giglIpyDe/9Hy/EJ4TIp25MFDdfNYzmjeLylp06e55PD5xg5Vc/kjVyNOn/fA177fKe9kX4i2Lnk1BKTQT+BsQCZ3EOrXEC+Fhr/aNXIiyB3ALrP7fslVZVL5stu/exaNW6Ass6tmzGgG6dMBgL/v4LmD0TQ3Y22VdPJCIytEqXy+Wo6u+Z/Dy9BbbESYcAlFJBQCRwTmudc/nhlR1JEv7zpi4tfyibhCPHWRy/IW+QQHB2yhvZrwcBFgsAhjNnqNG7C8bERLKHDMP41VckBYT5KuQKzR/eMxeU5aRDaK2ztNYnKlqCEMLfNYqty23jRtK8YWzesv2Hj/H9nIV5iSP0ndcxJiYCYN61E0KkOVF4zptjNwkhyoHZbGJU/550btMib9nZxGS+nTWfQ8dPkv7cS2TecTcOg4HUdz6AMKlFCM95MnZTZ28EIoS4dEajkf5XdGR4n26YXO0RWdk5TJ+3hK8WLGf9PQ+TuGo9uQMGXbSvef1a6YQn3PKkJvGEUireNSRHRLlHJIS4ZK2bNWbiyEGEBgflLTuXlMxvf65i4ZkUCrdBmvZoIsaOJGLEIEzbZWBmcbESk4TW+jpgBOAAflJKfaeU6l/egQkhLk3dmChuvGoYjWPrYcw3ttOmnXuYv3wNVqtrmHGbjfBHHsSQk4NlyybCn3oMPLiRRfgXT9skagFxQBTOW2EnKKW+LLeohBCXJTw0hHFD+nLf9eNo1qB+3vId+xJ478uf2LE3AYfdTvawETgCA50DBr41SQYMFBcpcWY6pdQaIAP4L/CC1jrbtXx+OccmhLhMwUGBXDmwN/NXrGHH3gQAUtLSmbc8Hp1Ql1H3P0zOyCsxb1iHrVXriw+QmwuuW2mFf/KkJvGQ1nqA1nqq1jpbKdUPQGsto8EKUQkYDAaG9e7GgG6dCMnXVpFw9Dg/zl3EqahaZF9340X7WVYso0bPzliWLvFmuKKCKW7SoT5AK5wzyL3rWmwCHtRaFzFGsW9IZzr/6fxTWlI2F8vJzWX99l2s3rQjb5nBYKCdakrPTm0JCQp0LszIoEb/HpgOOmsfqW+8S9btd/kiZK/yp/eMp53pirvclIhz0qFA4MJAMHbgycuOTgjhEwEWCyMHdKdaaBgLV67F7nDgcDjYsnsvuw8cpGfHtrRv2YxAvQtDkrMDnr1adXJGjvZx5MJXShyWQylVR2t9wkvxlJrUJPznl09pSdkU7UK5nDmfxJ9rNnL4xKkC6yOrhzOsdzfqGxyEP/M4OYOGkHXDzT6K1rv86T1z2WM3KaWma62vUUqdwHn7KzjnlHBoreuWWaSXSZKE/7ypS0vKpmj5y8XhcLD/yDGWrt1EUspf4z8ZjUaG9OxK62aNMMBFdz0FzPoFQ2Ym2ROvr1J3RPnTe6ZMB/iryCRJ+M+burSkbIpWVLlYbTY279zD6s07nJMYuTSNq8+gnl0ICwnOW2Y4fZoafbpiTEwkZ+BgUj74BEdMjNfiL0/+9J4pizYJAJRSfYEQnHdCfQg8r7WeetkRCiEqDLPJRJe2LWneKI5fFi7lbGIyAPsOH+XoydN0aq1oq5oQFhJM6Htv5g0YaNq3F4cMGFileXIL7JvAXuBhoBdwX7lGJITwmWphoVw/egjtVNO8ZVk5OazatI0vfprNivVbSH7meTLuuR+H0Ujq2+/LgIFVnCdJIhM4BVi11idx3u0khKiiAiwWhvTqyjXDB1A9LDRvudVmY83WnUxbtpbT/3iJ86s3ktt/4EX7m9eucXbCE1WCJ0kiBVgETHNNZ3q4fEMSQlQEDerW5rarRzGqf09iakbmLT9++ixTZs5jp8OEzWYrsI9J7yZi/Cgihg/EvG2Lt0MW5cCTJDERuENrPQVYClzcNVMIUSWZTSZaNG7ATVcNo2/XDnnLMzKz+O3PVfx32ix0wmHn6LI2G+GPPuQcMHDbFsKeflwGDKwCSmy4BmKA0Uqpa/Ite7mc4hFCVEAGg4GubVsSXSOC35fFk5GZBUB6ZhZzlqzEYjajYusyYvAwamzfCjabs72iCt0e6688SRI/4bzcdKScYxFCVHAN69Xh1rEj2LhzD9v37CfdlSxyrVa2JxxGRzVk4H+m0DH1HLaWrS4+QE4OBAR4OWpxOTxJEqla6+fKPRIhRKUQEhxE787t6NKmBcvWbWbPwcNk5zgbqnOtVuafSiQ+vBrd9x6gWYNYAgOco8hali8l/JEHSX1rErkDB/vyJYhS8GRYjveANcAmXD2vtdZ7yj80z0hnOv/p/FNaUjZFK+tycTgcHD5+ij/XbszrX3FBgMXCoB6daVW3FjX6dcd06CAAqa+/Q9Ydd5dZDGXFn94zZdaZDujg+rvAAVx835sQwi8ZDAYa1KvNzWOGs2X3PlZt3EZWTg7gHHX292XxJCSd4YakJEy4BgwcdaVvgxYe82hYDqVUdaABcEBrnVbS9t4kNQn/+eVTWlI2RSvvcsnMymbbnv1s2b2PlLT0vOWhaSlcv3g29oGDCb7vAYICK17bhD+9Z8ps7Cal1NXAczhrHdNwDvD3alkEWRYkSfjPm7q0pGyK5q1yycnNZenazezcl4A1f38KhwOT2UTLJg3p1aktYSEhBM6cAZmZzsmPfHhHlD+9ZzxNEp70k3gM6I5zbutXgXGXF5oQwh9c6Ll93/VjaRJX768VBgM2m53tew7w5fTf2LVqDWFPPkq1vz9A9YljMZw65f6gwus8SRJ217zWDq21A0gvaQchhLggMCCAqwb2ZmS/HnRqrQr03s61WjG9+S+MSUkAGA/sxyFjQVUonjRcL1dKfQ/UV0p9Aqwr55iEEFWM0WikZZOGtGzSEICDx06wJH4j55NTmDd8PHajiV4rF/PtyAlE7N5P13YtCbBYfBu0ADxvuB4OtAV2a61nl3tUpSBtEv5zDbW0pGyKVlHKJddqZfWm7Wzfe4DMrGwiz58lsUYUABHVwhjRtwd1Y6Iwx6/G2rkLeCFpVJSy8YYyuQVWKdUOuAaIAo4Cu8skOiGE37OYzfTt2oHendux79BRVm3aBkkpACSlpPH9nIV0JJcJzz+GvVlzUt//CGv7jj6O2v+4bZNQSk0AvsQ5HMc8IBWYrpQa46XYhBB+wGg00rxRHLeMHcGw3t2wmJ2/XQ12O90/fgdjbg7mndsJeUYGDPSF4moSfwf6aa3zGqqVUl8Ds4BfyzkuIYSfMRqNtGnemPq1o1m5cRt79h9kW9su1DlxFKPdzucDx9B89z7aqSYYjZ7ccyPKQnFJwpo/QQBorVOVUjZ3OwghxOWKqBbOqP49uaJdK/6sW5tJrdoTd/gAR2pEc2T1erbs3suAbp2Iq1sbsrMhUOZBK0/FpWP7JewjhBBlIrpGBNcMH0DvCeNI6PPXSEBnE5P5ad4SNk7+iMhuHQhYNN+HUVZ9xdUkWiulphZaZgCKGP9XCCHKnsFgoFnDWBrVr8uGHbtZs2UnuVYrlpxsukx+B/P5s1S/YQKJL/8b630P+jrcKqm4JDHRzfJPyiMQIYRwx2w20a19a1o3a8SydZtJXracwOxsADKDQvjcHkirrTvp0qaFtFeUMbdJQmu91JuBCCFEScJCQhjZryc76tZmclQMQ2dOZX8TxbmgUJav34JOOEyPDm1oElcPg8yKVyY86XHtVUqp64DRwDngOa11qo9DEkJUMK2bNaZZw1j0oL4c2LYLkp1fE6fPJfLr4uX02rud1rWjCX/47zKF6mWqiPWyq4DbgCnALb4NRQhRUQVYLLRt3oTbxo2kd+f2mEzOr7Ow1GQGTf2Cxq+9gGHoABzHj/k40srN5zUJpdQjwIW5DFcDHwL/xdmJT263FUIUy2g00q19K1o3a8SGHbup9++XCcl0Dq1hOHKYrxetYNDgfs5bZkWp+bwmobWepLUe7fp7DagD3AUsAw77NjohRGURFhJMv64dqfflV2wYOQ6b0cjP42/ivNXO9Pl/smX3Pl+HWCmVa01CKdUNeENr3V8pZQQ+BtoD2cBdWuui/teSgK9x1iLuLc/4hBBVT3BEBPW++JKlS5Zy/HQSZGXjcDhYtGodZ84nMiAnFVPP3hBQ8WbGq4g8GgX2UiilngRuBtK11t2VUuOBq7TWtymlugPPaK0vexwou93usNn8dzwXk8mIzeau36N/k7Ipmj+VS0pqOt/9upDjp88BEHPqOA9/8Co5DRsRMGUKjq5XFNjen8rGYjFd/iiwl2k/MB74xvW8N86BAtFaxyulSgzOEzabw2+G9i2KPw1tXFpSNkXzr3IxcM2wgfy+PJ59Bw5xzfQpmG02zPv3cfaue8hdtAyz5a+vQX8qm+jocI+2K7c2Ca31DCA336JqQHK+5zallM8bzoUQVZvFYubKAb0Y2a8He7r3IccSgNVkYsrwq/l29gJOnTvv6xArNG9+SacA+VOXUWtt9eL5hRB+ymAw0KJZY3ImvceMnn1hzWpO16oLSclMnbWAHh3bcEW7VpCV5etQKxxv3t20EhgJ4GqT2ObFcwshBAEWC4NumEDU/Q9gNpsAsDscrNy4jZVvvYtRNSNg3lwfR1mxeDNJ/AJkKaVWAe8Bj3rx3EIIAThrFR1aNuOWsSOoE+OcLtWSk02/rz7GdOIE1W+5juDPPvZxlBVHuV5u0lofBLq7HtuB+8rzfEII4anIauFcN3IQ67btImHOXAJynAMGZgSHsKFVRxnu2kUajoUQfsvZW7s1jerX5fM6dek7bQr7myg27jnImYAg+nbt4PcDBUqSEEL4vZiakYy7fjyza8Vw7NRZANZv383pc4l0adsCFb8CY2oKWbfeAX42FLl/vVohhHAjNDiYOyaOoklcvbxlh0+cYuH0WQQ98QjhTz1G9XGjMJ484cMovU+ShBBCuARYzFw1sDedWqu8Zf3/nEdgehoAhmNHsYdX81V4PiGXm4QQIh+j0ciAbp1o3bQRa7fu5PcR48kJCKDvsgVMHT0RdS6JJqGhvg7TayRJCCFEEWJqRjKqf0/W1YxkQUAA8d37kRxRg12LltGptaJPl/aYTSYsK5aR27UbBAb6OuRyIUlCCCHcMBgMXNGuFXVjovlt6SpId47rtHGHZveBQ/QwOxj8yN3YGjUm9d3JWK/o5vZYuVYrx06dITMrm8CAAGx2O+EhwRgMBpJT0zAYDNSJqUlYSIi3Xp5HJEkIIUQJ6teO5pYxw5m3PJ4DR44DkJmegfrPGxhyczHv0RifeQLHwj8x5Lv7yWa3czYxiZ17E9i6Zz9Wa/HzqAUGWBg/tD8Oh4OziUnUiY4ipmZkub62kkiSEEIIDwQHBTJ2cF8279rL2q07SU9LZ2On7tQ6dQKTzcYXg8eQ8/NvBAcGkpqegcVsJjU9A6vN8wk2s3Ny+X7OwrznZpOJe64dQ3CQ7y5lSZIQQggPGQwGOrZqTvsWTUk4eoJtDWN5v1UH6h/c5xwwMCWNJNLytjfn5oDlr8mNIquFUzOiOlabFaPRyPmkFGx2O9XDwzh68vRF57PabJw5n+jTqVclSQghRCkZjUaaxNWjSVw9Mnt3Y922XZh37ilQa2i6dxcTp33FgutuJ2XQUFo1bUTj2Lpue3AfPXmGXxcvIys7p8DkR+mZvh2ZVpKEEEJchuCgQPp27UCPjm04ceYcublWalhMNPzwVQJSk7nmv5NIi6tF5sDexR6nfu1o7rl2DHa7nVWbtrNxhwYgLSPTGy/DLUkSQghRBixmM3F1agFg2rYVU24OAPaICLLGXuPxMQBCg4PylqVLkhBCiKrF1rYdiSvWEvrSc+T26IWjVq1S7R8WEpz3WGoSQghRBTkiIkmb9FGR6wJ/+gFDSgpZt99V5ICBocF/JQlpkxBCCD9iPHmCsGeewJiSTNDMGaR89hX2OnULbBOarybh68tNMsCfEEJ4UfDkSRhTkgEwnjqJvXrERduEhuRrk8iUJCGEEH4j/fmXSX/sSRwWC6nvfghFDMMRFBCAyeT8es7JtZKTm+vtMPNIkhBCCG8KDCTj6ec4v34bub37XrTasuxPDNnZBdslMnzXLiFJQgghfKBwOwSAaecOql9/NZEDetLsaELecl9ecpIkIYQQFYHNRvijDzoHDNy/jz6//AAOB+Db22AlSQghREVgMJB1wy3Yw8JxBASw4aHHwTWER2JKqs/CkltghRCiIjAaybr1DnIGD8W8YR2OBs1hw1YAVm3cRnpGJn1bNyegenXvhuXVswkhhCiWvV59cq4aR7OGsQQF/DWCbPrsWVTv1Iakr770ajySJIQQogKqUb0at4wbQePYegRkZzH+528JS02m2VOPcPCJx8nIyvZKHHK5SQghKqjw0BDGDu7DUXsWZrtzGPL0kFDm12+G7effGN6nG41j65VrDFKTEEKICsxgMBA7bChJKzegB49kzuiJpIVXIzMrmzlLVpGVnVOu55eahBBCVALBMdEET/2BhkeOsXfVetLSMzCZjARYyvdrXJKEEEJUIo1j63H7+BiOnjxNdI1IjEWMIluWJEkIIUQlE2CxlHtbxAXSJiGEEMItSRJCCCHckiQhhBDCLUkSQggh3JIkIYQQwi1JEkIIIdySJCGEEMItg8M1qUUldgY45OsghBCikmkARJe0UVVIEkIIIcqJXG4SQgjhliQJIYQQbkmSEEII4ZYkCSGEEG5JkhBCCOGWJAkhhBBuyXwSFZRSqhvwhta6v1KqKfA14AC2Aw8CccAPWuvuru17A18BVwOPudbN80Xs5aGk8tBa25VSLwKjACvwiNZ6rVLqa1xloZQyA1OBs8BbwPcXyq+qKVReXwNjgVpa62zX+k7ABmAA0BBoobV+2kfhljsPP08JwNNa6zfy7TcLqOba7yDOcsrycvg+JTWJCkgp9STwORDkWvQu8JzWug9gAMYU2r4/8BkwSmu91YuheoUn5eH60usHdAOuAz4qdAwLMA04oLV+AOcXRJVURHkBnABG5Ht+I3DAm3H5Sik+T/uBa/LtVwNo5sVQKyRJEhXTfmB8vuedgaWux78Dgy+sUEoNBj4Ghmut9+Tb5wGl1GKl1FLXL6fKzJPy6A0s0Fo7tNaHAbNS6kJv0kDgZ2BzoV/L0UqpWUqpeKXU8+X7EryqcHkBfA9cD6CUMgKdgHX51vdwvV/WKaVGeSdMr/H083QWOK2Uaul6fi3wU6Fjfer6TE1TSgWXV8AViSSJCkhrPQPIzbfIoLW+8Ms3FajuetwEeA3nL6TCb9hVWutBwBvAm+UYbrnzsDyqAcn5tslfTh8AoUD9QocOA24GegEjlFLtyzh0nyiivADWAkopFQoMBJYUWp+O88tyFDDZlUiqhFJ8nsCZTK9zPR4DzCx0uP9orfsBB4G7yz7aiqfKvBGqOHu+x+FAkutxJs5LCA8DhX/ZLHP9uwpQ5R6hdxVVHimux4WXA3wIDAXaKqVuyrfNFq11stbahvNLtHn5hVwhzML5xXcD8F2hdStctbDTOJNtTW8H50XuPk/gTApXKaUaAieBjHzrcrTW8a7HVfFzVSRJEpXDJle7AziTwnLX4+Na6/Na6zmuZZPz7XOF698+OBvnqpKiymMlMEwpZVRKxQFGrfVZ1zbbtdZWnNfh38x3OaGlUirM1aDdDdjhvZfgE98BtwB1tNb7C63rCqCUqo2zhnWWqsvd5wmtdRqgcda+pxbaL0Ap1cH1uCp+rookSaJy+D/gn0qp1UAAMJegcFsAAANxSURBVL2IbR4HrlBK3ep63l0p9QfwCPCkd8L0movKQ2u9AeeHfTUwA+cdKwVorQ/gLIufgBDgPPAjzl+F07XWO70T/v+3d/+gdZVxGMe/hUToIkGCbUHEpTyIIDhYi0XI0tSlFF0qghAcmhAtERepU21qpm5CFUuHlBYSt7ZIzGC4OBSHoriE/BZRtNj4r9IlNqB2+L2XHE9ykt6kYuh9PsvlnPec854LSX685+R93v9HRASZ+nl1jead5eflCjBceRzzINro9+kSWQQ+r+2/AxyX9AXwKPky/IHnFFgzM2vkkYSZmTVykTAzs0YuEmZm1shFwszMGrlImJlZIxcJ6wqSBiRNrbF/StJD65x3s8N+dks6u5l7rFzjBUljTffccM5OSZOSdmylb7M6p8BaV4uIVzY+qqPr3QRGN3t++SN/kpzk9XwH/S5JukZOlpvcbP9mdS4S1tXa8c/AR+RkqSeAPcBQRHxVOW6CzPh5sz3RrAQITpMj8l5ghMwCmiIjuGfK6T3AfjLG4TEyb+svMnhuOCKquUIHgfmIWJYEsFfSLBmT8WFEnJfUAhbKfe8Ajpbi9AnwGS4Sdh/5cZPZiu8j4hCZ9XSsvVPSGaAnIt6ozUTeR+YctfOzHm43RMRSRAyQxeI7cnTxLXAOeLmExN0Ahmr3MABU4957gcPkDOB3Ksm218r1p4F3S5+3gH5J1cA6sy1xkTBb8XX5/IGVtQd2AU+TeUZ1M2Tk9GXgFP8Ojmv7gEzE+JiMxNhDhjG2yNDBx2vH9wOLle0vI2I5IpaAeXKkAzBXPutBc4vAI43f0KxDLhJmK9bKqFkEDgFPSXqx1jYA/BQRg8BpYKLaKGmcjKUeL7t+BX4EjpRRwPusjuz+GeirbD8jqadEfD9JPqKCXBMBMua8GkzYB/yyznc064jfSVg3GZR0vbL96r2cFBH/SHodmJX0XET8Vpq+AaYlvUW+YzjVPkfSs8AJoFVGDZT2MeDTsl7DbfJFc1ULeAm4ULb/JEcsfcDJiPi9vKsYkvQ2uQ7Ea6XPPuCPkmRqdl844M9sGynFYw4YjIjlhmNawEhELNT2jwK3I+Lif36j1jX8uMlsG4mIv4H36PDfaMuCUwdYvQaC2ZZ4JGFmZo08kjAzs0YuEmZm1shFwszMGrlImJlZIxcJMzNrdBf367fwrBwJlAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "A, B = 0.000347, -0.8042\n",
    "title = \"Rice (RE=GATC)\"\n",
    "label = r\"$Y=0.000347 \\times X^{-0.8042}$\"\n",
    "filename = \"Rice_RE_GATC.distribution.txt\"\n",
    "plot_distribution(A, B, title, label, filename)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAESCAYAAAABl4lHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3XV4VFf6wPHvWBRCEpLgLgd3irsVpxRqQLdu1Lfb37bbLXVl6wbbrXuhFC3uBPdAOUAgaIAAcR/5/TFDmhAbIMkkmffzPHnI9ffchHlz7zGDw+FACCGEcIfR0wEIIYSoOCRpCCGEcJskDSGEEG6TpCGEEMJtkjSEEEK4TZKGEEIIt5k9HYCoPJRSFuA4sEtrPdy1riEQDezNtasBeF9r/YVrn/uBRwEbcBS4W2t9XilVDTgLHMh17BNa61UFXLsO8CrQCXAAGcBrWuu5ru0OIFxrfd7NsvgCC4AZWutZl217BPgA6KG13lTMeb4F3gRSCrgPVYCTwF1a6yNKqTuA9133ILc9Wuvb3Yx7NDAPuEVr/fOVbFdKVQdeAIYC2YA/sBj4p9Y62bWPCXgMuA3n54cPMB943rX/6lxlqwNo1/IyoBHQ1LXcHojC+TOPBwYDc1334pw7ZRWeIUlDlKTxwC6gi1Kqpdb6T9f6dK11h0s7uT7go5RS24BknB/2Smt9QSn1PvAiMBXoDqzVWg8t6qJKqXAgEngOuFNr7VBKtQeWKaXStNbLrqQQSqkewMdAC2BGAbs8AHwPPA7cUsR5bgIStdZRruR5+X0w4Ew+rwK3ulav01qPupJ4L/OQK7YngHxJo7DtSqmqOO/h90AbrXW2UsoH+I9r3RjXrp8CIcAgrXWiUirQtf1zrfUUoIPrfP2Bj3KXNzdXEh+QO4krpd4CPgEmXF3RRVmQ11OiJD2I86/Fn3H+NVogrfUp4BDQHDABFqCqUsoIBOB8SgDoCYQqpTYppXYqpR4s5JQPAeu11t9qrR2ua+zG+eFzJtd+LyqltiulYpRSUwGUUncopdYppXYopS49wTwK/BPYcvmFXB+GocDTwFilVL0i7seLwGdFbPcDagMXi9gn97U/UErtuuxrc67tjYH+wJNAM6VU98uOL2r7fcAhrfVLWutsAK11FvAPYItSyuhKfJNwPgkmuvZJxZlEf3enDEXRWq8FWimlCkw0onyQJw1RIpRSrYAewI3AdmCNUupfhezbA+dris1a6xNKqbdxvsZIABJd5wGw4nz18QYQBqxSSsVqrS//gOoC/HH5dVwfQrkd0VpPVUp1BDYqpWa61rcGGmqtk1zH3eqK89kCwn8I+F5rfVoptRJ4GPi/AsrYBvDXWkflWu2vlNqFM1FG4EwWvwGv59qnj2uf3N7XWn+ptX60gHhyewBYqLU+p5T6CefTxM1ubu8DLLn8hFrrDOAVV5k6A/su3adc+5wBZhcTm7uWAzfgfGIV5ZAkDVFSHgQWaK0vABeUUkdx/vX6I399WILzd+48MMmVMIbiTDT1XOvfBL4CRmutX851/lNKqRk4P1AuTxp23Htq/sH17y7AFwhyLe+5/IOwIEqpmsA4oLNr1dfAZ0qpl1x/cefWAjh82bqc11NKqWHAd8B8rXVKrn0KfT2llPoA6HvZ6kytdTdXHcydwF25YtuglKrnus9FbsdZz5QzppBSahLOpwyAcGAE7t/na3EU6FbK1xDXQJKGuGau99pTgEylVIxrdRDOv8Jncdm7/MuMAeZdqvxUSn2Ms4L0UoXzXK31cde+BpwVtJfbhLP+46PL4rofCNRav+Nadem1i0Mpdel84Kykdse9OD9Y57uON7rK+Tec7+Jzc1DEB6zWeolS6h3gV6VUa3eSVjFPGjfhrGv4SCn1Ya4YHsH5Kq247ZE4X1195LrW9zjrKnD9TC3AZqClUqrqpYpx1/Y6wExggtY6vbhyFCMbZ+W4KKekTkOUhEnABaC21rqh1roh0BhnC5qJxRy7AxiplKriWr4RZxIA6I3rr12lVChwNwVX7s4A+iulJrkqly+9SnmJvK2Vrpqr1dC9wAOXyqi1rg+8Bjx26bq5aKBJMaedjrMhwIslEOKDwKta6wa5fgYPAPe6knpx2z/FWZ/wD9dTCa56jKE463BsWuvTOBPJF0qpINc+QTgT5oUSSBjgbGF1oNi9hMdI0hAl4UHgHa11zl+IWusEnC2Dnijm2C+BRcB2pdQenH/t3uHa9jBQRym1D2ci+bSgllBa64uu427E2SprL87WT3dfacupIozC+f/l+8vWvwvUxPn6JndMUUC6UqplYSd0VTg/DDzsqgMBV53GZV/bigrM1VKsA/DhZZu+wdmc9Z5itt/hetLpCdTAWfG9G4jB2YptotZ6p+uYh4D9QKTrleNm1/I9RcV4BYbifDoV5ZRBhkYXonQopW4DemutH/J0LBWBq2XaVK11cU+nwoPkSUOIUqK1/gFnk+G2no6lvHO9/nsaZ3NnUY7Jk4YQQgi3yZOGEEIIt0nSEEII4bZK10/Dbrc7bDbveOVmMhnwlrIWxtvvgZRfyl9S5bdYTOdxduQsUqVLGjabg4SENE+HUSaCgwO8pqyF8fZ7IOWX8pdU+cPDqx5zZz95PSWEEMJtkjSEEEK4TZKGEEIIt1W6Og0hKjObzUp8fBxWaxYAZ88a8Oa+VlL+Ky+/2exDSEg4JtPVffxL0hCiAomPj8PPL4DAwJoYDAZMJiM2m93TYXmMlP/Kyu9wOEhNTSI+Po6wsFpXdU15PSVEBWK1ZhEYGITBcPmgukIUz2AwEBgYlPOkejUkabhkW60kJrs7rYIQniMJQ1yLa/39kddTQGZWNl/OXkBqegY9OrahZ0cZX04IIQoiSQNISUsnNT0DgI07ozAZjbRTTfH38/VwZEIIUb5I0gBCq1Wlcb3aHDlxGoD12/cQuTOK+rVq0KxhXdo2byKvBIQQAjC98MILno6hRNntjhcyMgqaRrpwBoOBRnVrkxQZSdvIVRyv3xi7wUBCcgpHTpwmPTOTxvVql1LEV8/Pz8KVlrWy8bZ7kJKSSJUqwTnLRqN3NzktqfLv2LGNl19+nv37o/D396dWrbL5/56ens6bb77Cxo0bSEtLpUmTZjnbzpw5w8sv/5utWzcTExNDWFgYH3/8HkuWLGTgwCGuvRzMmPEx69atJiEhnmbNmpOens79999JWFg49es3KPC6l/8eAQQG+sbinOu9SFIR7uJnNjFpwU9cv2QuT3zyBjVjT+Zs2/XnIQ4ccWtYFiEqvZiYo0ycOAa73dnU026388QTU/njjwVXfK7XXnuRUaOGMGXKTYXuk5yczHPPPc1tt93IpEkTiIraA8DZs2eYOvU+Jk2awOTJN/HLLz9eXYFw/uHo7x9AVlYWERE1rvo8c+bMYvr0N3KWZ878hJdf/neh+69Zs5L+/Qfxf//3HOvXr82z7cSJY/To0Ztnn51GTMwR6tSpyzPPPJ9nn7VrV3P+fBxms5nw8AgAvv/+awYOHHzVZSiOvJ5y8f31J3x27wIg7PxZxo0cwryTcZw+dx6A1Zt30LhebXwsFk+GKYTHNWzYiPr1GxIZuY7evfsxY8bH1K/fgOHDR13xuUaMGM2NN97MK688X+g+778/nW7devLKK2+RnZ1NRoaz/tFkMvPoo0/QtKkiLS2Vu+6aQteu3WjUqHGh54qOPsyMGR/lWffMM8/Tvn1HOnbszMWLF/jww3eZNu2VnO3x8Rfx9fUlICAwZ93JkyeoW7devvMPHz6K2267keTkqezZs4uNG9fz2WdfFHrduLhzNGnSFACjMe/f8M2bK7755gtWrlzGsGF5pqDPcfz4MVq3bsu4cTfy3HNP43A4aNiwEVlZV9+ktjiSNFwyJ9xMSnw8gW+8TNpjfyewQwfGt8rmq98W5lSUf/jtLFo2aUDvzu0JqhJY/EmFqKRuvvk2fv75B6xWK3v37uaDDz67qvN06NCJ2NjThW5PTU1h9+6d/OtfLwBgsViwuP5wCwsLo0aNCGw2OwEBgTRs2JDz58/RqFFjTp8+xQcf/Ie4uDiMRgP//vdL1K/fkCZNmvLWW+8Ver2qVYPIzs77gbtz5w5+/30Wb7/9Pr6+vsybN4e1a1cxffoH+Y738/Nj8OChzJz5CZs2RfLeex/j6+tX6HXDwyM4d+4czZopHI68nfQWLpzP3XffT4cOnXjuuacZOXJMvuMjIiIwGp0f40ajiR07tpGRkUFMzBF8fX3p0aNXvmR0rSRpXGI2k/7gw2SOHI29prOnpK+PhV6d2rFk/WbqnDxGQnAof0YfI/r4Ka7v051mDfP/pSGEN7juuu589NG7fPbZx3z00UzM5rwfJQ89dA9pafmH7J469TG6du3m9nVOnz5FcHAwr732IocPH0Spljz22FP4+/vn2S829jQHD2patWqD1WrlzTdf4emn/0WdOnXZuHE93333Nc8+O63Q66xZs5LNmzeSkpLC+PF5X5UNHDiY2NhTvPDCswwYMJiFC+fx7rsfF3qukSPHMmnSBN544z/UqVO3yPL16zeQd999i40b19OrV18AoqL2cPCgpnv3HnzxxUyWLVtMzZq1SUxMYObMTzh4UPPtt18yZcqd9O8/kOnT32TPnp20b9+RG290xr5o0XyqVQsu8YQBlXCO8Oxsm6Mkx9e32+3MnbuIif/3MD5ZmSwYdRM7O3YDg4HOrRW9u7THbDKV2PWuhLfPJQDedw/OnDlGzZp/VW5Wmf46/m+9DkDqU/8k7eln8+wf+PyzBHzmfC2S8sKrpD/0SJ7tVf7+KP7ffgVA8vT3ybj9TrdjeeutV2neXDFu3ISrKUqO2NjTPP3043z77S/5th04sJ/777+TTz75H61bt+G996YTGBjIvfc+CDiH0UhOTuHhh+/jb3+7i379BrJy5XLee+9tQkOrA87xutq378hTTz1zTXFOm/YMGzdG8vPPvxMSElLofi+88C+2bdvCP/7xDP36DbymaxbnaodRufz3CCA8vOp2oEtxx8qTRjGMRiOTtq0lMOEiAKMX/sqBFm1JDwhk+z7NgSPH6dOlPa2bNfJwpEKUrZiYo4wYkf+VCZTck0Z4eATh4RG0bt0GgAEDBvHdd1/lbLdas3nuuacZOvT6nA/ow4cPct99DzJq1LgrKE3Rdu/eyZEj0fTt258vv5zJk0/+X4H7/fjjd2RlZfHSS6/zxRczSz1peIIkDTdkDxqCbfEiTCeOk/7aW9QOr0/08VMApKans3jdJuKTkunZsU2pPA4KUR4dPXqExo2bFLjtk08+L5FrVK8eRkREDY4fj6F+/YZs27aFhg2dFd0Oh4NXX32JBg0accstk/Mcs3nzJkaMGIPRaCQ6+jCNG199X6uDBw/w5puv8Pbb71OrVm1eeunfzJz5Cffd91Ce/bZv38qiRfOYMeNLAgIC+fDDdzh0SNOsmbr6G1AOyespd6Wk4DfrZzL+dhcOYOefh9iyex+p6Rn4ZqST6edPjbBQhvTqSkRoSJl0BvS2VzMF8bZ7cPlrBU+N8nr27Bkefvh+fv117jWdZ9q0Z9m1azsJCQmEhlbn7rvvY9SocTz11KP885//JiwsnEOHNG+88QpWaza1a9fhmWemERQUxO7du5g69R6aNGmKweD8Y+3++x+iU6cuvPzy80RHH8bX14/GjZvw/PMvX3WMe/bsIjCwSk4rJ6vVyqJF8xkz5oacfc6cOcMjj9zP22+/R8OGzrcOixbNZ+fO7TmV+KXBE6+nJGlcg6zsbFbOmsv4Zx5ha9derBg0mmwfH6oE+NPvuo60aFxwx5qS4m0fmAXxtntQXpJGeSHlL/ukIe9SroGP2czEP2YRmJZK/zVLmfT9DMA5ltWiNRvZf/ioV/fWFUJUPpI0rkVaGsbk5JzFAzfemvO9w+Hgj7Wb+OWPFVyIT/REdEIIUeIkaVyLwEASf1tA8rsfkfbQo/R6/BEeuGVcno5/J8/E8e3cxew/fNSDgQohRMmQpHGtDAYyJt1O6gvOYQcCA/y5bfQQ2rdoRu3Yk0z55hOqXDzPH2s38esfKzl3Id7DAYuKTl55imtxrb8/kjRKQaC/P4Ov68C9K+bSev9unnznBVpH7eR47Fl+XrSCU2fjPB2iqKDMZh9SU5MkcYircmmOcLPZ56rPIf00Soklcj3+rtE4zQ4752o6h1rOys7m50Ur6NKmBX26tJd5OsQVCQkJJz4+jpSUBMA5Oqs3JxAp/5WX32z2ISQk/KqvKUmjlGT37U/8/KVU/fsjZN54E8Om3M7spatJz8jE4XCwde+fJKemEVE9BF8fC80a1sPfV2YKFEUzmcyEhdXKWfa2JseXk/KXffmln0Zpy8wEoxEsFuKTklm6fgsnz5yjcbTGbjQS08g56UqVwAAmDO1P9ZBqbp/a2//DgNwDKb+Uv6TKL/00ygtfX3AN5RwSVJWJ1w+gda1wJv76FQ/MmM64Od/jm5FOSmoavy5eSUZm6Y2DL4QQ10qSRhkzGo2M37SSENcAiB3+3I2fq1ojNT2D9dt3ezA6IYQomiQND8i4fyqZw4YDkPX62/QbPiRn2+4Dh9mjD2Oz2TwVnhBCFEoqwj3AXrsOSd/8hGXDOrJ79aEp0KhuLY6ejAVg15x5rK5RmyF9u9OySUOPxiqEELmV+6ShlOoK3AkEAf/RWu/0cEglw2Agu7dzpi4DMKxPN35etALb8eNM/fgNTteux+9xt5M9dhQ1qocQGhyExVzuf1xCiEquInwKdQZaub4/4clASlOgvz83XT8A48Tx+GVm0PjoIW79bgbvVw8HgwEfi5khva4r9ZFzhRCiKOUuaSilHgcGuxY3AsuAL4CBwEjgaw+FVuqq+PkR2Ksnjl1bMVitzBtzM7g6/2VlW1m4OpLlkVsxYKB9i6aMHNTdwxELIbxNuUsaWuv3gPcuLSulfgImA3H89cRROZlMpP5rGhljx2NfuphqPQYQfjaOuIsJObtkp2dgN5nYvGc/RrOBnh3aeTBgIYS3KdOkoZTqBrypte6vlDICnwDtgUzgHq314QIO+wn4FsgC/lFmwXqQrU1baNOWQa7ljMwsVm/ZQdK69dzyzSfMG3MLukVbNu7YR2hQNXllJYQoM2XWI1wp9TQwBUjVWndXSo0Hxmit71BKdQee0VqPvdbr2O12h81WuXq5A2C1YuzVE9POHQAsHjaO1QOczXZr1whjwvB+hIcGezJCj5CZ26T8Uv6SKb/FYnKrR3hZPmlEA+NxPjUA9AYWA2itNymlig3WHTabo1IOK2A6qAk+fgwAh68vMdf1ytl2+ux5fpy3gpH9ehJSrSpGo/d0v5FhJKT8Uv4SG0bErf3K7NNFaz0byM61KgjIPaWdTSlV7upYygtbc8XF9dvImHgLqU/9k563TKBKgH/O9riLCXw1ZxE/LVxOemamByMVQlRmnvyTNAnIndqMWmurp4KpCBzVq5P88UzSH32SWuHVefr+W+ncWgHQXEfRd80Szp45x+zFq7FKj3IhRCnwZNLYAIwAcNVp7PVgLBWLqxmuwWCgV+d2NKgayPjfvmPEH78x9ePXsR34k6iDRzwcpBCiMvLk66A5wBClVCTOTtF3ejCWCstiNjPl1EGqJDqnka2WGE9aQCArNm7DZrPRvkUzzGaTh6MUQlQWZZo0tNYxQHfX93bggbK8fmWVPvUxMJoInP46S2+4jbTAKgCs3rKT47HnGDuot1dVjgshSo98klQGFgvpjz7BxU07Cbzrnjybjpw4xc6PP8V69pyHghNCVCaSNCoRe+06dGytGDOoDzXDqgMQlBhP/7depFrPTph//Qkq2UyNQoiyJUmjkjEYDDRrUJfbRg+hnWrK2Lk/4peZQUByEvaXprE36k8SkpI9HaYQooKSpFFJGQwGhvTqStptU0ioFgLAr+MmsXTrbr6YvZB123ZzMOYEqenpHo5UCFGRSGe6Sq7Bvfcyt3YDLMuWENOoGQAOh4Mte/YDEGbNZPykm6gaGODJMIUQFYQ8aVRyBoOBodcPpubUqXTv0CZPcqh96jiPvfB3Yu+/j8MyxawQwg2SNLyA2WyibfMm9OrUlinjric8NBijzcb4377FYs2m19J5WB5/mM9/nU9s3IV8x+/Yf5Cv5yySDoNCCEka3sbf15fJY4Zxz+A+BIQ46zqyzWZWDbielLR0Fq7akDMESWp6OifPnGPVpu2cj09k+catZGRmeTJ8IYSHFVunoZQKAXoB1YFzwHqttTS/qcCMRiNVGzXAtnwNpz/9mPOxZ7kQVgOAxJRUflm0AqPRyKmzcXmOs9nsHIw5TjvV1BNhCyHKgUKfNJRS4Uqp/wELcA5p3hwYDSxVSs1UStUooxhFaTEasUx9hFqvvMKAbp1yVsfGXaDq6hVM+eYTglzDk1yy/3BMGQcphChPinrSmIZzlr2Dl29QSrUEngemllZgomy1b9GUfYePcu5CPD6ZGYz7/QeCE+NpEq35dsoDRDdtCcCps3FcSEikenA1D0cshPCEQpOG1vrhIrb9iSSMSsVkMnHLyMGcibtA1QVzCb70hOHnh6FjR6rajSSnOid72br3T+rWjCD23Hm6tm1JcJB7k7cIISo+d+o0egEfAzWBk8C9WuudpR2YKHsWs5l6tWrAvfeR0LoVVf7+KFn/eIax42/g1Nk4flq4HIB9h46y79BRAOLiE7lt1BBPhi2EKEPutJ76ALhNa10TuANnAhGVXHbP3sSv3kjmDRMAqFMjnNoRYQC03bON+seiAYg9d57z8YmFnkcIUbm4kzQStNb7AbTWUYD3TsjrbXx9cyZ8AhjQrRO1MlK5cfa3PPDZ24yZ+yOWrEyWrt9MtlUmXRTCG7gzjMg5pdTnwEqgM2BUSt0HoLWeWZrBifKlZnh17tuxDr/MDACaHj6A3WgkNu4CX85eyNhBfQgNDsJsMmHIlWyEEJWHO0njgOvfpjjn9V4D1AJkjG0vlDbtZYxpqfguX8qSyfdhM1sASE5N47t5SwDwsZgZ0L0zbZo19mSoQohSUGzS0Fq/qJQajbOfRpTWeknphyXKK3u9+iR9/yvmXTvo0VThu0+zR0fn9CKveyKGk3UbsHrzDpo3rIePxZlUbDYbmVnZ+Pv5ylOIEBVYsXUaSqmPgUmAFbhHKfVWqUclyjeDAWvHzlSrWoUB3Ttzy6jB1K0ZQcNzsTz46Zvc+993qXL6FB99N5u9B50JZdaS1Xz64xw2bN/j6eiFENfA4ChmJjel1Hqtde9cy5Fa656lHtlVys62ORISvKOuPjg4gHJTVquV4OsHYtmzC4ADqg1f3flIgbs+NGk8/r6+JXLZcnUPPEDKL+UvqfKHh1fdDnQpbj93Wk8dV0rVBXANHXLiGmMTlZHVSnaffjiMRqxmM/NH31TorgePyq+QEBVVoXUaSqlYnJXdfsANSqljQD3gfBnFJioSPz9Sp71M5g03krRpM6FNW3MhJm9yMFmzsZktRB2MplXThljMMgeYEBVNUcOI1CrLQETlYG3XgYB2HRiDc4bAvTqaNVt30WT3NkbP/5k5N0ziEK35es4fDO7Rhdo1wriYkERE9RCMRhmpX4jyrtg6jYpG6jTKH2v8RcL798QcexqA2eOnsPW63nn2adqgLmMG9r7illUV5R6UFim/lL881mkIcU38og9jzEgHICs4hEMdu+bb5/Cxkxw5cZr4pGQuJiaVdYhCCDe50+T2qbIIRFRe1i7XcXH9NjLGTyDjjencPOVmfH0s+fb7fflavpi1gK9+W0T08VMeiFQIURx3njRGKKVMpR6JqNQc4eEkf/YFmeMnUiXAnx4d2+Rs67plHb3XLcdgtzv3dThYv303le3VqRCVgTvNV8KA00qpozhbUznKcz8NUTF0aNmc+MRkMg5HM2bRbCwZ6bTfvYVvpzxIUrUQzscn8s6XPzG4Zxfat2jm6XCFEC7uJI3RpR6F8Domo5HBPbsSuPBXLK76Dt/MTFIDq+TZb+XG7TSuV4eqgQGeCFMIcRl3Xk9ZgVeBmcBQnJMxCVEiUl94ldRnn8fh50fMv1+hT8/r8vTfsDscbI86UMQZhBBlyZ0njZnAf4B/A2uBr4HupRmU8CIWC2mPP0X65DtoFBZGI6B+7Zps2r2Pg0eP02bvdvanpWIwGOjTpb305RDCw9z5H+intV6Jsy5DAxmlHJPwQo6wsJzvw0ODGdW/Jy0TL3DbD//lyXdeIPOH71i4agNHT54mJS3dg5EK4d3cedLIVEoNA0xKqe5I0hBlwGC3c/PcHzA6HFRNSaL7prX8t21nDh47CUDHVs0Z0K0Tu/Yf5tz5BFo3a0Sgv7+Hoxai8nMnadwHTMfZiuop4MFSjUgIAJOJzKefxfLPv+O4cIHfxk/OM/Xszv0H2bn/YM7y5t37GTWgJ43q1vZEtEJ4jWJfT2mtTwJvAi8DL2qtj5Z6VEIAWSNGEb9+C+c+/xpjs+b4+fjkTOoEUC3h4l/7ZmezeO0msrNlrnIhSlOxTxpKqVeAAcAW4FGl1Byt9dulHpkQgCOoGpbhI7ndtWyz25m9ZBVVVixj0vczWN1/OKsGXI/NbCEtI5NdBw7RtW1Lj8YsRGXmTkX49UAfrfUTQB9gYumGJEThTEYjN3bvxM2LZmG22Ri8YgF3RG3J2b5+2272HTpKUkqqB6MUovJyJ2mcBKq6vrcAZ0svHCGKZ0lJxtSgAeAcniTo39MIqhIIOPt1LF63iS9mLZDxq4QoBe5M97oFaADsBloBWcAZgLIYTkQpdQswCrgAPKe1Ti5qfxka3UvYbPh99Tn+DesTP2g4Z85f4OdFK7BabeBwgMGAxWzm5pGDqFE91NPRlhqv/h1Ayu+JodHdaT3l6ddRY4Dbgfaufz/2bDiiXDCZyLj7fvyCAyAhjZph1Rk3uC+rNm2n+aK5NDp6kPmjb2HW4lVMvH4gEdVDPB2xEJVCsUlDa32sLAK5RCn1ODDYtbgR+BD4L865yW1lGYuoWBrUrsld3doT/OQ9mNJSaXboTz6/5wl+X27m9nHD8fP18XSIQlR45W6SZq31e8B7l5aVUuOBe3C24KrnqbhExeCzeBGmNGcleHK1YM7UrI0tNY1fF69kSK/kVkC5AAAgAElEQVSu1Ayr7uEIhajYrmggH6XUNQ38o5TqppRafelcSqnPlFIblVKrlVJNCzksAfgKmAz8cC3XF5Vfxt33kTB7PtbGTTj7ylvYzM5+HecuxPPj/GWciD3n4QiFqNjcqQifCJgAX+Bt4C2t9fQrvZBS6mlgCpCqte7ueoIYo7W+wzU8yTNa67FXXILL2O12h83mHZP3mExGbDa7p8PwqELvgdUKZjPLN2xjzebdOasHbFtP3bGjiG3UlE6tmlG1SsUect3bfwek/CVXfovFVGIV4U8BI4CfcL4eWopzWJErFQ2MB751LfcGFgNorTcppYoN1h02m8NrWlN4e8sRKO4eZNGldSvqRtRg9tLVVI8+xJDZ38Hs74jvOYDPJkxh8o2j8PfzLdOYS5K3/w5I+Uu09ZRb+7nzuunSAIXJWutM/uqzcUW01rOB7FyrgoDEXMs2pVS5q2MRFV/N8OqM6t+L0Qt+wehwYHQ4qHnmFEkZmfyxdiM2m7SvEMJd7iSNo8A24Aul1DRgcwldO4m8CciotZaBg0SpaFCnJmmffk50i7Zkmy3MucE5AOLRk7HMW7ne2b9DCFEsdwYsvANoq7VeAMwAni2ha2/A+doLV53G3hI6rxAFCu/cicBVa7mwdDVNB/bPWX/kxGl+W7qarKVLnB0DhRCFKjZpKKU+1FqnuBbbU3JPGnOADKVUJPAu8EQJnVeIQplMJkxt2tKrczu6tWuVs77qiqXUmTyRjIF9ObdtmwcjFKJ8c6cOIUkp9QZQBWiDcwDDq6K1jsE1VazW2g48cLXnEuJaGAwGendpj8ViZuuGzYyd+yMA9fbtZvML09jy92fp0LIZdWtGeDhSIcoXd15P/Qtnk9umWuv+WusjpR+WEGWjW/vWjBvWn6M9+mI3GEiuUpXF19+APnqcnxet4GDMCU+HKES5UuiThlIqFrj0gtcA1FBKnQbQWsv0aKLSqNukMXz7HSfWr2P3+k2kBwTmbFseuZU6YSEEmkwg08kKUXjS0FrXKstAhPC0gN596N6rNzVOxjJ/lbNFVXpGJtHPPEvvjatIf+dDsvv293SYQniUOzP3dcA5T7jfpXVa67tKMyghPMVgMNC4Xm3GDurD7CWrqZZwkd5zfsQ3KxPfCWNI+PJ7skeO9nSYQniMO/00vgJ2AD/n+hKiUmtYpxbDencj4twZ7EYTAGcjarG0anWyrdKdSHgvd1pPndFaf17qkQhRzrRp3piqj07l29at6fb1DDb0GsTxw8c4cOosg3t0oVlDGXRZeB93kkaMUuqfwE5cFeNa66WlGpUQ5USD2jWpN+kmZodHcPy0c6bjtPQM5q1cT8dWzRm5+g8cYWFk3HkvmEwejlaI0udO0vAFlOsLnIlDkobwGkajkRuG9GPX/oNsi9KkpqcDcG7pUgI+no7RYcdv9q8k/vArjpDKO7WsEFB0k1uzayyo+8swHiHKJbPJRJe2LWnTvAlLN2zhUMwJ+q5dhtHhHJbaHhCAI1imlBWVX1EV4d+4/tXAAdfXpe+F8Ep+vj6MHtCLDi2b8ctNd7Bs8CjS/fz5ccRETp6N83R4QpS6Yidhqmiys20Obxlf39vnEgDP3QO73c6C1ZEcijmBb0Y6mX7Ojn8N69RiRP8e+Pv64vfNl2SOG48jqFqpxeHtvwNS/hKdT+PaJmFSSm3krx7heWite159aEJUfEajkdEDerF59z427ozKGR035lQsc5auoffJI3Se9jQBb79Oyhv/IUv6dohKoqiK8FvKLAohKiCDwUD3Dm1o0bghm3ZHse/QUQDOnzpN0/+8BoDp7BnMSxZJ0hCVRlHDiBwry0CEqKiCg6pwfZ/upGdkcuTEabItPswffRNj5v2EwQG/9LmeETYbZmmSKyoBmV5ViBLSrX1rjpw4DQYDUW07E92kBWHnz3EiOY0/1mxkZP+eGI1GsNkwHj+GvVFjT4csxBUrtPWUUupV179jyy4cISqu2hFhDOzemYZ1anHT8IF06t2DE/UbAXAw5gRfzVlEbNwF/L+YSWif6wh45y3IyvJw1EJcmaKeNMa6hkJ/RClVI/cGrfXM0g1LiIqpY6vmdGzVHIC6NSNIy8hg5/6DAMQnJrPy59m0mf4ShqwsAt94BXtQEBn3yFxkouIoqp/GnUAwzh7htS77EkIUw2AwMKBbJ3p2apuzzpqaSlwN53+hjKbN2dNnCDa73VMhCnHFiu2noZTqCkQDTYCjWuvzZRHY1ZJ+Gt6lotyDP6NjWLRmIwAGu52ekSs5Ua8Rxxs0oX6tGowf2g+TyQR2OxjdGXzaqaKUv7RI+cu+n4Y7v52NgE3Av4BNSqnJ1xibEF6nReMGNKzjfMJwGI1s6D2Y4w2aAHA89iwrNm4HIPCF56h6/50Y4qR3uSif3EkaTwCdtNbjgI7AY6UbkhCVj8FgYNyQvgzp1ZXwkOB826MOHSFu6RL8Z36C35zZhPbqjPHoEQ9EKkTR3Glya9dapwBorZOVUhmlHJMQlZLJaKSdakrrZo05ePQ4druDqENHOHnmHA6Hg1NffUMrV/2GtV1H7A0beThiIfJzJ2lEK6X+A6wF+uKs3xBCXCWT0UjLJg0B8PW1cPLMOQBWDB7F8fqNGbvkN3j7XTAYPBilEAVz5/XUXcARYIjr33tLNSIhvEjjurWpViUwZ/lQ81ZMf/hfRBZQuRnw2kuY9+wqy/CEyEdGua3AvL3lCFSOe3D2/EW27P2Tk2fOkZb+19vfIT270q5FUwB8Fi2g2h234TAaSX/gYVKnvQwGQ6Uo/7WQ8pejUW6FEGWjRlgoowf0wmaz8dvSNRyPdU4ru3zjNrKsVjo0b0zIi88Bzua6hsQEeXUlPKbY11NKqYDc/wohSofJZGLs4D7UqO6cMtbhcLBmy06+W7CMM1/+QFafftjDI5xPGUJ4SFFjT72mlBoHLHStWlA2IQnhvXwsFsYP7UdotaCcdRcSkvhJx3D2+1nE/7Ei/7SyNhs+S//ImdNDiNJU1JPGbKAm0FYptQxooZR6WCnVvWxCE8I7Bfj7MWnMUHp2aovB9RrqfHwic5avJbF6eL79/T//jGqTb6barTdiPHG8rMMVXqaopBEIfAfs0VoPAY4Cx4FhZRGYEN7Mx2KhR4c2DOvdLWfd6XPnmfnzXLbu/fOvHY8dI/D1V5zHrFyO368/lXWowssUlTR6Al8D7ZRSn+IcvPBPrfWLZRKZEILWzRpxXbuWedat3borZ5ZAQkNJv20yDoMBa4uWpD38uAeiFN7EnQELV+Hsq7EQWAE001pfXwaxXRVpcutdvOEeOBwO9uhoNu2KIiUtHXDOUT5hWH/atmxEQkIa5q2bwWLB2qFT3oPT0pwDIPr5eSDy0ucNP/+ilKsBC3O1lvq71vooMFlr/cilhCGtqYQoGwaDgfYtmnLH+JGEhVQDwG63M2/FeuIuJgBg7dotf8IAAl9/mZD+PbBEri/TmEXlVdTrqY+VUlOBYwBa6x0ASqkwpdTjwKdlEJ8QwsXXx8INQ/oR6O98asjIyuK7OUtJdT19XM68czv+//0U85FogseNwLxlc1mGKyqpQpOG1vpOIA74XSkVo5TappQ6iLNVVazW+m9lFaQQwimoSiDjhvTFbDIBcDExmZ8XreBCQlK+fY0nT+IIrAJAVr8BWLteV6axisrJrWFElFJ+QAhwQWtdric1ljoN7+Kt9+DwsZPMW7meS/9/TSYjg3p0oW3zJnn2M8aeJnDas6Q+O61SjprrrT//SzxRpyFjT1Vg3v4fBrz7HhyKOcHCNZHYbH9NFzu4ZxfaqaY5/TuKUvXRB8nu2JmMv911RbMFlife/POHclYRLoQo35o1rMcDt43NM6nT8sht/LhgGdHHTxV5rM/C+fj99D1V/+9JgscOh8zM0g5XVBLujD3VuSwCueyaA5VSn7u+76mU+tr1lX/KMyG8WM3wUCaOGJgnccTGXeD35WtZt20Xhb1J8P/fjJzvrU2bga9vqccqKgd3njT+oZTa5BpCpNQ/tJVSTYFOwKWG5fcB9wP/A24u7esLUdH4+/pyy6jBdGnTAmOu10xb9vzJio3bC0wciT/MIvWJp7DVqSsDIIorUmzS0FrfAgwHHMCvSqnvlVL9SyoApdTjSqkFrq9/aa0Pa62n59rFpLXOAGKBWiV1XSEqEx+LhX7XdeTem0bTqO5f/012HzjEH2s3Ybfb8x7g50faM89zMXJ7gQMg+n/4Hobk/C2yhHC3TqMGUB8IA84DE5VSX5REAFrr97TWo1xfrxawS5pSyhdnwjhTEtcUorKqEhDA2MF9UY3q56z7MzqGuSvWkZFZQMNHf//8qz7/jCovP09In274LFtcmuGKCqjYSZiUUpuBNOC/wPNa60zX+iVuHNsNeFNr3V8pZQQ+AdoDmcA9WuvDbsQ4E5gBWHC+phJCFMFkNDKiXw98fSzs0dEAHDlxmm/nLqZT6+aEhwRTr1aNAltYGc6ezRkA0XT6FOa9e8gaUm5HDRIe4M7YU1211ltzLffTWq8p7sRKqaeBKUCq1rq7Umo8MEZrfYdrePVntNZjrzH+fOx2u8Nmq1zNiAtjMhnzNLf0Rt5+D4oqv8PhYNn6bazbuiffts5tmjN2SO/8icPhwPDLL5iefBwiamDdvAV8fEoj9BIhP/+SK7/FYrq26V6VUn2AVsATSql3LsUITAXauBFDNDAe+Na13BtYDKC13qSUKja4q2GzObym3ba3t1EHuQfFlf+6tq2pVqUqyyO35nk9tT3qIAaM9OnSPk/lOQDDRmNY1xPjhQvY0qyQZv1rW0oKxnNnsTfO24nQU+TnX6L9NNzar6g6jXickzBdqk+ohbNO42l3Tqy1ng1k51oVBCTmWrYppWSOciFKmWpUnztvHEX39q3x9bHkrN8WdYA5y9aSbbXmO8YRWh1bs+b51ge+8TKh/Xvg/8E7kJ2db7uo/Ar90NZaRwFRSqmZWuvYErhWEpA7lRm11vl/W4UQJS7Az5dendvRvUNr5q1cz5ETpwGIORXLb0vX0KND60LrOS4x79iG/38/w+BwUOWVF7A1bkrWqDFlVAJRXhQ1NPos17c7lFKnXV+xSqnTV3mtDcAI17m7A3uv8jxCiKtkMpkYO6hPnomdTp45x6+LV/H5r/M5eLTw6WIdfv5Y27YHnAMgZo0cXerxivKnqCeNCa5/S6pvxBxgiFIqEjAAd5bQeYUQV8BoNNKnSwfMZjORO/762y0pJZX5qzbQ4GA03du3pm7NiDzH2Vq1JmHxSvxnfkrmyNFw+VNJdjZYLIjKzZ3WU32BAJxPJR8C/9Za/1AGsV0VGbDQu3j7PbjW8p+7EE/UoSP8GR2Trx9HlzYt6Na+NQaDIU9dSGGC7piEIyCAlJffwFG9+lXHdCXk51/2Axa6UxH9FjAJ+BjoBfwClNukIYRwX0T1EAZW70yvTm1ZtXnHX3OP46wo3xZ1AICOrZrTt2uHnHk8LuezYB6+i+Y7v1+9govrt+IILZvEIcqWOz3C04GzgFVrfQZnayohRCXi6+PD9X26c/eE0TSsUzPf9p37DzJr8SoyswqeTseyOTLn+8zhoyRhVGLuJI0kYDnwi2v618JryoQQFVpwUBXGD+3PkF7XUSUwIM+2U2fj+OWPlSSlpOY7LvXlN0j84VeyO3cl9fmXyipc4QHu1Gn4Ak201vuVUm2Ag+V59j6p0/Au3n4PSrv8VpuNHfs067btzlnn7+fLyP49aVA7/xMJDkf+CnKbjSpPPUb6Xfdha9uuROOTn3/5nIQpAhillHoeZw/vf15jbEKICsJsMnFdu1YM6dUVoysZpGdkMmvxKn6Yv5QjJy6b7KmAfh7+//0U/++/IWRoPwLeeq0swhalyJ2k8SvO3txnc30JIbxIO9WUicMHEuDvl7MuNu4Cc5atZXnk1kIneyI1lYB33wbAYLNJk9xKwJ3WU8la6+dKPRIhRLlWt2YEU8YOY3nkNo6cOJ2TKHYfOEzMyVhaNGlI3ZrhhIcGE3hpyPXAQBIWLqfK3x/FmBBP2tTHPFgCURLcqdN4F9gM7MQ5ERNa64OlH9rVkToN7+Lt98BT5U9Lz2D5xm0cijlR4PbWzRoxsHtnfC49WdjtGM6fxxGRt8MgKSn4rF7p7F1exBAmhZGff/ms0+iAcx6Lz3DOa/HZtYUmhKjoAvz9GNW/J03q1ylw+75DR/lh/lISk1OcK4zG/AkD5wCI1e6aTNDtt2A8fSrfdlH+FPukAaCUqgY0AI5orVNKPaprIE8a3sXb74Gny2+32zly4jTpmZmcuxDPyTPnOB//12DWAf5+DO7ZFbPJSFJKGs0b1cPXYsFoNGLesY3g4YMwuD6Dkj6eSebEW67o+p4uv6eVyx7hSqkbgedc+/6ilHJorV+59hCFEBWd0WikaYO6edZFHTrC8g1bsdntpKVnMG/FupxtyyO3YjaZCA8NZnTnNvjefhcBX/+PrAGDyJxwc1mHL66COxXhTwLdcU6g9AqwzfWvEELk06ZZY4KrVmHu8nVkFNCD3GqzERt3gZmL1xDQsQ/h/qHUbN+WXjhHMr3EEH8RR0Ag+MogFOWJO3Uadte84A6ttQPI3x1UCCFyqVszgltHD6FxvdpUvaxneW5pGZkca9iUzYnpRO7ci93umrrU4aDqY1MJGdQb86aNZRS1cIc7TxrrlFI/AnWVUp8BW4s7QAghQqsFccOQfjnLm3ZFEbkzCqDAfh2bdu3jUMwJxgzqQ821q/BdvBCA4HHDuRi5vdxMMevtik0aWutnlVLXAzuAA1rr+aUflhCisuneoU3OUOuZWdnEXYwn22pl7dZdOZXnFxKS+HH+MsaeOUq7gEBMaalkTLpdEkY5UmTSUEq1AybgnBv8JHCgLIISQlROl6aT9fWx5EzyVDsijE2797Fr/yGsNhsZWVn8HFqHxY8+x7h1Swl+dlr+D6qCxrgSZaKo6V4nAl8AJ3BWgicDs5RSY8soNiGEF/D18aFf147cPHJQnvqPxOBQvh59C3/s/jNX7/ND/DBnEcaRQ/H77mtn8hBlqtB+Gkqp9cAwrXVqrnVVgXla6wFlFN8Vk34a3sXb70FlK7/VamP/4aPsORjN2fMX82zz9/MlPSOT3uuWMWrhLADsQ4dy/ptfMBjdadNT+ZS3HuHW3AkDQGudDNiuMTYhhCiQ2WyiXYumTB4zjI6tmufZlp6RicFup9OOTTnr1vgE8dH3szl68nRZh+q1ikoa9qs4RgghSkT/6zrSqG7tPOscRiOfPPR/rO43jJN16rOsx0Cysq38sXYTaRmZZGRmYbMX9tElSkJRFeGtlVKXzwVuAFqVYjxCCAE4e5uPHtiLXX8eApzzmVvMZjbt2sfi4eMx2mzYXXOWp2dk8ukPv+GTmcGg9StIfGAqXXteV+ic5uLqFZU0bipkvQxYKIQoExazma5tW+ZZN6JfD/YejMZqtZGRlcGOfYdytg1bMpdekSuJ37qeyAcep/Njj2CSxFGiCk0aWus1ZRmIEEK4w8/XJyeRBAcH4GPxYdOufUScPU2PjasACEm4SPLRGJZFbmVY7245TX3FtXOnR7gQQpRbvTq1IywkmH0Ho4l65gXUR+8QU7MOuzpcB4eOcjEhibaqCapR/b/m9xBXza2h0SsSaXLrXbz9Hkj5Cyh/XBxr1m9mW3xyntUhqck0qxlGs/79qBlevQyjLD3lrcmtEEJUPOHh9Bx9PfVr1/hrncPB6FnfMurxezn+j3+weNV6EpKSC5/bXBRKXk8JISodi9nMhGEDOHchnuOxZ8n+6Uda/bkbgJGLZvFxw6b87+gJjAYDEdVD6NCyGa2aNpK6DzdI0hBCVEoGg4EaYaHUCAvFZBtH+prF+B88wObr+nCifiMA7A4HZ85fZPG6zUQdOkI71ZTG9erg6yN1H4WRpCGEqPRsHTqRsmoDtv/NwNGzP/VPn+Pc+ficSaLM2VmcPBPHyTNxhIcEc/PIQeijx8nIzKJTayX9PXKRpCGE8A4WC+kPPEwToEk756q09Aw2bttNr8fu4WyNOiwcOYE44KPvZucclpGZRd+uHTwScnkkFeFCCK8V4O/HqL1bqHfyGF22R/LoB69gsmbn2WePPkxWdnYhZ/A+kjSEEF7NfGB/zvfHR43HZs5bn5GZlc2+Q0cBSM/MJDvbWqbxlTfyekoI4dWSP/iUzJFj8Pv6f9R64w3uSE0nPTOLE2fOErljLwDrt+/mROxZjpw4jclkYnjf7jRtULdU4zpy4hQrNm6nVnh1Bvfsip+vT6lez12SNIQQXi9r2HCyhg3HAFT3cX44R1QP5sDuKIZ//RnLB4/ikOsJw2a3M2/lem4Y0jffKLwlxWa3s2zDVlLS0klKSSXmVCwtGjekV6e2+Pv5lso13SWvp4QQogA+Fgt37txA26gdPPrBK/RZuzRnm8PhYM2Wnfk6B8acimXV5h1cSEi6pmtHHz9FSlp6znJmVja7Dxxi+z59TectCZI0hBCiAIazZwn77WcATHY71vAa1K9dI6cD4IWEJHb9eQibzTkvXWpaOnOXr2PHPs1PC5dxMfHqEofD4WDn/oKTQ0JScoHry5IkDSGEKICjRg3il60lu2s3sgYOpudbrzPx+oF0adMiZ5+Vm7bz6Y9zWB65jT0Ho7G6EkhGZhZzlq4hM+vKW11FHTrCyTNxgLODYv/rOuZsKw+tuCRpCCFEIWyqBQnzl5A04wtwPWF0bqPw9bEQlBhP66idOa+OLlWaX5KQnMJ/f5nL3oPR+V5jFTTmlcPh4MCRY6zcuD1nXfsWTakRFpqzfDVJqKSV24pwpdRA4Dat9T0FLQshRJkwGnFUC85ZDPT359aRQ6h6+y3U2byBqNYdmDv2VpKDgvMdmpmVzdL1W7DZ7HRo2Yzk1DRWbNzGsVNn8PXx4cZh/QgPDSHuYjwrN+3g5JlzOceGVguib9cOxCf+9UqqPDxplMukoZRqCnQC/ApaFkIIT6q1YQ3VNm8AoM2+XWwZMDwnadSpEU5otSD2HozO2X/z7n3UrRnOglWRXEhIBMCans43vy/Gx2LJlwxCgqoybnBfLGZznnGwsspBH5FykTSUUo8Dg12LG7XWrwLTlVLfAWitD+deFkIIT8ru3Yf0yX/D/7uvSb/9Lno//jCJqyNJSkmld+d21I4Io16tCBat2QhASlo6X8/5o8Bz5U4YRoOBNs2b0KdL+5x+GbknjpInDRet9XvAe56OQwgh3OEIDiHlnQ/JnHAz1tZtCAmqyqTRQwFn5bXx9Cla1qvNuTYt2BZ1wK1z1q9Vg4E9ulA9OCjPeh/LXx/TWVnZOBwOjw7hXupJQynVDXhTa91fKWUEPgHaA5nAPa6nCCGEqHCye/bO+T7ng9xqJej2WzFkZdLjtekcCPDP0+eicb06DOzeie/nLyU9I5NmDesxvG93LOaCP45NJhMmkxGbzY7d4cBqsxW6b1ko1SsrpZ4GpgCprlXjAD+tdQ+lVHfgP8DYwo7XWk8ualkIIcob/5mfYtmzC4Dakydyx8adpPgHElqtKg6HA6PR2Wj1b+OGk5SaRs2w0GKfHHwsFtJtmQBkZ1srb9IAooHxwLeu5d7AYgCt9SalVLHz0V4pk8lAcHBASZ+2XDKZjF5T1sJ4+z2Q8pe/8hurV8MRGIghNRX7M89So2UTahSwX3BwAHVwb65yf18f0jOcScPXz5xTZk+Uv1SThtZ6tlKqYa5VQUBirmWbUsqstS6xJgE2m6PEJlov70pyUvmKytvvgZS/HJb/5tsx9uhHwKcfknLXg3B5fHY7GK+si5wp1yRQ5y8kYTI4P7pLsvzh4VXd2q+sO/clAbkjM5ZkwhBCiPLAXr8BKa9PB0veYdYNyUmEDOyN74/fQQEd/AqTt9mtZ1tQlXXS2ACMAHDVaewtenchhKg8Al97CfP+KIIee4iqD9zl9nG5m91mejhplHVtyhxgiFIqEjAAd5bx9YUQwjPS0vBZsSxnMWvYCLcPzZ00PD0JVKknDa11DNDd9b0deKC0rymEEOVOQAAXV0US+NZrmI4cJvOGCW4fmudJw8PjT5WLzn1CCOEVAgNJffFVZ2X4Zc1sjadO4v/l56Q++TQE5G0R5eOTq4Ofl9VpCCGEuLz1lMNBlf97koAP3iG0b3csmyLzbPbNPZSIh580JGkIIYSHWdauxnfpYgBMx2NwkPcpJM/4U1bP1mlI0hBCCA/L7tufpPc/wR4cTPodd2Pt3iPPdp9y9KQhdRpCCOFpBgOZt04ma9BQ8M8/A0S1U8cJjr9AQkh1r2tyK4QQohCOiIj8K61W1LR/0vpINEuHjuVkjZvLPrBc5PWUEEKUY/4zPiHwwH58szK5fvFv2GJiiE9KLv7AUiJJQwghyrHsbt3Jaq4AWD54NCf9q/DNnD/YvHs/Npu9zOOR11NCCFGOWbtcR+LKDcS+8TrrazdxrrPZWL99N4ePn2BIm+ZENGpUZvHIk4YQQpR3Pj7Uen4at44bQUT1kJzV8SdOUWdof5In34o17nyZhCJJQwghKogaYaFMGj2Uvl07YDaZGLbkd4IT42m8dCHGof05fTau1GOQpCGEEBWI0Wika9uW/G3MUGrasnLWr+wzhHkrN+C4giHXr4bUaQghRAUUHBJMtXVriJr+AVnzfmd3+65U97UUO3XstZKkIYQQFZTBYKDGvfeSPuV2xp6No07N8FK/piQNIYSo4Pz9fGnaoG6ZXEvqNIQQQrhNkoYQQgi3SdIQQgjhNkkaQggh3CZJQwghhNskaQghhHCbJA0hhBBuM5R2l3MPiAOOeToIIYSoYBoAxfYOrIxJQwghRCmR11NCCCHcJklDCCGE2yRpCCGEcJskDSGEEG6TpCGEEMJtkjSEEEK4TebTKKeUUt2AN7XW/ZVSTYGvAAcQBUwF6gM/aa27u/bvDXwJ3Ag86dq22O2kDiAAAAZjSURBVBOxX63iyvz/7d1rrB1VGcbxf5vTai3WhovQhJAmBh8aw9XQgtD0qNhCGq02JhQQaEiApl5owGAhGqFcEtD4QRAUNQFD762RCimYUA630lAQarThMWmtitLWQ6FFbSm08GGtHYfp2TC7Z+/x7PT9fdkzs+ay1pw1s2bNnHnH9n5J3wemA28D82w/K+lecnkl9QCLgH7gB8Dixj7qRqV9ci/wZeBo22/m9NOA54HPAuOBE2zP/z9lt20q1v+/APNt31ZYbhUwJi+3hbQ/9tSc/YMyyPpfW72InsYQJOla4BfAh/OkHwHftT0ZGAbMKM3fC9wDTLf9hxqz2jZVypwPhCnAJGAW8JPSOkYAy4DNtueSDriuNcA+AXgFOK8wfhGwuc58dVoL9X8T8NXCcocDx9eY1bZpQ/2vrV5EozE0bQJmFsY/DTyeh1cD5zQSJJ0D3AWca/vPhWXmSnpU0uP5qmWoq1Lms4Hf2X7H9t+AHkmNN1g/BPwaeLF0RXWUpFWS1kn6XmeL0HblfQKwGLgAQNJw4DRgfSH9zPx3Xy9pej3ZbLuq9b8f2C5pQh4/H1heWtfP8jGwTNKoTmW4DQZb/2urF9FoDEG2VwJvFSYNs924an4D+Fge/gRwC+nqpHxArLX9eeA24PYOZrctKpZ5DLCzME9xX/wYGA2Uv3l5GHAxcBZwnqST25z1jhlgnwA8C0jSaOBzwGOl9P+QTjDTgTvzCaSrtFD/IZ0sZ+XhGcBvSqu72/YUYAtweftz2x5tqP+11Yuuq1CHqP2F4Y8Cr+fh3aQu6beA8pXUE/l3LaCO57D9Birzrjxcng5wBzAVOFHS1wrzbLC90/Y+0oH1yc5luTarSCfIC4GFpbSn8pXodtIJ5oi6M9cBzeo/pEbiS5LGA1uB/xbS9tpel4e77Thotf5DTfUiGo3u8EJ+bgGpkXgyD//T9g7bD+ZpdxaWmZh/J5MepHWbgcr8NDBN0nBJxwHDbffnef5o+23SvdzbC7csJkg6LD8gnwT8qb4idMxC4BJgnO1NpbTTASQdQ+pl9dP9mtV/bP8bMKk3vai03EhJp+ThbjsOWq3/UFO9iEajO1wD3CjpGWAksGKAeb4NTJR0aR4/Q9IaYB5wbT3ZbKsDymz7edLB8wywkvRfNO9hezOpvMuBjwA7gKWkK80VtjfWk/3OsW1SNNLfDpA8Kv/dVwFXFm5xdLMPqv8LSY3Co6XpbwLflPQE8HHSg+Zu0XL9r6teRJTbEEIIlUVPI4QQQmXRaIQQQqgsGo0QQgiVRaMRQgihsmg0QgghVBaNRjgkSeqVtGSA6UskjXyf5ba2uJ1jJN11MHksrGOypKua5bnJMqMk3Sdp2GC2HUJZRLkNocD2rA+eq6X1bQXmHuzy+aR/A+kFr8+0sN3dktaSXva672C3H0JZNBohFDTCaQM/Jb0cNh4YB8y2/fvCfLeS4v58o/GiVA4et5TUgx8BzCHFB1pCClG9Oi/eA5xBCmtxLCl+2D5S0LorbRdjEH0B2Gh7rySA4yU9QgoDcbftX0rqA17K+R4GnJ8bq2XAw0SjEdoobk+F0NxfbU8jxbW6ojFR0g+BHttfL71ZO5EU16cRD2xMI8H2btu9pMZjC6n3sRn4OTAzB9X7BzC7lIdeoBjufgTwRdIb0N8pRDldm9e/FLg+b/M14EhJxQB/IQxKNBohNPdC/v07//vOwdHASaT4PWWrSeGsHwAW8N6gcw13kCI+3EMK+TCOFGyyjxRw8bjS/EcC2wrj62zvtb0b2EjqCQGsyb/lwHzbgMObljCEFkWjEUJzA8XY2QZMAz4l6dxSWi/wiu2pwM3ArcVESTeRQl7flCf1Ay8DM3Iv4RYODGm9HRhbGD9VUk8OgT2BdEsL0vcXIIWALwZlHAv8633KGEJL4plGOJRNlfRcYfzCKgvZfkfSZcAjkibZfjUnbQCWSppHekaxoLGMpNOB64C+3Ksgp18FPJS/b7CL9OC6qA/4CvCrPL6H1KMZC9xge0d+1jFb0tWk7yZcnLc5Fng9R4INoS0iYGEIQ1huTNYAU23vbTJPHzDH9kul6XOBXbbv73hGwyEjbk+FMITZ3g/cSIv/tps/yHUWB35jIoRBiZ5GCCGEyqKnEUIIobJoNEIIIVQWjUYIIYTKotEIIYRQWTQaIYQQKnsXptKr+qJJ6EwAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "A, B = 1.62e-5, -0.8164\n",
    "title = \"AP85 Chr1A (RE=AAGCTT)\"\n",
    "label = r\"$Y=1.62e^{-5} \\times X^{-0.8164}$\"\n",
    "filename = \"AP85_RE_AAGCTT.Chr1A.distribution.txt\"\n",
    "plot_distribution(A, B, title, label, filename)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 46,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAESCAYAAAABl4lHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3Xd4VFX6wPHvzKT3QBIIBAj1UEIHaSpIkyId7HVtP9euu65tdd3VtbGCBVfZXXsXRRERBVGQKr1z6BBqKAnpM5ny+2MmQyZ1gEwmybyf58lD7j3n3nkPA/POvefccwwOhwMhhBDCG0Z/ByCEEKLukKQhhBDCa5I0hBBCeE2ShhBCCK9J0hBCCOE1SRpCCCG8FuTvAETtpZRKBfYBS7TWA0uVvQfcBCRqrU9Wco73gC1a66mV1LkZmKy1vuIcYvsv8JnWemE5ZR2AZ4G2gAPIAp7QWi91tWmL1jrKy9d5DGc7g4CPgGe01o4S5cOBl7TW3UodFwwcBDZorUdW8RrdgPu11re4/r6GASdcxUYgCnhLa/2Sq/5+wAwUlDrVH7XWy71s1xogAuhUsj3elCulxgMPAo0AkyvWf2itfyhRp7L34FHgalfVNq7jz7i2/w3c5fq9sev8h13bz7va3Vlr/Q9v2imqnyQNUZVCQCmlWmitD+DciAQG+DMorfVt5e1XSingZ+AWrfWPrn1DgLlKqQFAnrevoZQaBVwJ9ARswI/ANuALpVQ48ARwN2c/1EqaCGwAeimlOmitt1fwGkbgf8DYErunlUyySqnmwHal1Byt9Q7X7uu01mu8bUup1+wDhAEW4HJgvrflSqk7gAeAK7XWW1z7ugLzlVJjtdarq3oPtNYvAC+49v8KvKG1nlUihLddZX8DErTW95SK726lVDet9Ybzab+4MJI0RFVswOfAdcA/XfsmAt8CDxdXcn2Y3Oeqfxy4R2u9s+SJlFKXAC/j/AZrAZ7UWpf+wJoMvAiMAg7h/ObZFmgI5ADXaq11BR82AI8C7xZ/WAForX9WSl3D2W/mJqXUW8BFQCzwiNb6K9eHVD+gCbDRVf8TrXWeK7Z3geuBL3B+mEbivAop/nsp6S7gM2APcD/wf+XUAWdS2qe1Li/xFEsBDK72V0op1RH4pJyiV7XW75aIbS5wEmcCmF+qbrnlSqkQnG0dUZwwALTWG5VSd+K8KgDv3oML8T/gaWBCNZxLnCNJGsIbH+C8NVP84XgTzg+ThwGUUoOBR4B+WusTrttN3yilOhWfQCnVEJgFjNVar3KVLVZK9S5R5xqc394Haa3TXQkkS2vdz1X+FnAPcG8lsfYC/lJ6Z/GtE9ftqTBggdb6/5RSE4CpwFeuqi2ANK21VSk1H+c35mKHcH6Ao7X+xtXGQaVfy/XB3Q+YBKx1tfMJrfWpcuKdjPMDuqQHlVLXAzGun6XA6FKJ5WOlVMkPYLPWuo/WehvgcausVGwNgKuA3jiTwvNKqY6u46oq7wgYyrvC0VrPKbFZ6XtQDX4C3lNKhWutqyMJiXMgSUNUSWu9VillU0r1BDKAaK31FuddCABGAJ9rrU+46r+nlHoVSC1xmj7Abq31KledrUqpZcAgnPe8e7vO84DWOt1VZ5ZSaq9S6l6c974HASuqCNdO1QM8LFrr4iSxAUgqUbZSa211/W50xVbMgPNKqip3AXNdSeKUUmofcAfOe/KltQeml9o3TWs91XUb8HOc9/F/LVWn3NtTXlxp3AJsK3FraQHOK6E7XfUqKzfg+feBUuo3IBoIB37XWt+Ad+/BedNan1ZKFeJM8Duqqi+qlyQN4a0Pcd6aOeH6vSQTzttNJRmA4FJ1Sne4Gl11LDg7Sq/B2V8wV2u9Xyl1F84P2zdwfhCeBlpWEedKoC+lvr0rpZ7CeatoGVBUosjhirVYbonfD+K8VVWsCc6rjQq5PuhvAMyuDmtwXi3co5SaqrUuKnWIgwo+YLXWeUqpG4DtODueX6nstV3HVHiloZQy4LxN1qBEbBHAIKXU4zj/fisr3w4YlVJpxUlFa32J69w347xqgireA631x1W1wwtWvEvgoprJkFvhrY+AKThvXZT+JjsfuFoplQiglLoFOAXsLlFnBdBeKXWRq04n4FLOfoPepbVeBLwOfODqIL4ceE9r/T9AA2M4e9+8Ii8Dt7tGNeF6rRE4vy1vPJcG4+y3uU4pFamUCgVuBr6p4pjrcLa9idY6VWudCrTCOQJqSjn1NdC6opNprTNx3gZ8RinV9BzjL20YzquqViViawIcxXklUWm51roQ522nj11XNAAopZJcxxZ/iFfne1CGUioW5y3Ggxd6LnHuJGkIr7jup2/H+eF+ulTZAmAasEgptRVnn8cVWmt7iToncX5ovq6U2owz8dxSurMceA5nB/OfcfY13KmU2gT8BqzDeZuqsjh3A1cAf1JKbXLF8xdgTMnOWy/b/B3wNfA7sAVn/8QHVRx2F/CK1tr9LVhrnQW8hvNqobRZOG/LVRbHx67XLjls+WOl1IZSP+WOKCsV20ytdfHwVly34v6JcxRYpeVKqWCt9X9wdnS/6nrNXcAiYC/wB9cx1fYeVGA4ztt/5mo4lzhHBpkaXQj/UUqZcCaE0h3dogJKqUU4+742+TuWQCRXGkL4keuK5HbKH7YrSnGNdvtNEob/yJWGEEIIr8mVhhBCCK9J0hBCCOG1evecht1ud9hsgXHLzWQyEChtLS2Q2w6B3f5Abjv4rv3BwaaTQGJV9epd0rDZHGRl5fs7jBoRFxcRMG0tLZDbDoHd/kBuO/iu/YmJ0Qe8qSe3p4QQQnhNkoYQQgivSdIQQgjhtXrXpyFEXWazWcnMPIHVWnr+R0/HjxsI1GesArntcOHtDwoKIT4+EZPp/D7+JWkIUYtkZp4gLCyCyMjGGAyGCuuZTEZsNnuF5fVZILcdLqz9DoeDvLxsMjNPkJCQfF7nkNtTQtQiVquFyMiYShOGEOfLYDAQGRlT5ZVsZSRplBbAl72idpCEIXzpQv99SdIoIWTeXGKunYzxULq/QxFCiFpJkoaL4UwWUX95iNCfFxB/SR+Cf1vs75CEEKLWkaThErJoIcaM486NiAisaZ39G5AQQtRCMnrKxTxhMrYmKUQ/eDd5jz2FI76Bv0MSQtRzBQUF/OtfLxAcHEz37j0ZPnyku+zYsWNMm/YisbFxNGvWgmHDRjBt2ovExcWTktKcG2642X2Ou+++nVtvvZMBAy4ps13d5ErDpchq5WDLNhyYuxDLmHFlykPmzMa0ZbMfIhOiZu3fv48pU8ZitzuHddrtdh588G5++GHuOZ9r5crlXHPNRK66ajwffvjeOdepqKyi/Wazmdtvv5GbbrqG66+/kv/9720Ajh8/xr333sl1103m+uuv5IsvPnUfM3nyGG688Spuvvlabr31hnNuY0mzZ89i6tQX3NszZ77JP/7x1wrrL168iEGDhvCXvzzJ0qVLPMrS0w/Qr9/FPP740+zfv9e9/eSTf2P//r3ueh9//D6DBw+tcLu6yZUGYLYU8c6s78gvdC45bDIZiQwPJyjIRIsmjbkkIZaG992FwWIh/74HyX/wEQgN9XPUQvhGampLmjdPZfny37j44oG8/fYMmjdvwciRV5zTeWw2G6+88iLTps0gKakRt912IxdffCktW7byqk7FZS0rPCYkJIRXX32LiIgIrFYrd911K3369Kdx42TuuedBlGpPfn4ef/jDDfTu3ccdy2uvvU1cXFyFbcnMPE1oaCgREZHufYcOpZOS0syj3siRV3DttZPIybmbTZs2sGLFUt566x327NnN22+/4VH3scee4sSJDFq3di57bzR6fodv107xwQfvsGjRAi6/fJR7+5dfFjB8+CgAVq9eRWpqSywWS7nbviBJA8grKHAnDACbzU52bh4Ap7Oy6fz+DJrkO2eVDP1+rjNpCFGPXXXVtXz++SdYrVY2b97Ia6+9dc7n2L59KykpzWjaNAWAoUOHs3TpYo+kUVmdisoKCvIqPMZgMBAREQGA1WrFZrNiMBhISEggISEBgIiISFJTUzl5MsMjlsqsX7+Ob76Zxcsvv0poaChz5sxmyZJfmDr1NY96YWFhDB06nJkz32TlyuVMnz6D0NAwWrduw0svTS9z3sTEJDIyMmjbVuFweD6w9/3333HrrXfSrVsPnnzyEc6cOcOtt95Jz569eOyxPzF69FjWrVtDYWEh+/fvJTQ0lNTUVpjNZvd2v34DyiSjCyVJA4iPiWZw355s272PYydPlymfN2oS4QX5tDiwh/k33klXDIT5IU4haspFF/XljTem8dZbM3jjjZkEBXl+VPzxj7eRn192eu67776f3r37AHDiRAZJSY3cZYmJSWzbtsWjfmV1Kio7ceJEpee12WzceusNHD6czoQJU+jUKc3jNY8ePcLOnZqOHZ37DQYDDz10N2Bg3LiJjBs3sUy7Bg8eytGjh/nb3x7nssuG8v33c5g2bUa5f3ejR4/juusm88IL/3IntooMHDiYadNeYsWKpQwYcCkAW7ZsYudOTd++/XjnnZksWDCfxo2buLcXLvyRxo2bAHDnnXcDMG/ed8TGxrn7MIq3qzthgCQNwPmPpnvHdnTv2I6MU5ns3HeQhvGxnMrKZtXGrZxMbMzMOx6m+cG9HDCFs/aruaQ2TSYuJoqenRThh9Kxt2rt72aIeijipX8S6bpHnvenR8l/5HGP8sinHifiLedtj9y/PUfBH+/1KI96+D7CXff8c6a+SuGNt3j92mlpXWjXTrm/oZf05pv/rfL48p6TLf1gWWV1Kiorb96lkuc1mUy8994n5OTk8Pjjf2Lv3t20auW8BZSfn88TTzzC/fc/TGRkFAD//vf/SEhIJDPzNA88cDctWqTSrVuPMq9x3XU38fTTjzF16gt8/vk37iua0t599z/ExcVjs9nKLS8pPDycxx9/2mNfWloX0tK6APDssy95lD377EvlTiMyatSYSrerkySNUpIaxpPUMN69ndo0maycHPalH2GnK2sXFJrZvmc/AOlzvuOPb75ExqixBE17FUdcfHmnFaLO2b9/H6NGjS23zJsrjaSkJDKKh7HjvHJISPBcGK6yOhWVeXNegOjoaLp378nKlSto1aoNVquVJ598hOHDRzBw4GB3veJj4+MbcOmlg9i2bWu5SWPjxvXs3buHSy8dxLvvzuShh/5Sps6nn36ExWLh739/nnfemenxOvWFjJ6qQkrjRNLatuKKywYwZvDFRIafvTFlslqZPOsDjA47jb//hpxbbmKT3sO+9CPukSdC1FX79u2lVQVX0G+++V/ee++TMj/FCQOgffuOpKenc+TIYYqKili48Cf3LRhv6lRU1qFDpwqPyczMJCcnBwCzuZA1a36nRYtUHA4Hzz//d1q0aMnVV1/vfv2CggLy8/Pcv69evarcNu/cuYMXX3yWF154hccff5rs7GxmznzTo87atauZN28OTz75N3r06EVeXi67dulz/Wuv9Qz1bYrhoiKbw5dLQRaaLWzdtZeNO3aTn5HBxK8/pOumtViCQ5j+wFOcbuj81tI4oQHDL+5DYoOKR2RcqEBe9rK+tv3YsQM0btyiynq+nun1+PFj3HPPnXz55bcXdJ4VK5by6quvYLfbGD16LDfddCsAf/rTfTz66F9JSEissE5Fx5tMRpYuXVLuMbt37+K5557Gbrdjt9sZPHgYt9xyOxs3buDuu2+jdes2GAzO78p33vlHmjdP5fHH/ww4+0KGDbvc4/WLbdq0gcjIKPdIJ6vVyrx53zF27ATA+UzFvffeycsvTyc1tSXg7FdYv34tTzzxtwv6OyytOt778v6dJSZGrwV6VXWsJI0LsHN/Orv2p2P87lsi8vNY03uAR7nRYKBj61R6dm5PQnz1J4/6+sHpjfra9tqSNGqzQG47+D9pSJ/GBWiX2ox2qc2wDOjN5p176JKVTW5+AfsPH8Nut9N+y3p6/286302+iWHXTqFJUtkORSGEqEskaVSDkOBgenZq794+lZXNkp8WMf6bj4nJySZ16lN8eiqDQ1Om0KpZUxrERnMy8wzxsdEEB8lbIISoO+QTywcaxsVwVYiN6FxXh1xoKAcaNWXHmo38tmaju150ZARXjhxCXEyUv0IVQohzIqOnfMQyfhJZ3y+gsHVb5k+5mYIS0w8Uy8nL53+zvuPjOT9xJifXD1GK2qi+9TOK2uVC/31J0vAha6+LyFmykkuefYbhF19E2xYpGF0PIXXZuJomhw8CcOzkKX5a+rt8WAiCgkLIy8uWfwvCJ4rXCA8KCjnvc8joqRqWnZvHqQ0b6HntBAwWC4svHc7CoVdgCwomrV0reqW1p0Gsd2tE19cRRN6or2232axkZp6ocg3nip6MDgSB3Ha48PYHBYUQH5+IyeTZOyGjp2qpmKhIUt5+DVNhIQDd92xn4VDn7KFbdu5ly869JMTH0rRRIvEx0XTv2M4n88eI2slkCiIhIbnKevU1aXojkNsO/m+/JA0/yHvmOQzZ2QSvWoF5xkwaHM3iRGaWu/xk5hlOZp4BILeggIG9u/srVCGE8CBJww9srdpwZvb3BK1ZjfGiPlxrtbHrQDo796ez/9BR4o4d4WSicxbPtVs00RERtGmRQkxU2c50IYSoSdKnUctYVyyn8YRRrOt2Ed+PnkJ+5NnhuKplc0Zc2pcgkwnw/2WqPwVy2yGw2x/IbQfftb/e9GkopR4AugFtgY+11m9WcUjdZbGQ+JcHMdjt9Fy3klCHg4+u+oO7WO87yO4DhwgPCwWDgZ5pbenarh3BwbX+bRRC1BO1vodVaz0duAPYCpz78mF1iKGwAGuHjgA4IiKIf+0NBvftSfMmZxecsdnt5OYXkJuXz+JVG/n2598CeiSJEKJm1bqvqK4ri+JV0VdorZ8DrgG+1lrX61nKHDGx5Lz9LuYJUzBkZRKiFN2Bbh3asnzdZlZt2obDbocSw3EPHDnGzn3pqFbN/Re4ECJg1Lqk4bqyKL2Y7iXA7X4Ixy8sI0Z5bBsMBgb07ELf/ZqI99/n+LMvsux0Ljv2HgDg+8XLyS8spFuHtl493yGEEOerRpOGUqoP8KLWepBSygi8CXQFzMBtWuvdFRwapLUO6HswhqxMkp56DOOJDKLHDif2vx9yMDyM/IJCHA4Hi1auJSQ4mE5tW/o7VCFEPVZjSUMp9QhwA5Dn2jUeCNNa91NK9QX+BYwr71it9c3evo7JZCAurvy1e+syw4LvMZw66dyIjSV+6EBuslj5/LtFnD7jnBhx/m8rycnPY1DfboSFnv80AXWByWSsl++ztwK5/YHcdvB/+2vySmMPMBH40LV9MTAfQGu9UilV5VAvb9hsjvo5HG/YaILmLST6gbvJe/SvWAihSVIc1465nI/n/EhmtjNxLFu7me17DjBh6EDyCwuJjYokMiLcz8FXPxl2GbjtD+S2g0+H3HpVr8ZGT2mtvwKKSuyKAc6U2LYppWpdH0ttYu3Ri8xFy7CMusK9LzQkmLFDLqGv3kzTQ84+jtNZ2fxv1nd8OncB73/zg8ygK4SoNv4ccpsNlExtRq211V/B1BnlLNqUlJ3JuM/e4Z43X2DU/K8JKjqbmwsKzSxYtlqG5QohqoU/k8YyYBSAq09jsx9jqdOinvwLhoJ850OBB/fgKDWC6sCRY/ywZAWF5spnThVCiKr483bQbGCYUmo5YABu8WMsdVru35/HkJdH8PKlWP79H/oZw1m3dQf5hWZ3ne17DrA3/QjDB1xEu5byTIcQ4vzI3FN1mEeHmMNB0Lo1WHv2dpfb7HZ+Xr6GY0uXcSKpMQBGg4HxQy8lJTmJvQePANCuZbM693yHdIYGbvsDue0gc0+J6mIweCQMAJPRyMhQiJ/+DFt69GX2yEnkR0bx9YLFHvUuzetG784dajJaIUQdVevnnhIXwGwm5qF7MNjtdF6znEnff1FutTWbt2O12Wo4OCFEXSRJox4zWMxYO3cFwBERienFqeWuAphfaGbVhq0UWWXwmhCicnJ7qh5zRMeQ8+Z/ME+YhOHUKRp27crlkTH8smodTZISiI+JZu2WHWAwsHLjVtZv30nfbmn0Smvv79CFELWUJI0AYBk2wv17xzYt6djGOT+VYfZXXPTe63w+9hpON0zEbCli8e/riY6MQMkIKyFEOeT2VIAyZGXS4Mm/0EJv5eHX/kG3Q3vdZQuXrya/oJBTWWdYtXGrPFEuhHCTpBGggn9bfHYCxLh4LrnrDqIjnZOgFZot/LxiDR/P+Ymlazcxe8FieaJcCAFI0ghYljHjyZq/CGvHNHJfmkZIYiKD+vRwl+/cn+7uGD+Vlc3cX5axN/2wv8IVQtQSkjQCmLVbDzIXLsEycjQAbVuk0DAuFoCuG34nJX2fu+7O/enMXrCE9KMZfolVCFE7SNIIdCUmQCxeIbDBqRNM+uoD/vjmi4z6fpbHBIhrtmz3R5RCiFqiytFTSql4YADQEMgAlmqtc3wdmPCPti1S6Pj7IkJciaL9nh38aBjvLt+bfoRDxzJIaZzkrxCFEH5U4ZWGUipRKfU/YC7OxZPaAWOAn5RSM5VSjWooRlHDbFOnY7lkEA6jkZzpb9C2XWuP8q9+/JVNO3aTnZuH3W73U5RCCH+o7ErjaZzree8sXaCU6gA8Bdztq8CE/9hbpHJm1rcEbVhHTPeejAZ6pbXns+8XYrXZiD96iAXLndOONGucxMiB/QgKMhEeGurfwIUQPiez3NZhNT3b54nTWax7532uffkpNnXpxXdjriIvyrmOlslo5JorhtEooUGNxCIznQZu+wO57eD/WW6r7AhXSg1QSm1QSh1TSq1RSnWvlghFnZMYGc5V33yE0eGg28bVjJ3zmbvMZrezfnuZi1IhRD3jzeip14BrtdaNgZuBGT6NSNRaBosZWzfnsxz2iEh+HnOlR/meg4fLPARothSRnZtXYzEKIXzLm7mnsrTW2wC01luUUoF7XRjgHNEx5LzxNoUTJ2M8eZK0tJ4sWrnWXV5YaObwsQyKrDaKbDaOnzzNum0aq9XG6EH9ad+qhR+jF0JUB2+SRoZS6r/AIqAnYFRK3QGgtZ7py+BE7VQ0eBgA3YGWKcn8smo9e9MPk7ZlHbEzX+GrSTdwKsFzSO76bTuJioigUcN4goNlnkwh6ipvbk/tANKBNkA2sBhIBhr7MC5RR8TFRNOxdSrh+XmM+/YzWu3byQPT/07bnds86h3JOMnn8xbyxQ+LZB4rIeqwKr/yaa2fUUqNwfmcxhat9Y++D0vUJW1aNOX08UNE5Dtnw82PiKSoR0/aNmzIrgOHPOoeO3mKw8dPyMOBQtRR3oyemgFcB1iB25RSL/k8KlGnmEwm1J8e4u37nuRIcjN+ueEOxk8ey9ghl9AutVmZ+is2bMFsKSrnTEKI2s6bm8tdtdYXu35/VSm13JcBibqpQWwMQ+6+g+3jx9CrTUtMrmVlk5MS2Lk/nW7rV3KqYRLpzVtx8Mhx3vpsNteMHkZSw3g/Ry6EOBfeJI2DSqkUrfUh19Qh6b4OStRNDWJjaBAb47GvSVICDU6dYOLXHxFktbK8/2DmjxiPlRA+/HY+HVun0q97GnEx0X6KWghxLipMGkqpo4ADCAMmKKUOAM2AkzUUm6gHkhMbcuXCOe4JEFvv1diNJnf5tj37yTidyY3jR2IwGAA4k5PL6axsWjRtjNEoEzELUZtUmDS01sk1GYionwwGA9H/fYfCh+8ndMkvhH3wMQMMofy2ZqO7zsnMMxw4cozUpsmcPpPNx3N+wlJURPPkRowfdinBQTJEV4jaQv43Cp+zN2tOzudfU7BpA9au3ellt2MwGFiyegMAjY4d5qsffyUkOAhLkdV93MGjx/lhyUrGXDbAfRUihPAvufYXNcNgwNrVOW2Z0Wikd+cO3Dr5Cpof2MP9r/6Daz+eScjp02UO27U/nb3pR2o6WiFEBbwZcvunmghEBJ640BBunPclRoeDLpvXMua7z8ut98vKtViKZIiuELWBN1cao5RSpqqrCXFuDEUWTP0HAGANj8D04lSmjLiMKSMu44/XTiQsNASAM7l5vP7hLH78bZUkDyH8zJs+jQTgiFJqH87RVA6tdX/fhiUCgSMqmtxX38Q8fhLGkydo1LOHR/ngPj344dflOFwjqLbs2ovJZOKSXl39Ea4QAu+SxhifRyECWtFlQ8rd33X7Rtp/9CYfjJzMyUTnVGcbd+xi445d9OjUjoG9u7uH5C5ZvYEdew/Qv0dn0tq2qrHYhQg03tyesgLPATOB4chEhaIGGDJPE/2Xh0natpmHXn+OdnqLR/m6rTv58bdVWG02jp44xerN28nJy2fxqvWybrkQPuRN0pgJvAOEAEuAV30akRBA8PJlGLIyAXA0aEjLKZPL1Nm2Zz+ff7+Q1ZvOzqhbaLFwIjOrxuIUItB4kzTCtNaLcPZlaKDQxzEJgWX0GDJ//JWiLt3IfXk67bt1oXmTRoSGBHvUO3bydJmZdA8fO1GToQoRULzp0zArpS4HTEqpvkjSEDXE1rkLWT/+AiYTQcCUEYPda3GceG06S07lciC1TZnj0o9l0KOTquFohQgM3iSNO4CpOEdR/Qm4y6cRlaKU6g3cAsQA/9Jar6/J1xd+ZvIc7W0wGDDu3UPHV16gY2Ehy/sN4oeRE7EGh7jrHDqWgd1ul3mrhPCBKv9Xaa0PAS8C/wCe0Vrv83lUnnoCHYEUZIZdAUQ9/TiGggIMDgcdDh+gTatUbp08hvCwUAAKzRa27q7pf6ZCBIYqrzSUUs8ClwG/A/cppWZrrV/2VUBKqQeAoa7NFcACnB3xg4HRwPu+em1RN+Q+P5Vg7BgWLsTw33cZ7ZqepFv7tqzY4BxltXTNRponNyI2OsqfoQpR73hze2oEcJHW2u56MnwF4LOkobWeDkwv3lZKfQZcD5zAecUhApw9pRm2OXPJXboSa+ezD/r17tyBzTv3ELVnF8eSU/jku5/ootqQlZNLr7T2NEpo4MeohagfvEkah4Bo4AwQDBw/3xdTSvUBXtRaD1JKGYE3ga6AGbhNa727nMM+Az4ELMCfz/e1RT1jMHgkDIDg4CAmRgbR/rVn2dqpO9+0dIQ+AAAgAElEQVSOu5qVhWYA9h86yg3jRxATFemPaIWoNwzFo1EqopT6HWgBbMT5Td8CHAM4l+lElFKPADcAeVrrvkqpicBYrfXNrlFZj2mtx51fM86y2+0Om63yNtUXJpMRmy0wH2Qrt+2FhQT17oVB7wBgY5defHrt7e7ihvGxXDlqEE0aJdRkqD4h731gth181/7gYNNaoFdV9by50phy4eEAsAeYiPOqAeBiYD6A1nqlUqrKYL1hsznIysqvjlPVenFxEQHT1tLKa7shN4fIXhcRrndgjYhg3qhJHuWnMs8w87PvGD/0UlKb1u01xuS9D8y2g+/an5jo3ZLLVSYNrfWBC47GeZ6vlFKpJXbF4LzlVcymlArSWlsR4jw4oqLJnfYG5gmTMZ7IYHj/gWRl5wIOflm5DqvVir3IyrcLf2PS5YNIaZzk75CFqHP8uXJfNs6+kmJGSRiiOhRdOgiA5kDz5EYApDROYus/X6DHoh/4atKNfP3TYq4aPYRGDaVzXIhzcU5PP7k6r6vLMmCU67x9gc3VeG4hPDS0FTH2209Jda0U2HLrBhb/Ls+JCnGuvHlOYwpgAkKBl5VSL2mtp1bDa88GhimllgMGnE99C+ETwSuWY8p23g3Ni4ziQIvWmI9msHz9Zrp3aOd+MFAIUTlvRk+twnlF8BlwBfCT1npgDcR2XoqKbI5A6SQL5A7B82m7aesWoh+8m1+Hj+PnBmc7wkNDghk24CJUy+bVHabPyHsfmG0Hn3aEezV6ypvbTcUTFOZorc149kMIUWfYOqWRNf8XGt98k8d+s6WIo/+aSt5PP/opMiHqDm+Sxj5gDfCOUuppYJVvQxLCh4xGGic0ZGj/Xu7RUw1PZjD+649JvX4KIQ/fT8ahw1R1BS5EoPJmyO3NSqkorXWuUmoNzqe3hajTurZvS9f2bTl0LIPQqyYSbC0CIG/hQj5u35uunTswuG9PP0cpRO1T5ZWGUup1rXWua7MrcqUh6pGmjRJZevNdbG/fGZvRyKzJN2I3mVi/bSerN2/HUlTk7xCFqFW8eU4jWyn1AhAFpOGcwFCIesFgMNBt2GV84DDS+Eg6R5s0c5ctWb2BDdt3cXPLZBxdu2JzOAgPlVFWIrB5c3vqCaXUy0AbrfUg34ckRM1qmdKEO68ZT6G5CIMBPvxmPlabDYD4LRtJvvcGtnXsyjfjrqX38MF079jOzxEL4T8VJg2l1FGguDfQADRSSh0B0Fo3qYHYhKgxkeHhRIaHAzBm8AB+WbWO3FOZTPrqQwwOB522bsBqCuLrhATS2rUiOMifkykI4T8V/svXWtftGd2EOE+tmjWlVbOm/LJwCftatiXpxDHMIaHMGzUJs6WIRSvWckmvrkSEh/k7VCFqnDdPhHfDuU64+3+I1voPvgxKiNpg4OCL2dY8hXXrfqco/RBn4pzzVG3ZtZctu/bSsmkyYy7rT3BISBVnEqL+8OY5jfeAdcDnJX6EqPeMRiNp7VrR7OqraXzPPZiMnv9doufNIXzYILJX/+6nCIWoed7cmD2mtf6vzyMRohaLCA9j1MB+bNixi/SjGUTk5TJ2zudE5eUQPXYE6555iZZ33ObvMIXwOW+Sxn6l1KPAelwd41rrn3walRC1ULuWzWnXsjlncnJZ/rd/EFbonP8nNyqaX4IjaOnn+ISoCd4kjVBAuX7AmTgkaYiAFRsdxcGL+vNaQhKTZ33Az0NGk+kwUGA2y3Mcot6rbMht8Sp6d9ZgPELUCWltW7EyL59/3/UIDldfx9GMk+QXmtmwbRedlv2CGnoZEcMv93OkQlSvyjrCP3D9qYEdrp/i34UIaL3SFM2TG7kTBsDsBUv48bdVWPV2LvvoP7S4fgphD90H+YE7jbeofyp7TuNa159yq1aIUkJDQpgycjDb9+xn3uIVHmWjv5/lngCxaNlSCA72R4hC+ERlt6dWcPaJcA9a6/4+i0iIOqRJUkKZffOm3ARfvo/SW/h6yk2Mk6Qh6pHKOsKvrrEohKijYqIiad28KXsOHiYyPJyeaYoenRQzomNpcPggxxskcSrzDA3iYjAYDAAEbdqAtXNXcG0LUZdUdnvqQE0GIkRdZDAYGDfkEvILCokID3MnhuRGCRx0TXr43ux5pDZNZvzQSwhbvYq4sSMwjxhF7ouvYE+WadxE3eLNE+FCiEoYDAYiI8LdCQMgpVGSR539h4+yctUaoh68B4DQ+fOIfPrxGo1TiOpQYdJQSj3n+nNczYUjRP3QtFHZvo61m7ZzQKUBYI+KJu9vz9V0WEJcsMr6NMa5pkK/VynVqGSB1nqmb8MSom5LTkogKjKC3Lyzw22LQkL5z6Uj6dC6I90bRJPUOBlTyYMcDrDbwWQqcz4haovKbk/dAsThfCI8udSPEKISwUFBTL78Mob068Wtk8eQ2vTsf5vtTVvwSXgDvvxhEQWFZvf+0G++Im70UEzbt/kjZCG8YnA4yh1V66aU6g3sAVoD+7TWJ2sisPNVVGRzZGUFxsNUcXERBEpbS6trbXc4HGzWe1ixYQu5+QXu/Q3iYlCpzYnMz2XwbddiPH0KR3Aw2e9/gmVoxU+T17X2V6dAbjv4rv2JidFrgV5V1fOmI7wlsBJ4AliplLr+AmMTIuAYDAa6tG/DHyZfQb/uae79p7OyWbFhC0dmf4MjOxsAe1Ijivr081eoQlTKm6TxINBDaz0e6A7c79uQhKi/goOC6N+9M/17dPbYv61TN16993HyuvYg9+VpOKJj/BShEJXzJmnYtda5AFrrHKDQtyEJUf/16tSeiDDPGXEzGjXhP/c8hnnI8DL1wz75kOClS2oqPCEq5E3S2KOU+pdSapxS6l84+zeEEBcgODiIyy/pQ1hoCHHRUe79x06dZtf+dI+6pj27iHr0YeImXkHUw/dBXl5NhyuEmzdJ4w/AXmCY68/bfRqREAGiVbOm3H3dJG6dMoaeae3d+xcsX82pzDPu7chn/oqh0HmBH7RxA8iaHcKPqlyEybWmxowaiEWIgNWnS0e27d5HQaGZQrOFL374mYnDB9EooQG5L08Ho4mQH+eRM+0NCPJm7TQhfEOmERGiFggPC2XCsIEEuxJCfqGZL+cvIjc/H3ujxmS/9zGZv67A1rlLmWODNqxzPhgoRA2oMmkopSJK/imE8I3kxIZMunwQoSHOqdTNliIW/76B7xYt5a1PZ7MlKKzMMcErlxM/fBAx11+J8fChmg5ZBKDK5p76p1JqPPC9a9fcmglJiMDVtFEiVwwa4N7esfcAO/enk1dQyPwlKzl07MTZyoWFZydAXPAjkc88WdPhigBU2ZXGV0BjoLNSagHQXil1j1Kqb82EJkRgatG0Mc0aJ5XZb3c4+PbnJew7dIScvHyw2SgaNNhZJhMgihpSWdKIBD4CNmmthwH7gINAxXMbCCEumMFgYEj/3jSIK/uAX6HZwtc/LWb6O7M4ll9I7vNTyZzzI7nTXsfepKlnZYcDXGt6CFFdKpx7Sin1KNAbGAh8CVwKjNda76q58EApdTVwBXAKeNL1gGGFZO6pwBAobTdbijAY4NjJ08ya/wsl/7+ajEbGDrmE1KaNMRrLfv8L/fpLwv/9BjnT3sCW1rlMeV0VKO99RWrt3FNa6xe01pOAzcBLgAG4Tyk1v9qi9M5Y4GbgA+DGGn5tIfwqNCSYkOBgmic3YviAizzKbHY7sxcs5ucVa8scZzh5kqgnHiF443rihw8k5KcfaipkUc9VOOBbKRWhtc4HHtZa71NKXa+1XldOebVSSj0ADHVtrgBeB/4DpANyrS0CVlq7VqS1a8WG7bv4ecUa9/5NejdtW6SQmnJ2+vXgNb9jyM0FwN6oMUX9L67xeEX9VNlTQjOUUmuAzwCKE4ZSKgG4HufkhTdVd0Ba6+nA9OJtpdRE4DbgMqBZdb+eEHVN1/ZtiI2NYM7CZVitzu9RC5av5uaJo9zPeVhGjCJz0TKiH7yHvIf+jCMq2p8hi3qkwqShtb5FKXUl8I1SqhlwEogBjgJvuj7cz4lSqg/wotZ6kFLKCLwJdAXMwG1a693lHJYFvIfzKuPOc31NIeobg8FA945taRTfkHe//p5Cs4Xs3DzWbNlBv25np123tW1H1nc/Qom1y4uFffQ+tmbNKRp4WU2GLuqBKhdhAlBKhQHxwCmtteV8Xkgp9QhwA5Cnte7ruoIYq7W+2TWM9zGt9QWvR2632x02W2A8HWsyGbHZ7P4Owy8Cue1wtv2rN+1gzsJlAASZTNwwYTitmjep/OCdOwnq2R2D2Yz95luwTZsOkZE1EHX1kPfeN+0PDjZ51RHu1SQ2WutCnFcYF2IPMBH40LV9MTDfdf6VSqkqg/WGzeYImJEVgTyKJJDbDmfb3zolhcQGcZw4nYXVZuPD2T8xYfhAmic3qvDYmIcfxmB2LjNrW7eerAIbFNWdv0t57302esqrejU295TW+iugqMSuGOBMiW2bUkpmYhPiHBiNRsYMvpioiHAArDYbs39azPY9+zmTk8uydZvYuGMXmdlnR6rnTH0N85jxOIKCyJk+QyZAFOekyn8tSqmeWuuyY/ouXDZQMrUZXTPqCiHOQXxMNFeOHMIXP/xMbn4BVpuNeYtXlKnXLDmJ0YP6E9moEdn/+wDTTo2tnSpTL2jtaqw9epXbFyKEN1caf1ZKrXRNIRJXja+9DBgF4OrT2FyN5xYioMTHOhNHbIkFnUpLP5rByg1b3dvlJYzgFcuIHzmE2GsmYTyUXqZciCqThtb6amAk4AC+VEp9rJQaVA2vPRsoVEotB6bhXItcCHGe4mOjuXH8CLqo1h77w0JD3L9v37OfImsFF/QFBe4JEEMWLSTy73/1Wayi7vL2ZmYjoDmQAGwDpiilbtRa/+FcXkxrvR/o6/rdDvzfuRwvhKhcSHAwwwZcRI+OisMZJ2jdrCkR4WG8M2suWTm5mC1F6H0HSWvbquzBDgeWocMx/ectHFHR5D3zz5pvgKj1vOnTWAXk43wq+ymttdm1/0cfxyaEOE8N42NpGB/r3u6sWvPbmo0ALFu7iaiIcDbt2E3DuFj6dU/DYDBQaDTBsy9iHjcR47Fj2JNLDd0tngBROs4Dmjfv/j1a69XFG0qpgVrrxVprme1WiDqic7vWrN68nUKzhdz8Ar768VcAdh04RFCQiQOHj5F+LIO+XTsxoHefcs/hMQFiOSsIisBQ2dxTlwAdgQeVUq+4dpuAu4G0io4TQtQ+4WGhXH5xH779+bcyZUvXbnL/vnLjVnp36UBIcLBHneIJEI2nTxM/fCDZ73+CZfhIn8ctap/KOsIzcS7CFAoku34SgEdqIC4hRDVr0yKFYf17ExkeXmm9+b+t4sTpTI99wetWY8h3PlBmT24iEyAGsCqnEVFKJWutL/Rp8Boj62kEhkBuO1xY+x0OB3kFhRSazbw/u/wp04NMJq4bezkJJfpFTHt2EfXQfeTf/xBFg4ed12tXB3nva+l6GkqpWa5f1ymljrh+jiqljlRXkEKImmcwGIiKCCchPo6BvbuVW8dqs7F83SaPfbbWbTnzzbxyE0bYB+8SvGihT+IVtUtls9xOdv2ZXFEdIUTd1jOtPQCZ2TkcPHKcrJxcd9muA4c4fvI0jRIanD2gnKfETbt3EfXEIxjMZgqvupac56dCVMUPGYq6zZsht5cCETivSl4H/qq1/sTXgQkhfM9gMNCrcwcA8vIL2HXgEBt37OJkpnNauGXrNnFZ356cOJ1FSqNEIsLDypwj8u9/dU+AaNqxHcLK1hH1hzfTiLwE7ALuAwYgD+QJUS9FRoTTrUNbrhg0wL1v36GjvDNrLt8tWspn3y/0eJq8wGxm14FDnHrhX2SOuMI5AeK0N+Q5jnrOm3e3ADgOWLXWx5RSoT6OSQjhRw3jY+nQOpXte/Z77M/MzmGz3kOPTgq73c6XPyzixOksZ+GgMSSm9WFEcgpJpc4XtHoV1l4XyQSI9YQ3VxrZwELgC6XU3cBB34YkhPC3AT06l3sratWmbeTk5bNt9/6zCcPlREIS2/fu99gXvGIZ8aOHEXvleIwHD/gyZFFDvLnSuBJorbXeppRKwzmdiBCiHouNjuK2KWPIzS8gIiyU976eR25+AfkFhcz8/NsKjzt45PjZjZITIC7+hchnnyZn5ns+jlz4mjdXGknAFUqpp3CuvPeob0MSQtQGwUFBxMdEExoSwtD+vTFWcHupaaNE9+8ZpzI5knGSnfvTsVjMWIaPxGE0Yo+OkQkQ6wlvrjS+xHl7SibXFyJAtW7elPHDLuWnpb+Tm18AOKcmGdqvF+1aNufjOT9x7OQpAD6duwCA5MSGXPv3f2IeNwHjkSMyAWI94c27laO1ftLnkQgharWWKU24/cqx5OYX4HA4iImKxOC6+mjeJMmdNIodPXGKzOwc4nv2hp5lzxf61RdEzHiNnOlvYO3avSaaIKqBN7entiilrlZO7ZRS7XwelRCiVjIajcRERRIbHeVOGOCcRTckuOx30P2Hyp+ByHDiBFFP/oWgrZuJGzGYkPnzfBazqF7eXGl0c/0UcwCDfROOEKIuiouJ5o6rxpGbV8D+I0f5ddV6ADbs2EXHNqnuWXOLE03whrUeEyBaLr7UP4GLc1Zl0tBaX6aUigVaAHu11rlVHSOECDyhISGEhoRgMhndSeN0VjZvfPQVRqMRk9HI4L49CQ0NIaHPAKy/riD6T/eTf++DMu1IHeLNNCKTgCdddb9QSjm01s/6PDIhRJ0UFxNNfEw0mdk57n12ux273c6PS1e593Vu15qhX36L0WQqc46w99/BnpKCZcjwGolZeM+bPo2HcK7rfRJ4Fpjg04iEEHXekP69SGwQR0xUZIVDdTfv3MOGHbvL7Dft2knUE48Qe81kou++A3Ll5kZt4k2fhl1rbXZdYTiUUnk+j0oIUae1aNKYG8c7V/azWm0cP3WaL+cvwmaze9T7ZdU6tuzaS89OiuSkBGIiI4n5x1MYLBYATLs0hIVht9sxGAwene/CP7xJGr8ppT4FUpRSbwGrqzpACCGKBQWZaNookUnDB7FzfzqqZQvmLV5OTp6zI/zE6Szm/+a8bRVkMtFlys2MiYgk7LtvyJk2gx0HDzN/yUoS4uOYMOxS4uIi/NmcgFflyn0ASqkRQGdgh9b6O59HdQFk5b7AEMhth7rf/m279/HDkpUVlg/q052LwoPISU7hna/mYrYUAZDUMJ7/S21IXpdeATsBor9X7qv0SkMp1QWYjHNt8EPAjmqJTggR0Dq2aUlEeBins7I5knGSoydOkZ179s73r6vWsywoiCKr542NqN9XEnH7K+R170XB9BmEd+hQ06EHvMqWe50CvINz+pD5QA4wSyk1roZiE0LUY6lNk+nRSXHFZQO4/cqx3H/TlcRERbrLS67dARBssTDx648ASFy/hpwH7uNUVjbgXPfcm7sm4sJVNnrqfmCg1vo/Wus5WuvXgYuBB2omNCFEIAkymRjSr2e5nd3dO7YjpUkSWzt1xW4wUBAWzjcjJ/L1T79y/ORpPvx2Pm9/9i1HT5wq58yiOlV2e8qqtfYYKaW1zlFK2XwckxAiQLVq1pQ/TBpNgdlCVnYOJ05n0Ta1GcmJDcnMzuHT3AI2d+5F7JlMcmLiIDePj+b86DzY4WDxst+52jVqS/hGZUnDXsF+b57tEEKI8xIXE00czllyO7Q+uz8+JprbpowhLHwC67fsYetyz/6ObhtWMXDxArJiQokb7JzpqNBswW63l7uglDg/lSWNTkqpT0rtMwAdfRiPEEJUKCQ4mJioCDq3a8XarTvIPON86jwyN5sx331BZH4e9msmcuo/73O8/yXMmv8LFquVUQP7oVo293P09UNlSePKCva/5YtAhBDCW0ajkTGXXczStZtIiI+hy96dBBc5Hwg8ExvPW8ezySu+bQX8tHQVjRMaEBstc1xdqAqThtZ6cU0GIoQQ5yKxgfNhPwB6dWNFfEOaPP0YSy69nDyT50ebpcjKopVrmTBsoB8irV9kySwhRL3QcdgQdrWZxb4lK50rApbQZ+VizmyP55NCM21TmxETFcn2PfvJys5FtWxOzzTlnr4dICs7h/RjGbRunkJEWGhNN6VWk6QhhKg32rZsztXRUazevJ2Q4GBaNG3M8aVLGfPdFwTZrKzfuJrZE67DEnq2Y3z5+s3sP3yUq0YNwWg0Yikq4vN5P5ObX0Dz5ANMGSnLB5UkSUMIUa80SmjAFZcNcG93efwrgmzOBwUTTmZgDQouc8yRjJNMe+9zoiMjMBgM7nXQDx49TkGhmXC52nCrlcNnlVKDlVL/df3eXyn1vusnzt+xCSHqFsuMmRwYPByrycT3191BXMN4khrG071jOxonNPSom5OX7zGdCcDh4ydqMtxar9ZdaSil2gA9gOLrxzuAO4GLgKuAt/0UmhCiDnI0aEjEZ7PI2reXSS1beZRZior46YWp6CYtKpwA8dDxE7RpkVITodYJfk8aSqkHgKGuzRVa6+eAqUqpj1z7TFrrQqXUUWRtciHEeXKUShgAkatWcMvrz3O4Q2e2//lJTjZMZMP2XR51Dh/LwOFw8MuqdZzKOsOQfr1oEBtTU2HXOn5PGlrr6cD0SqrkK6VCgWTgWM1EJYSo9/LziX7oXgCabt9Mwzmfk/32uzRJSuDw8ZNs3OFMHsdPZTJv8Qp27D0AwOLf1wf00F2fJg2lVB/gRa31IKWUEXgT6AqYgdu01mXXeixrJs5bUsE4b1MJIcSFMxoxjxlP+Juv4YiMIu/vz2MwGOjQOpUOrVM5nXWGdNdVRnHCANibfgRLUZHHEN1A4tUiTOdDKfUIcAOQp7Xuq5SaCIzVWt+slOoLPKa1rvZp1u12u8NmC4wpkk0mY5nlMwNFILcdArv91d72dWsxHDqMY+xYj92ZWdm8//lcTuUVlDnkqisGk9auZfXFcA589d4HB5sufBGmC7QHmAh86Nq+GOe6HGitVyqlqgzufNhsjjq9otm5qOurt12IQG47BHb7q73trTo4f0qdM+yLL7nvtVf477AJpDf3TBCfz13E8uQkGjVswIAeXQgKMlVfPFXw4cp9XtXz2ZBbrfVXQFGJXTHAmRLbNqWU3/tUhBCiNENGBlF/fZTQXZo//vtFOm7dUKZO+tEM1mzZwVufzeanpb+TV1D2iqQ+qskP7WygZCozaq2tFVUWQgh/Cd68AYPZDIC9WXM63HEb7SMj2bH3ANt27/eoa7YUsXnnHhw4uPziPn6ItmbVZNJYBowBvnD1aWyuwdcWQgivWYYM5/TilUT/6X7y776fVNUWcC5R2zihIVt37+P4ydMex+w9eNhjydkdew+Ql19AWrvWhIWG1Gj8vlSTSWM2MEwptRznuhy31OBrCyHEObG3SOXMF994PPRnMBjo3rEd/Vb8ir1pCqcuGcTMz78FIL/QzLK1m1i/fReWorN35jfu2M3EywcRH+Ndn0Gx02eyWbh8NfEx0Qzp1wujsXZM4OGz0VP+UlRkcwRKB6F0hgZm2yGw2+/vtpv0DuKHXIzBYqFw/ERmTbieLUcrn2qkZUoyE4cPOqfXmfvLMvS+gwCMuKQPndo6H070YUe4V6OnakfqEkKIOiLyuWcwWJwLPpkOHiClZYsqj0k/loHdfm7DZIsTBsC6bTvPLUgfkqQhhBDnIOfVGRRedS2O4GByps2gZYtmBAd53ukPDgpi0EXd3bPjWq02MrNzy5wrv6CQTXoPJzPPeOwvfQcov6CQLTv3Mmv+L+w5cLiaW3RuJGkIIcQ5cMQ3IOf1tzi9fC22Dh2JCA9j9KD+7sTRIX0v/3fVWHqmtfeYRfe9r79n9oIl2FwLRG3fs593vprLgmW/8+ncBe7p2AEKzRaP18zNL+DHpas4cOQY3y5cVgOtrJgkDSGEOA/2Fqnu31s3b8oN40YwzljETTNeJPHK8Zj27iapgedqDnvTD7Ntz34OHcvghyUrMVucHeaWoiL3XFcAWTllr0qKZZ7JKXMlUpPk4TohhKgG8UFGWs+YCkDIimVEvPgcSY8+U6ber6vWExRkKvPBv3HHblo0SSYrO8c5vrQSVquN4GD/fHxL0hBCiOpgMlE4YTIRr0/DERVN3t+fJzE8skw1S1GRe0huaEgwliIrDoeDgkIzn89b6NVLFVoskjSEEKJOCwsj//GnsIwZh/HwYeyNGhPncBAdGUFOXj44HJhsNmwlOs1HD+rPidNZ/LZm4zm9VKHZTHRkRHW3wCvSpyGEENXI2rkrlhGjAOfDgOOHXkqvzu35oz2XR956keYH9gDQt2snWqY0obNqXen5So/MAig0F5VTs2ZI0hBCCB9KahjPoBZNSXn5WWIPHeSut17mjlAbA3p2ASA8NJQuqk2Fx19x2QCMpZaiLXTNi+UPkjSEEMLHgrZsgiLn/Kz2Zs2JHj/Bo7x/9zSSGsYTFRHusb935w60ataEK0cN8dhfaPEckluTpE9DCCF8rGjwUDKXrCT6zw+Qf9e9EOnZQR4ZEc4N40YAcDIziz0Hj9CqWTKJDeIBaNookZ6dFGu3aqDscxw1SZKGEELUAHuz5pz57Otyy8L++xb25KZYRo8hIT6OhPi4MnVCS8yUK0lDCCEClEnvIOrpJzAUFWEeM56c6W/giI4pUy8spHYkDenTEEIIP4r8598xuJ7bMB45hCOi7LMdgMeaHP7s05CkIYQQfpTz6gwKrr0BR0gIOdNmgKn89cY9koaMnhJCiMDkiIsnd/oMTq9Yh619hzLlwb8tBrudsNBQ9z65PSWEEAHO3qx5mX3BS34lbtIY4saNJPrQ2fU1zJI0hBBCeMjLI/qh+wAIXrWCpBnT3UXSpyGEEMJTUBCFk6/EERSEPTaO/GdfcBeZLUXnvBJgdZGkIYQQtVFoKPmPPknmgiXkzHgbGicTGhLsLl66ZiPWvJpfK12ShhBC1GK2TmlYho8EoCDdQ9IAAAfPSURBVGFcrHu//b13COnTnawffqjReCRpCCFEHTF6UH9aNG1EVM4Zrpj7JQ0zjtL2pqvQr0xzr9Hha5I0hBCijoiJiuQPV45mZHwEBsf/t3f/wVaUdRzH34fhpxgRoULDoE7IJ2rKsgRKGW9loMMURc1IpsQ0UzpYyVhj1tTkj3QGa/ojTc1qgiYU8EdFNYgzEBohaWo2xvhtBoIK5ceVlDR+qNAfz3OndblX9nDPueee+Lz+Obv77J59du9z9rvP7t3vpnsaz44azarj3sDin6/k6Z2dTa+Dg4aZWRsZUKsxfu5ctq9ex9bT38O9sy/ipcFD2PPCi6xYva7p7w937ikzszY0fOJEDt2/mombtrDtD4+xb/8Bhg4ZRK12hBeM95KDhplZm6rVarx1wqmcOu5NbNuxi3FjTmz6Oh00zMza3LChQ5hw8rg+WZfvaZiZWWUOGmZmVpmDhpmZVeagYWZmlTlomJlZZQ4aZmZWmYOGmZlVVmv2I+ctsAvY2upKmJm1mZOBE4400/9j0DAzsybx5SkzM6vMQcPMzCpz0DAzs8ocNMzMrDIHDTMzq8xBw8zMKvP7NPohSVOAhRHRIWkCsAg4BDwJXAaMB5ZGxNQ8/9nAT4CPA1fksvtaUfejdaRtjoiDkr4JzAReBhZExMOSFpG3V9JA4A6gE/g2cGfXPmo3pf2xCPgocFJE7M/lZwCPAu8HTgHeEhFXtai6DVGx3f8NuCoiFhaWWwGMyMttIe2LfX1c/aPSy3bfkjbhnkY/I+lK4EfA0Dzpu8DXI2IaUANmlebvAG4HZkbEn/uwqg1TZZvzD+IcYAowB/h+6TsGAcuBzRExn/TDa0vd7A+AZ4DzC+OfAjb3Zb2aqY52vwn4RGG5UcBpfVjVhmlAu29Jm3DQ6H82AbML4+8GHsjDK4FzuwoknQvcApwXEX8tLDNf0mpJD+Szl/6uyjafDdwfEYci4u/AQEldT68OAe4F/lQ6szpB0gpJGyR9o7mb0FDl/QFwJ/BJAEkDgDOARwrl781/80ckzeybajZU1XbfCeyUNCmPXwDcVfquH+S2v1zSsGZVuAF62+5b0iYcNPqZiLgHeKkwqRYRXWfN/wZen4ffDFxPOksp/zDWR8QHgYXAjU2sbkNU3OYRwPOFeYr74nvAcKD8vsvjgYuBs4DzJZ3e4Ko3RTf7A+BhQJKGAx8Aflsqf5F0kJkJ3JwPIm2jjnYP6WA5Jw/PAn5R+rpbI+IcYAvw2cbXtjEa0O5b0ibaqmEdow4Whl8HPJeH95K6pl8EymdUD+bP9YCaXsPG626b9+Th8nSAm4DpwNslXVSY54mIeD4iXiH9wCY2r8p9YgXpIHkhsKRUti6fje4kHWTe2NeVa7Ce2j2kIPERSacA24H/FMoORMSGPNxu7b/edg8taBMOGv3f4/m+BaQg8bs8/HRE7I6IX+dpNxeWmZw/p5FuqLWb7rb598AMSQMkjQcGRERnnufJiHiZdE33xsKli0mSjs83yKcAf+m7TWiKJcBcYGxEbCqVnQkgaQyph9VJe+up3RMRLwBB6kXfUVpusKR35uF2a//1tntoQZtw0Oj/vgRcI+khYDBwdzfzfBmYLOnTeXyqpDXAAuDKvqlmQx22zRHxKOlH9BBwD+m/aV4lIjaTtvcu4DhgN7CMdMZ5d0Rs7JvqN0dEBCkL6a+6KR6W/+YrgEsKlzna1ZHa/RJSUFhdmr4f+IKkB4ETSTea20Xd7b4VbcJZbs3MrDL3NMzMrDIHDTMzq8xBw8zMKnPQMDOzyhw0zMysMgcNOyZJ6pC0tJvpSyUNfo3ltte5njGSbjmaOha+Y5qky3uqcw/LDJO0WFKtN+s2K3OWW7OCiJhz5Lnq+r7twPyjXT4f9K8mPez1vjrWu1fSetKDX4uPdv1mZQ4aZgVdqbWB20gPip0CjAXmRcRjhfluIOUA+nzXQ1M5kdwyUg9+EHApKVfQUlK66pV58YHAVFKKi3GkHGKvkBLYXRIRxXxEHwI2RsQBSQCnSVpFSglxa0T8WNJa4Klc7xpwQQ5Wy4H7cNCwBvLlKbOebY2IGaTcVp/rmijpO8DAiLis9JTtZFKOn66cYCO6CiJib0R0kILHFlLvYzPwQ2B2TrC3DZhXqkMHUEx5Pwj4MOlp6K8UMp6uz9+/DPhaXue/gNGSisn+zHrFQcOsZ4/nz3/wv3cenAS8g5TLp2wlKbX1L4FreXUCui43kbI/3E5K/zCWlHByLSnp4vjS/KOBHYXxDRFxICL2AhtJPSGANfmznKRvBzCqxy00q5ODhlnPusuxswOYAbxN0nmlsg7gmYiYDnwLuKFYKOk6Uvrr6/KkTuCfwKzcS7iew9Nb7wRGFsbfJWlgToc9iXRJC9K7GCClgS8mZhwJ7HqNbTSri+9p2LFsuqQ/FsYvrLJQRByS9BlglaQpEfFsLnoCWCZpAekexbVdy0g6E/gqsDb3KsjllwO/ye862EO6cV20FvgY8NM8vo/UoxkJXB0Ru/O9jnmSriC9Q+HivM6RwHM5K6xZQzhhoVk/loPJGmB6RBzoYZ61wKUR8VRp+nxgT0T8rOkVtWOGL0+Z9WMRcRC4hjr/bTe/lOssDn/fhFmvuKdhZmaVuadhZmaVOWiYmVllDhpmZlaZg4aZmVXmoGFmZpX9FxlFSaM0Wxg4AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "A, B = 0.000325, -0.8844\n",
    "title = \"Molokai Chr01A (RE=AAGCTT)\"\n",
    "label = r\"$Y=0.000325 \\times X^{-0.8844}$\"\n",
    "filename = \"Molokai_RE_AAGCTT.Chr01A.distribution.txt\"\n",
    "plot_distribution(A, B, title, label, filename)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python [default]",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.15"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
